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PEEFAGE. 


THE  object  of  the  following  treatise  is  to  exhibit,  in  a  concise 
form,  the  elementary  properties  of  the  expressions  known  by 
the  name  of  Laplace's  Functions,  or  Spherical  Harmonics, 
and  'of  some  other  expressions  of  a  similar  nature.  I  do  not, 
of  course,  profess  to  have  produced  a  complete  treatise  on 
these  functions,  but  merely  to  have  given  such  an  introduc 
tory  sketch  as  may  facilitate  the  study  of  the  numerous 
works  and  memoirs  in  which  they  are  employed.  As 
Spherical  Harmonics  derive  their  chief  interest  and  utility 
from  their  physical  applications,  I  have  endeavoured  from 
the  outset  to  keep  these  applications  in  view. 

I  must  express  my  acknowledgments  to  the  Rev.  C.  H. 
Prior,  Fellow  of  Pembroke  College,  for  his  kind  revision  of 
the  proof-sheets  as  they  passed  through  the  press. 


N.  M.  FERRERS, 
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August,  1877. 
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CHAPTER  I. 

INTRODUCTORY.      DEFINITION  OF  SPHERICAL  HARMONICS. 

1.  IF  V  be  the  potential  of  an  attracting  mass,  at  any 
point  x,  ?/,  z,  not  forming  a  part  of  the  mass  itself,  it  is 
known  that  Fmust  satisfy  the  differential  equation 


da?      dy        dz 
or,  as  we  shall  write  it  for  shortness,  V2F  =  0. 

The  general  solution  of  this  equation  cannot  be  obtained 
in  finite  terms.  We  can,  however,  determine  an  expression 
which  we  shall  call  Vit  an  homogeneous  function  of  x,  y,  z 
of  the  degree  i,  i  being  any  positive  integer,  which  will 
satisfy  the  equation  ;  and  we  may  prove  that  to  every  such 
solution  Vt  there  corresponds  another,  of  the  degree  -  (i  +  1), 

expressed  by  -^  ,  where  r*  =  x*  +  yz  +  z*. 

For  the  equation  (1)  when  transformed  to  polar  co-ordi 
nates  by  writing  x=-r  sin  6  cos  <t>,y=r  sin  0  sin  <£,  z  =  r  cos  0, 

becomes 


d"(rV)        1     *  f.gdr\        I  ™ 

r  ^!?—  +  si^  d9  (SU  6Wj  +  ffi*fl  dp  ~ 

And  since  V  satisfies  this  equation,  and  is  an  homo 
geneous  function  of  the  degree  /,  F£  must  satisfy  the  equa 
tion 


F.  H. 
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since  this  is  the  form  which  equation  (2)  assumes  when 
is  an  homogeneous  function  of  the  degree  i. 

Now,  put  Vi  =  r2*1  Uit  and  this  becomes 


or 


Now,  since   U{  is  a  homogeneous  function  of  the  degree 


therefore  equation  (2)  becomes 

_    /  .    QdU\         1     dzU, 
+  - — ^  -TX  f  sin  0  -~~    +  -T-J-V  - ,    *  =  0, 


shewing  that   t/i  is  an  admissible  value  of  V,  as  satisfying 
equation  (2). 

^  It  appears  therefore  that  every  form  of  U.  can  be  ob 
tained  from  Vit  by  dividing  by  rim,  and  conversely,  that 
every  form  of  Vi  can  be  obtained  from  Ut  by  multiplying 
by  r2m.  Such  an  expression  as  Vi  we  shall  call  a  Solid 
Spherical  Harmonic  of  the  degree  i.  The  result  obtained 
by  dividing  Vt  by  ?•',  which  will  be  a  function  of  two  inde 
pendent  variables  6  and  $  only,  we  shall  call  a  Surface 
Spherical  Harmonic  of  the  same  degree.  A  very  important 
form  of  spherical  harmonics  is  that  which  is  independent 
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of  <f>.  The  solid  harmonics  of  this  form  will  involve  two  of 
the  variables,  x  and  y,  only  in  the  form  a;2  +  2/2,  or  will  be 
functions  of  x*  +  y*  and  z.  Harmonics  independent  of  <f>  are 
called  Zonal  Harmonics,  and  are  distinguished,  like  spherical 
harmonics  in  general,  into  Solid  and  Surface  Harmonics. 
The  investigation  of  their  properties  will  be  the  subject  of 
the  following  chapter. 

The  name  of  Spherical  Harmonics  was  first  applied  to 
these  functions  by  Sir  W.  Thomson  and  Professor  Tait,  in 
their  Treatise  on  Natural  Philosophy.  The  name  "  Laplace's 
Coefficients"  was  employed  by  Whewell,  on  account  of 
Laplace  having  discussed  their  properties,  and  employed 
them  largely  in  the  Mecanique  Celeste.  Pratt,  in  his 
Treatise  on  the  Figure  of  the  Earth,  limits  the  name  of 
Laplace's  Coefficients  to  Zonal  Harmonics,  and  designates 
all  other  spherical  harmonics  by  the  name  of  Laplace's 
Functions.  The  Zonal  Harmonic  in  the  case  which  we  shall 
consider  in  the  following  chapter,  i.e.,  in  which  the  system 
is  symmetrical  about  the  line  from  which  6  is  measured, 
was  really,  however,  first  introduced  by  Legendre,  although 
the  properties  of  spherical  harmonics  in  general  were  first 
discussed  by  Laplace;  and  Mr  Todhunter,  in  his  Treatise, 
on  this  account  calls  them  by  the  name  of  "Legendre's 
Coefficients,"  applying  the  name  of  "Laplace's  Coefficients" 
to  the  form  which  the  Zonal  Harmonic  assumes  when  in 
place  of  cos  6,  we  write  cos  0  cos  &  +  sin  0  sin  6'  cos  (</>  —  <£>'). 
The  name  "  Kugelfunctionen "  is  employed  by  Heine, 
in  his  standard  treatise  on  these  functions,  to  designate 
Spherical  Harmonics  in  general. 


1—2 


CHAPTER  II. 

ZOXAL   HARMONICS. 


1.  WE  shall  in  this  chapter  regard  a*  Zonal  Solid  Har 
monic,  of  the  degree  i,  as  a  homogeneous  function  of 
(x*  +y2)*,  and  z,  of  the  degree  ^*  which  satisfies  the  equation 

d?V    <FV    d*V_ 
<Lx?  +  dy'+'d?' 

Now,  if  this  be  transformed  to  polar  co-ordinates,  by 
writing  r  sin  0  cos  </>  for  x,  r  sin  9  sin  <$>  for  y,  r  cos  6  for  z,  we 
observe,  in  the  first  place,  that  x2  4-  yz  =  r2  sin2  0.  Hence 
V  will  be  independent  of  <£,  or  will  be  a  function  of  r 
and  0  only.  The  differential  equation  between  r  and  6 
which  it  must  therefore  satisfy  will  be 

d*(rV)  ,      1      d    f  .    adV\ 
r  —  ra~  +  -T—  >.  -70    sm  6  -T-    =  0. 
dr          sin  0  dO  \          dO  J 

Now  F",  being  a  function  of  r  of  the  f  degree  i,  may  be 
expressed  in  the  form  ?**P.,  where  Pt  is  a  function  of  ^  only. 
Hence  this  equation  becomes 


or,  putting  cos  6  =  //-, 


In  accordance  with  our  definition  of  spherical   surface 
harmonics,  Pi  will  be  the  zonal  surface   harmonic   of  the 

YC<U    &t  f>.  un  , 
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degree  i.     When  it  is  necessary  to  particularise  the  variable 
involved  in  it,  we  shall  write  it  Pt  (p). 

The  line  from  which  0  is  measured,  or  in  other  words 
for  which  ft  =  1,  is  called  the  Axis  of  the  system  of  Zonal 
Harmonics  ;  and  the  point  in  which  the  positive  direction 
of  the  axis  meets  a  sphere  whose  centre  is  the  origin  of 
co-ordinates,  and  radius  unity,  is  called  the  Pole  of  the 
system. 

Any  constant  multiple  of  a  zonal  harmonic  (solid  or 
surface]  is  itself  a  zonal  harmonic  of  the  same  order. 

2.  The  zonal  harmonic  of  the  degree  i,  of  which  the 
line  p  —  1  is  the  axis,  is  a  perfectly  determinate  function  of 
JJL,  having  nothing  arbitrary  but  this  constant.  For  the 
expression  r*Pt  may  be  expressed  as  a  rational  integral 
homogeneous  function  of  r  and  z,  and.  therefore  P.  will  be 
a  rational  integral  function  of  cos  6,  that  is  of  p,  of  the 
degree  i,  and  will  involve  none  but  positive  .  integral  powers 

Of  yLt. 

But  Pt  is  a  particular  integral  of  the  equation 


and  the  most  general  form  of  f(p)  must  involve  two  ar 
bitrary  constants.      {Suppose  then   that   the   most  general 

form  of  f(p)  is  represented  by  P{  \vdp.     We  then  have 


_ 

}      dp  j;          '  dp 


Hence,  adding  these   two  equations   together,    and    ob 
serving  that,  since  Pi  satisfies  the  equation  (3),  the  coefficient 
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of  |w7yu,  will  be  identically  equal  to  0,  we  obtain,  for  the  de 
termination  of  v,  the  equation 


whence      P,  (1-^+2   (l  -^2)        -/^j*  =  0, 

cZ»      /2  dP,        2/A   \  7 
or  _  +  f  _  __c_    A*  =  0, 


the  integral  of  which  is 

log  v  +  log  P,2  (1  -  yu,2)  =  log  Cl  =  a  constant  ; 


Hence          /«fr  =  C+  GI  j  ^^  ; 
and  we  obtain,  for  the  most  general  form  o 


Now,  P{  being  a  rational  integral  function  of  //,  of  i 
dimensions,  it  may  be  seen  that    r  -  ^—  ^  will  assume  the 

»  \      —  f*  )  *   i 

form  of  the  sum  of  i  +  2  logarithms   and   i  fractions,  and 
therefore  cannot  be  expressed  as  a  rational  integral  function 

of  fju.    Expressions  of  the  form  Pi  I          ^         are  called  Kugel- 

functionen  der  zweiter  Art  by  Heine,  who  has  investigated 
their  properties  at  great  length.  They  have,  as  will  hereafter 
be  seen,  interesting  applications  to  the  attraction  of  a  sphe 
roid  on  an  external  point.  We  shall  discuss  their  properties 
more  fully  hereafter. 

3.     We  have  thus  shewn  that  the  most  general  solution 
of  equation  (2)  of  the  form  of  a  rational  integral  function  of  a 
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involves  but  one  arbitrary  constant,  and  that  as  a  factor. 
We  shall  henceforth  denote  by  Pt,  or  Pt  (p),  that  particular 
form  of  the  integral  which  assumes  the  value  unity  when  p 
is  put  equal  to  unity. 

'  We  shall  next  prove  the  following  important  proposition. 

I  fb.be  less  than  unity,  and  if  (I  -  2yu,h  +  h1)""  be  expanded 
in  a  series  proceeding  by  ascending  powers  of  h,  the  coefficient 
of  h1  will  be  P,. 


Or,  (l 

We  shall  prove  this  by  shewing  that,  if  //  be  written  for 
(1  —  2/zA  +  h*)~^,  H  will  satisfy  the  differential  equation 


For,  since  H=  (1  -  2/^  +  A,8)'*, 


1  dff_ 

-"•  H9  d,      ' 


And  H, 


d 

•*'  dh 
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„  Id  f  ,,        „.  d/f  ,,  „. 

Hence  f  -y-  -U  1  -  /r)  -7-  [  +  ^  (MT) 

"  *  v      y 


dh* 

=  -  S/,//3  +  3  {(1  -  ^2)  7t  +  (1  -  M)  0  -  70] 
=  -  3/z//3  +  3  {^  (1  +  A2)  -  2  A}  #5 

=  -  3/*  ;/^3  -  (i  -  2fjLh  +  A2)  /y5} 

=  0,  since  1  -  2yu/i  +  h*  =  2I^. 
Therefore          |{d-  ^)  f 

This  may  also  be  shewn  as  follows. 

If  #,  ?/,  «  be  the  co-ordinates  of  any  point,  /  the  distance 
of  a  fixed  point,  situated  on  the  axis  of  z,  from  the  origin, 
and  JR  be  the  distance  between  these  points,  we  know  that, 


and  that 


Now,  transform  these  expressions  to  polar  co-ordinates, 
by  writing 

x  =  r  sin  0  cos  0,     y=-rs'm0  sin  0,     z  =  r  cos  0, 
and  we    et 


and  the  differential  equation  becomes 

1      d  ( .    „  d  fl 
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or,  putting  cos  6  =  p, 


r-r-« 


2N   d 


Now,  putting  r  =  zh,  we  see  that 


1      H 

£  =  7' 


and  r  1 

d 

.'.  the  above  equation  becomes 
h  d2  ,,  ,-,.x       d 


or 

4.     Having  established  this  proposition,  we  may  proceed 
as  follows: 

If  p.  be  the  coefficient  of  h*  in  the  expansion  of  H, 
/.  liH—h 

Also,  the  coefficient  of  h*  in  the  expansion  of 
d 


is  T"  1 U -*l" 

rfyU,   (^  ^/ 

Hence  equating  to  zero  the  coefficient  of  A4, 
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Also  ;;.  is  a  rational  integral  function  of  p. 
And,  when  p  =  1,  11=  (I  -  2h  +  tfy* 

=  1+A  +  A2  +  .  ..+/*'•  +  ... 
Or  when  p  =  1,  pt  =  1. 

Therefore  pt  is  what  we  have  already  denoted  by  Pit 
We  have  thus  shewn  that,  if  h  be  less  than  1, 


If  li  be  greater  than  1,  this  series  becomes  divergent. 
But  we  may  write 


since  ,  is  less  than  1, 


P      P 


Hence  Pi  is  also  the  coefficient  of  h~(i+V  in  the  expan 
sion  of  (1  —  2/j,h  +  h*)  ~*  in  ascending  powers  of  -.-  when  h  is 

ii 


greater  than  1.     We  may  express  this  in  a  "notation"  which  is 
strictly  continuous,  by  saying  that 

P,-P-««). 

This  might  have  been  anticipated,  from  the  fact  that  the 
fundamental  differential  equation  for  P,  is  unaltered  if 
—  (i  +  1)  be  written  in  place  of  i ;  for  the  only  way  in 
which  i  appears  in  that  equation  is  in  the  coefficient  of  Pf, 
which  is  i  (i  -f  1).  Writing  -  (i  +  1)  in  place  of  i,  this  be 
comes  -(i  +  l)  [-  (i +  !)+.!}  or  (i  +  l)i,  and  is  therefore 
unaltered. 
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5.     We  shall  next  prove  that 


where  r*  =  a?  +  y*  +  z*. 

Let  lr 

and  let  k  be  any  quantity  less  than  r. 

Then  {*•*  +  f  +  (z  -  &)'}4  =/  (*  -  Q, 

and,  developing  by  Taylor's  Theorem,  the  coefficient  of  k'  is 


Also       ai8  +    2  +  ^  -  7;~  =  (f  ~  2&*  +  ^T 


since  2  =  /ir, 

in  the  expansion  of  which,  the  coefficient  of  k*  is 

p, 
?•• 

Equating  these  results,  we  get 


The  value  of  P,  might  be  calculated,  either  by  expanding 
(1  -  2u/i  4-  &2)~*  by  the  Binomial  Theorem,  or  by  effecting  the 

rm         d*  /IN 
differentiations  in  the  expression  (—  l^f—  5—  «  --  --71  (-  )> 

1  .  Z  .  o  .  .  .  1  \LZ  \i  J 

and  in  the  result  putting  -  =  p.     Both  these  methods  how 

ever  would  be  somewhat  laborious  ;  we  proceed  therefore  to 
investigate  more  convenient  expressions. 
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6.     The   first   process   shews,  by  the  aid  of  Lagrange's 
Theorem,  that 


Let  y  denote  a  quantity,  such  that 


h  being  less  than  1. 
Then 


(i  -  2M  +  /O 
A 

\i  2 

Also 


Hence,  by  Lagrange's  Theorem, 


f  i1  V 

2    J  +"" 

therefore,   differentiating  with  respect  to  /-t  and  observing 
that 
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1  ^7* 

Hence          P'  =  2M.2. 

7.  From  this  form  of  Pi  it  may  be  readily  shewn  that 
the  values  of  p,  which  satisfy  the  equation  Pt  =  0,  are  all  real, 
and  all  lie  between  —  1  and  1. 

For  the  equation 

(JL*  —  1)'  =  0  has  i  roots  =  1,  and  i  roots  =  —  1, 

...   lL  fa*  _  l)<  =  0  has  i  -  1  roots  =  1,  (i  -  1)  roots  =  - 1,  and 

d/jju 
one  root  =  0, 

**LQM»-.1)«SBO  has  (t-2)  roots  =  1,   one   root   between   1 

and  0,   one   between  0  and  =-1,   and   (i-2)    roots  =  -1, 
and  so  on.     Hence  it  follows  that 

— .  (u?  —  1)*  =  0  has  5  roots  between  1  and  0,  and  -  roots  be- 
a/A  ^2  & 

tween  0  and  —  1,  if  i  be  even, 

and  ^—  roots  between  1  and  0,  — ^—  roots   between  0  and 
1,  and  one  root  =  0,  if  i  be  odd. 

It  is  hardly  necessary  to  observe  that  the  positive  roots  of 
each  of  these  equations  are  severally  equal  in  absolute  mag 
nitude  to  the  negative  roots. 

8.  We  may  take  this  opportunity  of  introducing  an  im 
portant  theorem,   due  to  Kodrigues,  properly  belonging  to 
the  Differential  Calculus,  but  which  is  of  great  use  in  this 
subject. 

The  theorem  in  question  is  as  follows: 
Jfrfile  any  integer  less  than  i, 
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It  may  be  proved  in  the  following  manner. 

If  (x~-  I)1  be  differentiated  i  -m  times,  then,  since  the 
equation 

(^_iy  =  Q 

has  i  roots  each  equal  to  1,  and  i  roots  each  equal  =  -  1   it 
follows  that  the  equation 


has  i-(i-m]  roots   (i.e.  m)  roots  each  =1,  and  m  roots 
each  =  —  1,  in  other  words  that  (x*  —  ljw  is  a  factor  of 

dl~m 

(3*  -I)'- 


We  proceed  to  calculate  the  other  factor. 
For  this  purpose  consider  the  expression 


Conceive  this  differentiated  (I)  i  -  m  times,  (II)  i  4  m 
times.  The  two  expressions  thus  obtained  will  consist  of  an 
equal  number  of  terms,  and  to  any  term  in  (I)  will  corre 
spond  one  term  in  (II),  such  that  their  product  will  be 

term  in  (II)  is  the  product  of  all  the  factors  omitted  from 
the  corresponding  term  in  (I)  and  of  those  factors  only. 
Two  such  terms  may  be  said  to  be  complementary  to  each 
other. 

Now,  conceive  a  term  in  (II)  the  product  of  p  factors  of 
the  form  as  4  «,  say  x  +  af,  x  4  a"  ...  x  4  o.(p\  and  of  q  factors 
of  the  form  x  4  &  say  x  4  £,,  x  4  $„ . . .  x  4  /3[y}.  We  must 
have  p-\-  q  =  i  —  m. 

The  complementary  term  in  (I)  will  involve 
p  factors  x  4  ft ,  x  4  @"  . . .  x  4  fi(p\ 
q  factors  x  4  «,,  x  4  a,,  . . .  x  4  «(7). 
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Now,  every  term  in  (I)  is  of  i  +  m  dimensions.  We  have 
accounted  for  p  -f  q  (or  i  —  m)  factors  in  the  particular  term 
\ve  are  considering.  There  remain  therefore  2m  factors  to 
be  accounted  for.  None  of  the  letters 

a',    a"  ...a"",  /?„    £„.  ..#,,, 
IB',   /?"...  /5<*>,  «„    «„...«„„ 

can  appear  there.     Hence  the  remaining  factor  must  involve 
m  a's  and  m  /3's,  —  say, 


There  will  be  another  term  in  (II)  containing 

(x  +  P)  (x  +  /3")  ...(»+  £<")  (a  +  a)  (x  +  aj  .  .  . 

The  corresponding  term  in  (I)  will  be,  as  shewn  above, 

(x  +  a')  (a  +  a")  ...(«+  ^')  (a  +  ft)  (a  +  ft,)  .  .  .  (x  +  fta)) 

(«  +  ,«)  (aj  +  ^)  ...  (x  +  ma)  (a  +  1/3)  (^  +  2^)  ...  (a  +  JS). 

Hence,  the  sum  of  these  two  terms  of  (I)  divided  by  the 
sum  of  the  complementary  two  terms  of  (II)  is 

)  (a  +  ^  (a;  +  2/9)  ...  (^  +  m/9). 


Now,  let  each  of  the  a's  be  equal  to  1,  and  each  of  the  fit's 
equal  to  —  1,  then  this  becomes  (a?  —  l)m.  The  same  factor 
enters  into  every  such  pair  of  the  terms  of  (I).  Hence 


c-       ,  2    ^mx- 

"  ~^m     =  ^   "  ^  —  s^^  —  '  to  a  numencal 

factor  pres. 

The  factor  may  easily  be  calculated,  by  considering  that 
the  coefficient  of  xi+m  in  d          ~ 


_ 

(IX 

Ji+m  /    2  _   -I  \i 

and  that  the  coefficient  of  x^  in  -  —  {'  i+ltl     }  is 

U.X 
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Hence  the  factor  is 

(i+  m)  (i  +  m  —  1)  . . .  (i  —  m  + 1)  '       1.2...  (i  +  m)  * 

d* 

9.    This  theorem  affords  a  direct  proof  that  C -y-.  (V  —  1)', 

VIM 

C  being  any  constant,  is  a  value  of  f  (p)  which  satisfies  the 
equation 


For  ^-i).>'-i}'=^- 

=»(i+ 

from  above, 

d      f/     2          IN     ^      f^*      /     2  -I  Nil!  •/« 

••  ^  r  1}  ^  te  (/t  -  1}  I  J  ~  *  (* 


or 


Hence,  the  given  differential  equation  is  satisfied  by  put- 


Introducing  the  condition  that  Pi  is  that  value  of  /  (JJL) 
which  is  equal  to  1,  when  /*  =  1,  we  get 

1  d* 


10.  We  shall  now  establish  two  very  important  proper 
ties  of  the  function  Pi  ;  and  apply  them  to  obtain  the  develop 
ment  of  Pi  in  a  series. 
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The  properties  in  question  are  as  follows  : 
If  i  and  m  be  unequal  positive  integers, 


And 

The  following  is  a  proof  of  the  first  property, 
We  have 


Multiplying  the  first  of  these  equations  by  Pm,  the  second 
by  Pit  subtracting  and  integrating,  we  get 


4-  [i  (i  +  1)  -  m  (m  +  1)}   P4P«^*  =  0. 

Hence,  transforming  the  first  two  integrals  by  integration 
by  parts,  and  remarking  that 

i(i+l)  -w(m  +  l)=(i-w)  (i 
we  get 

fn     ,rt/_ 

""(    M)" 


4-  (i  -  m)  (i  +  m  + 1)  \  Pfmd^,  =  0, 
or 

t-yTi-Pi-7fl)  +  (*-w)(*  +  m  + 


since  the  second  term  vanishes  identically. 
F.E. 
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Hence,  taking  the  integral  between  the  limits  -  1  and 
+  1,  we  remark  that  the  factor  1  -//,'  vanishes  at  both  limits, 
and  therefore,  except  when  i  —  m,  or  i+m+I  =  0, 


We  may  remark,  also  that  we  have  in  general 


r 

\ 

J 


P  *P<     p  dP 

*-  J-     -- 


a  result  which  will  be  useful  hereafter. 

11.     We  will  now  consider  the  cases  in  which 
i  -  m,  or  i  +  m  +  1  =  0. 


We  see  that  i+  m  +  1  cannot  be  equal  to  0,  if  i  and  m  are 
both  positive  integers.  Hence  we  need  only  discuss  the 
case  m  which  m  =  i.  We  may  remark,  however,  that  since 

Pi  =  P-d+u,  the  determination  of  the  value  of  f  Pz.  dp  will  also 

r1 

give  the  value  of  I     P.P^^  dp. 

rl 

The  value  of  I     Pfdp  may  be  calculated  as  follows  : 


Integrate  both  sides  with  respect  to  p  ;  then  since 
j  (1  -  2f,h  +  FT'  ^  =  -  1  log  (1  -  2/.A  +  A"), 
we  get,  taking  this  integral  between  the  limits  —  1  and  +  1, 


all  the  other  terms  vanishing,  by  the  theorem  just  pio/ed. 
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=  f1  p;dp  +  h*r  p&  +  ...  +  h"  I1  P 

J  -1  J  -1  .  J  -1 

Hence,  equating  coefficients  of  A2i, 


. 

12.     From  the  equation   I     Pf  'wtcfyt  =  0,  combined  witji 

the  fact  that,  when  yu,  =  1,  Pt  =  1,  and  that  Pt  is  a  rational 
integral  function  of  p,  of  the  degree  i,  Pt  may  be  expressed 
in  a  series  by  the  following  method. 

We  may  observe   in  the  first  place  that,  if  m  be  any 

f1 
integer  less  than  i,  I     /j^P^fjb  =  0. 

For  as  PTO,  PTO_1  .  .  .  may  all  be  expressed  as  rational  in 
tegral  functions  of  /LC,  of  the  degrees  m,  m  —  1  ...  respectively, 
it  follows  that  f  will  be  a  linear  function  of  Pm  and  zonal 
harmonics  of  lower  orders,  ^m~l  of  Pm-1  and  zonal  harmonics  of 


lower  orders,  and  so  on.     Hence    /^P/fyt  will  be  the  sum  of 
a  series   of   multiples  of  quantities  of  the  form   I  PmPidfi>t 

m  being  less  than  it  and  therefore  |     A^P^  =  0,  if  m  be  any 
integer  less  than  i. 
Again,  since 


it  follows,  writing  —  h  for  h,  that 

i-P.-P1A+  ...  +  (-1W+  - 

2—2 
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And  writing  —  p  for  JJL  in  the  first  equation, 


P0',  P/...P/...  denoting  the  values  which  P0,  P,  ...  P,, 
respectively  assume,  when  —  p,  is  written  for  p.  Hence 
P/  =  Pf  or  —  P^  according  as  i  is  even  or  odd.  That  is, 
Pi  involves  only  odd,  or  only  even,  powers  of  L^  according 
as  i  is  odd  or  even*. 
Assume  then 


Our  object  is  to  determine  At,  At_z.... 

Then,  multiplying  successively  by  /z,'"2,  fjL*~*,  ...  and  inte 
grating  from  —  1  to  + 1,  we  obtain  the  following  system  of 
equations : 

L_     _L  *'-2       i  i 

o  *  O      *     £^  •  p     I      ...     | 


2i_2s-3 


+            -*"S-2 i              _i               -4-L/-2»  t  fv 

— r — • -p  »  •  •  "r  — : ~^ i"  •  •  •  -11  W. 


And  lastly,  since  P4  =  1,  when  p  =  I, 


the  last  terms  of  the  first  members  of  these  several  equa 
tions  being 


A         A        A 

>  —   ~  '  ^o»  if  *  be 


13.  The  mode  of  solving  the  class  of  systems  of  equa 
tions  to  which  this  system  belongs  will  be  best  seen  by 
considering  a  particular  example. 

*  This  is  also  evident,  from  the  fact  that  Pt  is  a  constant  multiple  of 
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Suppose  then  that  we  have 


6-fa      c+  a 


From  this  system  of  equations  we  deduce  the  following, 
6  being  any  quantity  whatever, 


a+tf      ft  +00+0      «  (o)-a)(ft)-ye) 

For  this  expression  is  of  —  1  dimension  in  a,  b,  c,  a,  &  7, 
6,  a)  ;  it  vanishes  when  0  =  a,  or  0  =  /3,  and  for  no  other 

finite  value  of  6,  and  it  becomes  =  —  ,  when  6  —  to. 

0) 

We  hence  obtain 


-i.f.<,-ifFi  I  - 

~^&+0     c+0/  "  o>   o>-a«- 


and  therefore,  putting  6  =  —  a, 


(a  —  b)  (a  —  c)         (to  —  a)  (a>  —  /5) 
with  similar  values  for  y  and  z. 

And,   if  to  be  infinitely  great,  in  which  case  the   last 
equation  assumes  the  form  x+y  +  z=l,  we  have 

(I  ^M\  f  I  Q\ 

a  +  a)  (a  +  p) 

iC  =  j ~~ ^  , 

with  similar  values  for  y  and  z. 

14.     Now  consider  the  general  system 

_*t-4 

<*i  +  *i 

Xi 
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»  i  '«'-2  i  _j_  »'-g« 

T  •••  T    ' 


a,  +  ft)      a._2  +  a)  af_2s  +  &)  ft) 

the  number  of  equations,  and   therefore  of  letters  of  the 

<?       I       "I  n 

forms  x  and  «,  being  — —  if  i  be  odd,  5  + 1  if  t*  be  even ;  and 
the  number  of  letters  of  the  form  a  being  — ; - —  if  i  be  odd, 

n 

and  -  —  1  if  i  be  even. 
We  obtain,  as  before, 


...fe_2s-f  (9)... 

and,  multiplying  by  e^_2,  +  ^,  and  then  putting  6  —  —  a,._21, 
gj-y  +  Q  •  •  •  fa-»  +  gi-2>)  •  >  • 


= 


^     a,        («-«i)(ft,-«,2)...(a>-af.2J... 

,_,+  w)...  (a,_2,-f  ft))... 


-  a  or 


15.  To  apply  this  to  the  case  of  zonal  harmonics,  we  see, 
by  comparing  the  equations  for  x  with  the  equations  for  -4, 
that  we  must  suppose  w  =  GO  ;  and 

«i  =  i  «,-_2  =  i  -  S,...^  =  i  -  2s  ... 
a,  =  i-l,  a,._2  =  i  -  3,...af_2<  =  i  -  2s-  1... 
Hence 

_  (2i-2g-l)  (2»'-2g-3)...(2(t'-2g)-.l{... 
*-2'  "  (-2s)  {-(2*-2)}...{(t-2*-l)  or  (»-2»)} 


.4..,-25-l    or  t-2»)' 
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Or,  generally,  if  i  be  odd, 

i  =  O     /i        T»         i'S ~    » 


2,4...(t-3)x2  ' 


2.4...(f-5)x2.4     ' 


And,  if  i  be  even, 


2..4...(i-2)  x2 

_(2i-S)  (2i-7)...(^-3) 
f-4~     2. 4... (i-i-4)  x  2.4      ' 


2.4...* 


We  give  the  values  of  the  several  zonal  harmonics,  from 
P0  to  P10  inclusive,  calculated  by  this  formula, 


2      J 
5    3     3 
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7.5    4     5.3   2     3.1 


9.7   5     7.5         5.3 
ST*1*      272*+  274** 


2.4.6 

231/Lt6  -  315//,4  + 


8 
6     9.7.5     4  t   7.5.3     2     5.3.1 


P  ^13.11.9   7     11.9.7   5      9.7.5    3     7.5.3 

'    2.4.6   ^  ""2.  4x2^  +2T274/"  ""27476^ 
_  429//7  -  693/^5  +  315//.3  -  35/^ 
16 

=  15.  13.  11.  9    8    13.11.9.7  6     11.9.7.5   4 
2.4.6.8    ^  "  2.4.6  x2  ^      2.  4x2.  4^ 

9.7.5.3    2  f  7.5.3.1 
2x2.4.6^   '  2.4.6.8 

_  6435^8  -  12012/  +  6930/^4  -  1260/>62  +  35 
"T."28~~  ~  ' 

p  _17.15.13.11   9     15.13.11.9   »  ,  13_.11_.9._7   5 
2.4.6.8     ^    "   2.4.6x2  **  f  2.  4x2.4  ^ 

11.9.7.5   8     9.7.5.3f 
^"2x2.  4.  6^      2.4.6.8^ 
_  12155/*9-  25740/^7  +  18018^-  4620^+  315/x 
128  ~' 

P  =19.17.15.13.11    10    17.15.13.11.9   8    15.13.11.9.7  6 
10      2.4.6.8.10    *  •     2.4.6.8x2  "^+2.4.  6x2.  4^ 
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13.11.9.7.5   4     11.9.7.5.3   ,      9.7.5.3.1 
""2.4x2.  4.  6^2x2.  4.  6.  8.^      2.4.6.8.10 

_  46189//.10-  109395/*8+  9009(V;-  30030/4-  3465/u.2-  63 

256 

It  will  be  observed  that,  when  these  fractions  are  reduced 
to  their  lowest  terms,  the  denominators  are  in  all  cases 
powers  of  2,  the  other  factors  being  cancelled  by  correspond 
ing  factors  in  the  numerator.  The  power  of  2,  in  the 
denominator  of  P4,  is  that  which  enters  as  a  factor  into  the 
continued  product  1  .  2...i. 

16.     We   have  seen  that   /    p,™  P.  .  dp  =  0,  if  m  be  any 

J  -i 
integer  less  than  i. 

It  will  easily  be  seen  that  if  m  +  i  be  an  odd  number,  the 
values  of  lpfmPi.3ft  are  the  same,  whether  p,  be  put  =  l  or 

—  1  ;  but  ifm  +  i  be  an  even  number,  the  values  of  I/*"1  Pi  .  dp 

corresponding  to  these  limits  are  equal  and  opposite.   Hence, 
(m  +  i  being  even) 


J 


-1 
and 


ri      • 

then          p,mP..dp,  =  Q,  if  w  =  /-2,  i- 
Jo 


/•i 
We  may  proceed  to  investigate  the  value  of  I    p,m  P.  .  dp, 

Jo 

if  m  have  any  other  value.  For  this  purpose,  resuming  the 
notation  of  the  equations  of  Art.  13,  we  see  that,  putting 
9  —  m  +  1,  and  o>  =  oo  ,  we  have 


^i j ~i-2 i  I •*'*-** ^_ 

at  +  m  +  1      a._2  +  m  +  1  af_2,  +  m  +  1 

-  O  ...  (m  +  1  -  Q  ...  . 
1)...' 
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and    therefore,    putting    xt=At...,    «,=  *...,    a,=  i-l..., 

we  get 


o  f-t-m-i-i      ^— 

.  i  4.  $\       (m  —  l  \ 

if  i  be 


(m-i+  2)  (m  -  1+  4) .. 


.      . 

lf  *  be  even" 


In  the  particular  case  in  which  772  =  i,  we  get 

--   (i'odd)' 


24       * 
and  =—:  -  !  —  — 

' 


1  7.  We  may  apply  these  formal  03  to  develope  any  positive 
integral  power  of  /^  in  a  series  of  zonal  harmonics,  as  we 
proceed  to  shew. 

Suppose  that  m  is  a  positive  integer,  and  that  X*  is  de 
veloped  in  such  a  series,  the  coefficient  of  Pi  being  C{t  so 
that 


then,  multiplying  both  sides  of  this  equation  by  P.  and  inte 
grating  between  the  limits  -  1  and  1,  all  the  terms  on  the 

r1 

right-hand  side  will  disappear  except       Ci  P2t  dp,  which  will 


2 
become  equal  to  ^  —  -  Ct. 

Zil  ~T~  J. 

Hence  C  = 


which  is  equal  to  0,  if  m  +  i  be  odd.     Hence  no  terms  appear 
unless  m  4-  i  be  even.     In  this  case  we  have 
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Hence  the  formula  just  investigated  gives 

(m-t  +  2x)(m-^  +  4^(m-l) 
Ct  =  (2t  +  1)  (m  +  4-4Tjy  (m  +  i-l)7Ti(m  +  4)  (ro  +  2) 

if  /  be  odd,  and 

(wi-*  +  2)(»»-i  +  4)...J»  _ 
Ct  =  (2i  +  1)  ^  +  ^  +  jy^jTJITi)  .  .  .  (m  +~3)  (m  +  1) 

if  i  be  even. 

Therefore  if  m  be  odd, 

2.4.6...(ra-l)  p 


,    rr  ^-1  p  _A_p 

'  (m  +  4)(m  +  2j     3^m  +  2 


If  m  be  even, 

2.4. 

pT  =  (2m  +  1)  (2m+l)  (2m  -1) 


, 


P 


Hence,  putting  for  w  successively  0,  1,  2  ...  10,  we  get 


2  1          2          lp 

K O  J-2*~rtJ-Q  —   QZ2~rQO 

O  .  O  O  O  O 


/3  -  7  P  4-       P 

"77.5    3  +  5     » 
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42!)     7     39     5     33    3     3     » 

9  2.4.6.8  4.0.8 

17.15.  13.  11.  g^"4  'J  15.13.  lO7' 


128  64  48  40          1 

6435  ^  +  495  ^'  +  143  P«  +  9S  P*  +  9 


=  J?Lp  +  1»2  p      16  p       56  3 

12155     "^2431     7  +  C5     £  +  143     »  +  ri" 


=  21-   ------_-P   +17        -.. 

21.19.17.15.13.11    M+    '  19.17.15.  IS.ll^8 

113     6-8'10     P  |3     8-10    P+5    10    P+*P 
17.15.  IS.!!-*^9  15.  13Tn^  +  5i3jI7!i  +  lIP'> 

—  P  +—  8  P  +  —  P  +  48  P  +  5°  P  4    1  P 
46189    W+2717S     187    6  +  143    4  +  143    2  +  n° 
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18.  Any  zonal  harmonic  Pt  may  be  expressed  in  a  finite 
series  of  cosines  of  multiples  of  6,  these  multiples  being 
id,  (i-2)0....  Thus 


therefore,  writing  cos  0  for  /x,  and  observing  that 

1  -  2  cos  0h  +  ti*  =  (1  -  Ae^  «)  (1  -  he-  v-^), 
we  obtain 


-  he-  vriT*  =  P0 


or 


whence,  equating  coefficients  of 


(1.3...  (i-l)}2  .f  ., 
the  last  term  being  -j  —     ^  —  :  —  j-  if  %  be  even,  and 

1.3...(i  +  l)  1.3...  (i-  2)  ,,.,., 

iTElVnj  2.4...(i-i)  2  cos  ft  lf  ^  be  odd' 

19.     Let  us  next  proceed  to  investigate  the  value  of 

r  P  cos  mQ  sin 
Jo 
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This  might  be  done,  by  direct  integration,  from  the  above 
expression.  Or  we  may  proceed  as  follows. 

The  above  value  of  P.  when  multiplied  by  cos  m6  sin  0 

(that  is  by  ^  {sin  (ra  +  1)  0  —  sin  (m  —  1)  0j)  will  consist  of  a 

series  of  sines  of  angles  of  the  form  (i  —  2n  ±  (m  ±  1)}  6,  that 
is  of  even  or  odd  multiples  of  #,  as  i  +  m  is  odd  or  even. 
Therefore,  when  integrated  between  the  limits  0  and  TT  it 
will  vanish,  if  i  +  m  be  odd.  We  may  therefore  limit  our 
selves  to  the  case  in  which  i  +  m  is  even. 

Again,  since  cos  mO  can  be  expressed  in  a  series  of  powers 
of  cos  0,  and  the  highest  power  involved  in  such  an  expression 
is  coswt#,  it  follows  that  the  highest  zonal  harmonic  in  the 

development  of  cos  m6  will  be  Pm.  Hence  /  P.  cos  m6  sin  6  dd 
will  be  =0,  if  m  be  less  than  i. 


Now,  writing 


<7f_2cos 


we  see  that  P{  cos  m0  sin  0  dO  will  consist  of  a  series  of  sines 
of  angles  of  the  forms  (m  +  t+V)  0,  (m  +  i  —  1)  0  ...  down  to 
(m  —  i—V)0,  thqre  being  no.  term  involving  m0,  since  the 
coefficient  of  such  a  term  must  be  zero.  Hence 

I   Pi  cos  m0  sin  0  d0, 

Jo 

will  consist  of  a  series  of  fractions  whose  denominators  in 
volve  the  factors  m  +  i  +  1,  m  +  i  —  1  ...  m  —  i  —  1  respectively. 
Therefore  when  reduced  to  a  common  denominator,  the  result 
will  involve  in  its  denominator  the  factor 

(m  +  i+1)  (m+i-V)  ..,  (w  +  1)  (m-1)  ...  (m-/-l) 
if  m  be  even,  and 

(m .+  i  +  1)  (m  -J-  z  —  1)  ...  (m  +  2)  (m  —  2)  ...  (m  —  i  —  1) 
if  m  be  odd. 

For  the  numerator  we  may  observe  that  since 

r 

Pi  cos  m  0  sin  OdB 
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vanishes  if  m  be  less  than  i,  it  must  involve  the  factors 
m  -  (i  -2),  m  -  (i  -  4)  .  .  .  m  +  (i  -  2),  and  that  it  does  not 
change  sign  with  m.  Hence  it  will  involve  the  factor 

[m  -  (t  _  2)}  {m  -  (i  -  4)}  ...  (m  -  2)  m2  (m  +  2)  ...  (m  +  i-  2) 
if  m  be  even,  and 

{w  -  (t'-  2)}  {m-  (i-  4)}  ...  (m  -  1)  (m  +  1)  ...  (m  +  i-  2) 
if  m  be  odd. 

To  determine  the  factor  independent  of  m.  we  may  pro 
ceed  as  follows  : 

Pt  =  C.COS  i6  +  0;._2  COS  (i-  2)0+  ...  ; 

..'.  Pi  cos  7710  =  ^0,  {cos  (m  +  1)  0  +  cos  (m  -  i}  9} 

+  g  Q_2  {cos  (m  +  i-  2)  0  +  cos  (m  -  i  +  2)  0}  +  ...  ; 

.*.  P.cosm^sin  6  =  j  Q  {sin  (m  +  i  +  1)  0-  sin  (n»  +  i-  1)  0 

+  sin  (m  -  i  +  l)  0  -  sin  (m  -  i  -  1)  6} 
+  ^  CL2  (sin  (??z  +  ?;—  1)  6  —  sin  (w  +  i-  3)  0 


»  — t  +  3)0-sin  (m  — *  +  l)  0}  +  ... ; 

/•IT 

.'.   I    P^osmO  sind  dO 

=^{_J_       J_  . ,  _  _j_      _i i 

%  \m  +  i  +  l     m  +  i-1     m-i+l     m  —  {—  i  j 


__       ___4 


__  ___ 

2    [??i  +  *  -  1     m  +  i  -  3     m-i+3  m  —  i  +  1 

|  __  i^L_.  +  -    ^'-1     1 

']     »!»-(*>  I)1  "mM-.(*-  1)1} 

] 

'  *" 


m2  -  (t-  I)2     wi"  -  (i  -  3)' 
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Now,  when  m  is  very  large  as  compared  with  i,  this  be 
comes 


since  Ci  4-  C{_z  4-  ...  =  1,  as  may  be  seen  by  putting  6  =  0. 

[*  2 

Hence   I    Pi  cos  m6  sin  6  dd  tends  to  the  limit  -     .  ,  as  m 
Jo  m 

is  indefinitely  increased. 

The  value  of  the  factor  involving  m  has   been  shewn 
above  to  be 

{m  -  (i  -  2)}  {m  -  (i-  4)}  ...  (m  -  2)  m2  (m  +  2)  ...  (m  4-  1  -  2) 

{m  -  (»'+!)}  {m  -  (i  -  1)}  ...  (771  -  1)  (m  +  1)  ...  (m  +  1  +  1) 
if  m  be  even,  and 


{m  -  (i+  1)J  m  -  (i-  1)J  ...  (m  -  2)  (m  +  2)  ...  (m  +  %+  1) 
if  m  be  odd. 

Each  of  these  factors  contains  in  its  numerator  two  factors 
less  than  in  its  denominator.     It  approaches,  therefore,  when 

m  is  indefinitely  increased,  to  the  value  —  .     Hence 


J 


if  m  and  i  be  even,  and 
_    \ 


-  (»  +  1)}  {m-  (4-  1)}  ...  (m  -2)  (m+2)  ...  {m+ 
if  m  and  i  be  odd. 

In  each  of  these  expressions  i  may  be  any  integer  such 
that  m  —  i  is  even,  i  being  no£  greater  than  m.  Hence  they 
will  always  be  negative,  except  when  i  is  equal  to  m. 
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20.     We  may  apply  these  expressions  to  develop  cos  mO 
in  a  series  of  zonal  harmonics.  ; 

Assume 


Multiply  by  Pt  sin  0,  and  integrate  between  the  limits  0 
and  TT,.  and  we  get 

_     {m-(t 


(;+!)}  jm     (,-_  i)}  ...  [W  4.  (^ 
Hence 


=  -  (2i  4- 


Hence,  putting  m  successively  =  0,  1,  2,  ...  10, 
cosO<9  =  P 


4  1 

=-  P  --  P 
3    2     3    0 


-86*    -21^-15 

F.  H. 
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„„  2.4.6.8  4.6 


cos  60  =  -1.3 


0 


"  63     5     9s     7    1J 

2. 4. 62. 8. 10 
-1.1.3.5.7.9.11.13    6 

4.62.8 


____ 
1.3.5.7.9.11    4       3.5.7.9    2     5.7 


1! 


___ 
231    6      385    4     21    2     35    0 

2.4.6.8.10.12 
5-1.  1.8.  5.9.11.13.15  P' 

4.6.8.10  6.8 


1.3.5.9.11.13    5        3.5.9.11    3     5.9 
1024          128          112          ^ 
429     7      117    5     495    3     15    l' 

2.4.6.8s.  10. 12. 14 
-  1 . 1.3.5.7.9.11.13.15 . 17 


on  -,_  Z.  ^.  O  .0'  .  1U  .  1Z.  l^fc  „ 

cos  80  =  - 17  — .    .    „    „   ^   „   ^    ^    ^   ^  P8 


____          ___ 

1  .  3  .  5  .  7  .  9  .  11  .  13  .  15    6     y  3  .  5  .  7  .  9  .  11  .  13    4 


2 


_   %  P    _  P 

5.7.4.11    "     7.9    ° 


=  16384    4096    J>56     64  p_^p 
"  6435     3465    1001  4  693  2  63  ° ' 

2.4.6.8.10.12.14.16 
COS9^"19-1.1.3.5.7.11.13.15.17.19P- 

4.6.8.10.12.14          6.8.10.12 
° 1.3.5.7.11.13.15.17  7    3.5.7.11.13.15  5 

_7___8_i!°_  P^s-1  P 
5.7.11.13  3   7.11  : 
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32768    3072  p_  128    _16__S 
12155  9  2431  7  455  6  143  s  77  *' 

2.4.6.8.10M2.14.16.18 


-     21 
^ 


-1.1.3.5.7.9.11.13.15.17.19.21 
4.6.8.10M2.14.16 


__17 


1.3.5.7.9.11.13.15.17.19 
6.8.10M2.14  8.10M2 

^3.5.7.9.11.13.15.17    6     y  5  .  7.9.  11  .  13.T5-    4 

_5       iQ2      P_  _1    P 

7.9.11.13    2     9.11    ° 

131072       _  32768  p  __^12  _      JL28  p       500   p 
"  46189     10     24453    8     1683    6  ~"  1001    4  ""  9009    2 

-^P 

99    °* 

21.  The  present  will  be  a  convenient  opportunity  for 
investigating  the  development  of  sin  6  in  a  series  of  zonal 
harmonics.  Since  sin  6  =  (1  -  yu,2)i,  it  will  be  seen  that  the 
series  must  be  infinite,  and  that  no  zonal  harmonic  of  an  odd 
order  can  enter.  Assume  then 

sin0=  C0P0+  02P2  +  ...  +  <74P,+  ... 
i  being  any  even  integer. 

Multiplying  ^by  Pi?  and  integrating  with  respect  to  ^ 
between  the  limits  —  1  and  +  1,  we  get 


supposing  Pi  expressed  in  terms  of  the  cosines  of  6  and  its 
multiples 


3—2 
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Hence,  putting  i  =  0, 


I  O 

Putting  i  =  2,  and  observing  that  P2  =  -  +  -  cos  20, 

cos20)^ 


~  J*  jl  +  2  cos  20  - 1  (1  +  cos  40)1 


—a*- 


For  values  of  i  exceeding  2,  we  observe,  that  if  we  write 
for  Pi  the  expression  investigated  in  Art.  18,  the  only  part 

of  the  expression  /  Pi  (1  -  cos  20)  d9  which  does  not  vanish 

J  o 

will  arise  either  from  the  terms  in  Pt  which  involve  cos  20,  or 
from  those  which  are  independent  of  0.     We  have  therefore 


f 

Jo 


-2)  2.4  ...  (»— 2) 

077^}]  P-0082^ 


4      "2.4...      t       2,4...  (t-  2) 
* 


1.3...  (»-l)  1.3...  (i-3) 


4  .  2.4...      i       2.4...(t-2)      V    t 
2t  +  l  1.3...  (t'-l)  1.3...  (t-3) 
2      2.4...»(*  +  2)  2.4...(»-2)«" 

Hence    sin0  =  ^P0-gp2- ... 

__  (2t'+l)ir  1.3...  (t-1)  1.3...   (f-3) 
2         2.4...i(»  +  2)  2.4...  (t-2)i 
2*  being  any  even  integer. 
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22.     It  will  be  seen  that  -~  ,  being  a  rational  and  integral 

function   of  //T1,   //~3...,   must   be  expressible  in  terms  of 
P,_j,  Pi_3...     To  determine  this  expression,  assume 


then  multiplying  by  Pm,  and  integrating  with  respect  to 
from  —  1  to  +  1, 


And 

Now,  since  i  >  m, 


fm        ^  =  [^.PJ1  -  [^-PJ4  =  2, 

since  either  m  or  i  must  be  odd,  and  therefore  either  Pm  or 
Pt  =  —  1,  when  p  =  —  1  ; 


.•.'=(2»-l)PM  +  (2*  -5)  P<_  +  ( 
Hence  ^  -  ^  =  (2«-  1) 

rfft  ^ 

23.     From  this  equation  we  deduce 


the  limits  /u-  and  1  being  taken,  in  order  that  Pt  —  Pf_2  may 
be.  equal  to  0  at  the  superior  limit. 
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Now,  recurring  to  the  fundamental  equation  for  a  zonal 
harmonic,  we  see  that 


24.     We  have  already  seen  that  I    Pfm  dp,  =  0,  i  and  m 

J  -i 
being  different  positive  integers.     Suppose  now  that  it  is 

required  to  find  the  value  of  I    Pfm  dp. 

J  it. 


We  have  already  seen  (Art.  10)  that 

a  -•!* 

(    M 


And,  from  above, 

3cZP,_ 
1       /ij~ 


m(m+I),p         p 
2m  +1     (    m+l 


p 
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25.  We  will  next  proceed  to  give  two  modes  of  ex 
pressing  Zonal  Harmonics,  by  means  of  Definite  Integrals. 
The  two  expressions  are  as  follows : 


These  we  proceed  to  establish. 
Consider  the  equation 


-T 

TTJo  a  — 


The  only  limitation  upon  the  quantities  denoted  by  a 
and  b  in  this  equation  is  that  62  should  not  be  greater  than 
a2.  For,  if  62  be  not  greater  than  a2,  cos  S-  cannot  become 

equal  to  ^  while  Sr  increases  from  0  to  TT,  and  therefore  the 
expression  under  the  integral  sign  cannot  become  infinite. 

Supposing  then  that  we  write  z  for  a,  and  V—  1  p  for  6, 
we  get 

d^  1 


We  may  remark,  in  .passing,  that 

f         ^ 

Jo  s  —  V—  l 


and  is  therefore  wholly  real. 

Supposing  that  p2  =  #2  +  y*,  and  that  x*+  y*  +  t?  =  ^2, 
thus  obtain 

cfe  __  1 

z  —  V—lpcos1^     r 
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Differentiate  i  times  with  respect  to  zt  and  there  results 


Hence  P.  - 


TT     o  (s- 
In  this,  write  fir  for  z,  and  (1  -  ^  r  for  p,  and  we  get 


w 


Inch,  writing  TT  —  Sr  for  S,  gives 
P  —  - 

JL   j     — — 


*  +&**-!)  cos*}"1 

26.    Again,  we  have 

1  1 


In  this  write  1  -^  for  a,  and  ±  (ya2  -  1)4  h  for  6,  which  is 
admissible  for  all  values  of  h  from  0  up  to  ^~  (//,2-l)^,  and 
we  obtain,  since  a2  -  62  becomes  1  -  2/xfc  H-  7i2, 

1 


ZONAL  HARMONICS. 

Hence,  equating  coefficients  of  h*, 


41 


The  equality  of  the  two  expressions  thus  obtained  for  P.  is 
in  harmony  with  the  fact  to  which  attention  has  already 
been  directed,  that  the  value  of  Pi  is  unaltered  if  —  (i  +  1) 
be  written  for  i. 


27.  The  equality  of  the,  two  definite  integrals  which 
thus  present  themselves  may  be  illustrated  by  the  following 
geometrical  considerations. 


CQ*  =  a2  +  62  -  2a&  cos 


Hence 


Let   0  be  the  centre  of  a  circle,  radius  a,  C  any  point 
within  the  circle,  POQ  any  chord  drawn  through  (7,  and  let 


=(a2- 

=  0. 


-    a   cos 


—  zao cos 
sin 


»  _  2a6  cos  S- 
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Again,  since  the  angles   OP C}   OQG  are  equal  to  one 
another, 

sinj?^  sin  OPC _  sin  OQG     sinifr 
CP  =       OO  00        ''~OQ  > 


(a-  +  62  -  2a&  cosS-)^     (a2  +  62  -  2a&  cos  ^  ' 


whence 


(a2  +  62  -  2ab  cos  $•)*  +  (a2  +  62  -  2a&  cos 
And  7  =  (a2  +  62-2a5cos 


_Q 


In  this,  write  a2  +  62  =  At,  2a&=  +  (yu,2-  1)*,  which  gives 
az  —  bz  =  1,  and  we  get 


=  "  {fl  ±  ^  ~  1)f  cos  tN^ 


"We  also  see,  by  reference  to  the  figure,  that  as  ^  in 
creases  from  0  to  TT,  ->|r  diminishes  from  TT  to  0.     Hence 


28.  From  the  last  definite  integral,  we  may  obtain  an  ex 
pansion  of  Pt  in  terms  of  cos  6  and  sin  0.  Putting  JJL  =  cos  6, 
we  get 


+  (cos  0  -  V^I  cos  ^  sin 
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)..>W—  ;  ^cog  ^2yn  £cog  0y-2Kt  £gjn  Q\jlm 


Now 
And 


1.2..  .2m  2m(2m-2)...2 

_  i*  (*'-  1)  •••.(*'-  2m  +  1) 

(2.4..,2m)a 

,  =  (cos  ^)f  -  ^^  (cos  0f*  (sin  ^)2  +  .  .  . 


CHAPTER  III. 

APPLICATION  OF  ZONAL  HARMONICS  TO  THE  THEORY  OF 
ATTRACTION.  REPRESENTATION  OF  DISCONTINUOUS 
FUNCTIONS  BY  SERIES  OF  ZONAL  HARMONICS. 

1.  WE  shall,  in  this  chapter,  give  some  applications  of 
Zonal  Harmonics  to  the  determination  of  the  potential  of  a 
solid  of  revolution,  symmetrical  about  an  axis.  When  the 
value  of  this  potential,  at  every  point  of  the  axis,  is  known, 
we  can  obtain,  by  means  of  these  functions,  an  expression 
for  the  potential  at  any  point  which  can  be  reached  from 
the  axis  without  passing  through  the  attracting  mass. 

The  simplest  case  of  this  kind  is  that  in  which  the 
attracting  mass  is  an  uniform  circular  wire,  of  indefinitely 
small  transverse  section. 

Let  c  be  the  radius  of  such  a  wire,  p  its  density,  k  its 
transverse  section.  Then  its  mass,  Mt  will  be  equal  to  2-Trpc/J, 
and  if  its  centre  be  taken  as  the  origin,  its  potential  at  any 

point  of  its  axis,  distant  z  from  its  centre,  will  be  —  -  ^. 

Now,  this  expression  may  be  developed  into  either  of  the 
following  series  : 

Vra?"1"^?"-"1"^"1^  l!"L2Ti7."^r~c"5I+-'T''(1)> 

ML     lc*     1.3c«  '  .         1.3...(2t--l)c" 


We  must  employ  the  series  (1)  if  z  be  less  than  c,  or 
if  the  attracted  point  lie  within  the  sphere  of  which  the  ring 
is  a  great  circle,  and  the  series  (2)  if  z  be  greater  than  c, 
or  if  the  attracted  point  lie  without  this  sphere. 
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Now,  take  any  point  whose  distance  from  the  centre  is  r, 
and  let  the  inclination  of  this  distance  to  the  axis  of  the 
ring  be  6.  In  accordance  with  the  notation  already  em 
ployed,  let  cos  0  ~  pt  Then,  the  potential  at  this  point  will 
be  given  by  one  of  the  following  series  : 


j  x  o 


- 

2  .  4    V 

1.8.5...2»-l        c" 


2.  4.  6..  .2i 


For  each  of  these  expressions,  when  substituted  for  V, 
satisfies  the  equation  V2  V  =  0,  and  they  become  respectively 
equal  to  (1)  and  (2)  when  6  is  put  =  0,  and  consequently 
r  —  z.  The  expression  (2')  also  vanishes  when  r  is  infinitely 
great,  and  must  therefore  be  employed  for  values  of  r  greater 
than  c,  while  (!')  becomes  equal  to  (2')  when  r  =  c,  and  will 
therefore  denote  the  required  potential  for  all  values  of  r 
less  than  c. 

These  expressions  may  be  reduced  to  other  forms  by 
means  of  the  expressions  investigated  in  Chap.  2,  Art.  25,  viz. 


i  =  i  f  (^  +  j~js=i  cos  ^y  a-, 

TTJo 


or  P.  =  -       (u,  +  l  cos  ^~(i+1)  cZik 

TTJo 

Substitute  the  first  of  these  in  (!')  and  (observing  that 
fj,r  =  z]  we  see  that  it  assumes  the  form 

M    «  (      1 
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which  is  equivalent  to 

Mr 


The  substitution  of  the  last  form  of  Pi  in  the  series  (2') 
brings  it  into  the  form 


which  is  equivalent  to 


2.  Suppose  next  that  the  attracting  mass  is  a  hollow  shell 
of  uniform  density,  whose  exterior  and  interior  bounding 
surfaces  are  both  surfaces  of  revolution,  their  common  axis 
being  the  axis  of  z.  Let  the  origin  be  taken  within  the 
interior  bounding  surface  ;  and  suppose  the  potential,  at  any 
point  of  the  axis  within  this  surface,  to  be 


Then  the  potential  at  any  point  lying  within  the  inner 
bounding  surface  will  be 

A0PQ  +  A^r  +  Aff  +...+  Aff  +  ... 

For  this  expression,  when  substituted  for  V,  satisfies  the 
equation  V2F=0;  it  also  agrees  with  the  given  value  of 
the  potential  for  every  point  of  the  axis,  lying  within  the 
inner  bounding  surface,  and  does  not  become  infinite  at  any 
point  within  that  surface. 

Again,  suppose  the  potential  at  any  point  of  the  axis 
without  the  outer  bounding  surface  to  be 
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Then  the  potential  at  any  point  lying  without  the  outer 
bounding  surface  will  be 


r  r  r 

For  this  expression,  when  substituted  for  V,  satisfies  the 
equation  V2  V—  0  ;  it  also  agrees  with  the  given  value  of  the 
potential  for  every  point  of  the  axis,  lying  without  the  outer 
bounding  surface,  and  it  does  not  become  infinite  at  any 
point  within  that  surface. 

By  the  introduction  of  the  expressions  for  zonal  har 
monics  in  the  form  of  definite  integrals,  it  will  be  found  that 
if  the  value  of  either  of  these  potentials  for  any  point  in  the 
axis  be  denoted  by  <f>  (z),  the  corresponding  value  for  any 
other  point,  which  can  be  reached  without  passing  through 
any  portion  of  the  attracting  mass,  will  be 


-  r 
IT  J  o 


o 

3.  We  may  next  shew  how  to  obtain,  in  terms  of  a  series  of 
zonal  harmonics,  an  expression  for  the  solid  angle  subtended 
by  a  circle  at  any  point.  We  must  first  prove  the  following 
theorem. 

The  solid  angle,  subtended  by  a  closed  plane  curve  at  any 
point,  is  proportional  to  the  component  attraction  perpendicular 
to  the  plane  of  the  curve,  exercised  upon  the  point  by  a  lamina, 
of  uniform  density  and  thickness,  bounded  by  the  closed  plane 
curve. 

For,  if  dS  be  any  element  of  such  a  lamina,  r  its  distance 
from  the  attracted  point,  6  the  inclination  of  r  to  the  line 
perpendicular  to  the  plane  of  the  lamina,  the  elementary 
solid  angle  subtended  by  d-8  at  the  point  will  be 

dS  cos  0 

'-V  ~~^       ' 

And  the  component  attraction  of  the  element  of  the 
lamina  corresponding  to  dS  in  the  direction  perpendicular 
to  its  plane  will  be 

7  dS       n 
p/c  —  „-  cos  0, 
r 
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p  being  the  density  of  the  lamina,  7c  its  thickness.  Hence, 
for  this  element,  the  component  attraction  is  to  the  solid 
angle  as  pk  to  1,  and  the  same  relation  holding  for  every 
element  of  the  lamina,  we  see  that  the  component  attraction 
of  the  whole  lamina  is  to  the  solid  angle  subtended  by  the 
whole  curve  as  pk  to  1. 

Now,  if  the  plane  of  scy  be  taken  parallel  to  the  plane 
of  the  lamina,  and  V  be  the  potential  of  the  lamina,  its 
component  attraction  perpendicular  to  its  plane  will  be 

cJ  T/ 

—  -=—  .     Now  since  V  is  a  potential  we  have  V2  V—  0,  whence 


=,  or  Hence  is  itself  a  potential, 

dz  \dz  )  dz 

and  satisfies  all  the  analytical  conditions  to  which  a  potential 
is  subject.  It  follows  that,  if  the  solid  angle  subtended  by 
a  closed  plane  curve  at  any  point  (x,  ?/,  z)  be  denoted  by 
ft),  a)  will  be  a  function  of  x,  y,  zy  satisfying  the  equation 
^ft)  =  0.  Hence,  if  the  closed  plane  curve  be  a  circle  it 
follows  that  the  magnitude  of  the  solid  angle  which  it  sub 
tends  at  any  point  may  be  obtained  by  first  determining 
the  solid  angle  which  it  subtends  at  any  point  of  a  line 
drawn  through  its  centre  perpendicular  to  its  plane,  and 
then  deducing  the  general  expression  by  the  employment 
of  zonal  harmonics. 

Now  let  0  be  the  centre  of  the  circle,  Q  any  point  on  the 
line  drawn  through  0  perpendicular  to  the  plane  of  the 
circle,  E  any  point  in  the  circumference  of  the  circle.  With 
centre  Q,  and  radius  QO,  describe  a  circle,  cutting  QE  in  L. 
From  L  draw  LN,  perpendicular  to  Q  0. 

Let  OE=c,  OQ  =  z. 


Then  EL  =  (c*  +  rf  -  z,  ON=  ~        {(c°  +  *>)*- 
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And  the  solid  angle  subtended  by  the  circle  at  G( 

ON 
=  4-77  -sr-r 


. 

To  obtain  the  general  expression  for  the  solid  angle  sub 
tended  at  any  point,  distant  r  from  the  centre,  we  first 
develope  this  expression  in  a  converging  series,  proceeding 
by  powers  of  •&.  This  will  be 

zlz*     1.3s-  1.3...(2i-l)gwi 

-^^-ra^-"  2.4...*      ?*+- 

if  2-  be  less  than  c,  and 


s2~~274s4  v       '        2.4...2i      2Tl       ") 

if  z  be  greater  than  c. 

Hence,  by  similar  reasoning  to  that  already  employed, 
we  get,  for  the  solid  angle  subtended  at  a  point  distant  r 
from  the  centre, 


c     '  2    c3        2 . 4    c5 

-  (-  1)'  ~2 . 4 . ,.2t   ;  "^r~~  +  ' ' ' 
if  r  be  less  than  c,  and 

2?r  (2  I?'"  ""  274  ~pr+'"~(~1^      2 747.72r~ 
if  r  be  greater  than  c. 

4.  We  may  deduce  from  this,  expressions  for  the  potential 
of  a  circular  lamina,  of  uniform  thickness  and  density,  at  an 
external  point.  For  we  see  that,  if  V  be  the  potential  of 
such  a  lamina,  k  its  thickness,  and  p  its  density,  we  have  for 
a  point  on  the  axis, 


F.  H. 
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whence  V=  Zirpk  {(c2  +  z*)*  -  z} 

=  f{('2+*2)*-4 

if  IT  be  the  mass  of  the  lamina. 

Expanding  this  in  a  converging  series,  we  get 
Iz*     l.l4      l/U3/ 


.....- 
1       '       2. 4. 6... 2*      a*'1 


if  z  be  less  than  c,  and 


ira  c2 _ LJ.c4    i.i.s 

~  + 


if  z  be  greater  than  c. 

Hence  we  obtain  the  following  expressions  for  the  po 
tential  of  an  uniform  circular  lamina  at  a  point  distant  r 
from  the  centre  of  the  lamina  : 

1    1  P  r4 


2.4.6.. .2i       cai'1 
if  r  be  less  than  c,  and 

IT  f  1  P/  __  1^1  P/     IJLp  ^ __ 

~  C2(2    r       2.4  r3    +2.4.6    r5  " 

V1.1.3...(2t-3)l^ 
"l"^       2.4.G..T2?"     r"-1   '     '" 

if  r  be  greater  than  c. 
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It  may  be  shewn  that  the  solid  angle  may  be  expressed 
in  the  form 


27T-2 


and  the  potential  of  the  lamina  in  the  form 


-  *      cos 


5.  As  another  example,  let  it  be  required  to  determine 
the  potential  of  a  solid  sphere,  whose  density  varies  inversely 
as  the  fifth  power  of  the  distance  from  a  given  external  point 
0  at  any  point  of  its  mass. 

It  is  proved  by  the  method  of  inversion  (see  Thomson 
and  Tait's  Natural  Philosophy,  Vol.  1,  Art.  518)  that  the 

potential  at  any  external  point  P'  will  be  equal  to  —  -  —  T,  0' 

being  the  image  of  0  in  the  surface  of  the  sphere,  and  M 
the  mass  of  the  sphere.  We  shall  avail  ourselves  of  this 
result  to  determine  the  potential  at  a  given  internal  point. 

^  Let  G  be  the  centre  of  the  sphere,  0  the  given  external 
point.  Join  CO,  and  let  it  cut  the  surface  of  the  sphere  in  A, 
and  in  GA  take  a  point  0',  such  that  CO  .  GO  =  GA\  Then 
0'  is  the  image  of  0. 

Let  P  be  any  point  in  the  body  of  the  sphere,  then  we 
wish  to  find  the  potential  of  the  sphere  at  P. 

Take  0  as  pole,  and  OC  as  prime  radius,  let  OP  =  r, 
POC  =  6.  Also  let  GA  =  a,  CO  =  c. 

Let  the  density  of  the  sphere  at  its  centre  be  p,  then  its 

c5 
density  at  P  will  be  p  —5  .     Hence 
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the  limits  of  r  being  the  two  values  of  r  which  satisfy  the 
equation  of  the  surface  of  the  sphere,  viz. 

r2  +  c2  -  2er  cos  0  =  a2, 

and  those  of  6  being  0  and  sin"1  -  . 

c 

Hence,  if  rl5  r2  be  the  two  limiting  values  of  r,  we  have 


Now 
Also 


r,  '  r,       c2  -  a*  ' 


.    /j.  _  1\  _     ^cos^-Cc'- a8)]' 
' 


Svrpc5       2c  2 


2^     (^TflS ' <?"_""*  j81     6  cos  0sin  ^ (a2- c2sm*0)^'  dd 

^47rpc6     pin-^ 
(c2-a2)2J0 


cos  e  sin  ^  ^2  -  <?  si 


Now,  if  F  be  the  potential  at  P,  we  have  (see  Chap.  i. 
Art.  1) 
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This  is  satisfied  by  V=  —  ^  — 3-. 

Assume  then,  as  the  complete  solution  of  the  equation, 


It  remains  to  determine  the  coefficients  A0,  A1...Ai...Bi), 
2?t....ffj,  so  that  this  expression  may  not  become  infinite  for 
any  value  of  r  corresponding  to  a  point  within  the  sphere, 
and  that  at  any  point  P  on  the  surface  of  the  sphere  it  may 

be  equal  to  jr-p,  where  O'P :  OP  ::  a  :  c,  and  therefore,  at 
the  surface, 


a    OP      3(c'-(Ay- 
And,  at  the  surface,  we  have 


_1_  /I 
~  ^-~tf  \    r 


1      2ccos 


C    -6 
P.  P. 


Pl  +  ...  identically. 


*-(-*&**)*• 


and  7?    B...B...A    A..A  all  =  0. 
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Hence  since  P0  =  1, 

2    vpf    f    M  } 

°     3c*-a2  l(c'-u')       j 


3  (cy 

and  A  =  |-^. 

3  c"  -  a 

whence  we  obtain,  as  the  expression  for  the  potential  at  any 
internal  point, 

v=  2     Trpc5      3a2  -  c2      4   TT/OC^  cos  0  _  2  Tr/oc5 
"3(c»"-asj*       r       *3cP-a-V        3  V    ' 

G.  We  shall  next  proceed  to  establish  the  proposition  that 
{/"  the  density  of  a  spherical  shell,  of  indefinitely  small  thick 
ness,  be  a  zonal  surface  harmonic,  its  potential  at  any  internal 
point  will  be  proportional  to  the  corresponding  solid  har 
monic  of  positive  degree,  and  its  potential  at  any  external 
point  will  be  proportional  to  the  corresponding  solid  harmonic 
of  negative  degree. 

Take  the  centre  of  the  sphere  as  origin,  and  the  axis  of 
the  system  of  zonal  harmonics  as  the  axis  of  z.  Let  b  be  the 
radius  of  the  sphere,  Bb  its  thickness,  U  its  volume,  so  that 
U=  4>7rb28b.  Let  CPi  be  the  density  of  the  sphere,  Pt  being 
the  zonal  surface  harmonic  of  the  degree  i,  and  G  any  con 
stant. 

Draw  two  planes  cutting  the  sphere  perpendicular  to  the 
axis  of  z,  at  distances  from  the  centre  equal  to  f,  f  +  c?? 
respectively.  The  volume  of  the  strip  of  the  sphere  inter 

cepted  between  these  planes  will  be  —^  U,  and  its  mass  will  be 

<7P,  ^U.  V  :   l/  -  '/at 

Now  f  =  bfjij  hence  d%=bd/ji,  and  this  mass  becomes 
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Hence  the  potential  of  this  strip  at  a  point  on  the  axis  of  z, 
distant  z  from  the  centre,  will  be 


CU  Pt  -, 

T-tf  +  y-  ^    "  ^ 


which  may  be  expanded  into 


§2>:  (_4r 
**    'Vr-v>  + 
<rx- 


£=£---- 


and    • 


To  obtain  the  potential  of  the  whole  shell,  we  must  inte 
grate  these  expressions  with  respect  to  //,  between  the  limits 
-  1  and  +  1.  Hence  by  the  fundamental  property  of  Zonal 
Harmonics,  proved  in  Chap.  II.  Art.  10,  we  get  for  the  po 
tential  of  the  whole  shell 


--  —  -  ^+1  at  an  internal  point, 
_  i  ~\~  J.  o 

CU    V 
—  T  ~i+i  a^  an  external  point. 

JL  Z 


From  these  expressions  for  the  potential  at  a  point  on 
the  axis  we  deduce,  by  the  method  of  Art.  1  of  the  present 
Chapter,  the  following  expressions  for  the  potential  at  any 
point  whatever  : 

CU  P/  .  .        .      .  . 

Vl  =  —.  —  ^  jjfi  at  an  internal  point, 

£1  ~F  L    0 

CU  PV 

Y  =    .        -~  at  an  external  point. 
2     2i  +  1  r 

From  hence  we  deduce  the  following  expressions  for  the 
normal  component  of  the  attraction  on  the  point. 

Normal  component  of  the  attraction  on  an  internal  point, 
measured  towards  the  centre  of  the  sphere, 

dv,        i    rTTP/-1 

(j  u  —  7       • 


dr         ^i  +  l         b* 
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Normal  component  of  the  attraction  on  an  external  point, 
measured  towards  the  sphere, 


=  _,_ 

dr  L  U         ' 


In  the  immediate  neighbourhood  of  the  sphere,  where  r  is 
indefinitely  nearly  equal  to  b,  these  normal  component  at 
tractions  become  respectively 


P         i  4-  1  P 

*        ^  * 

-F*. 


and  their  difference  is  therefore 


And  writing  for   Z7  its  value,  4nrl*S\  this  expression  be 
comes 


Or,  the  density  may  be  obtained  by  dividing  the  alge 
braic  sum  of  the  normal  component  attractions  on  two  points, 
one  external  and  the  other  internal,  indefinitely  near  the 
sphere,  and  situated  on  the  same  normal,  by  4?r  x  thickness 
of  the  shell. 

7.  It  follows  from  this  that  if  the  density  of  a  spherical 
shell  be  expressed  by  the  series 


<70,  C^,  (72  ...  Gi  ...  being  any  constants,  its  potential  (FJ  at 
an  internal  point  will  be 


and  its  potential  ( P8)  at  an  external  point  will  be 

(Cjp,   i  cw    i  cnpp  i 

J(    r     ^3~^~  +  5~f'i''"+2i^ 

In  the  last  two  Articles,  by  the  word  "density"  is  meant 
"volume  density,"  i.e.  the  mass  of  an  indefinitely  small 
element  of  the  attracting  sphere,  divided  by  the  volume  of 
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the  same  element.  The  product  of  the  volume  density  of 
any  element  of  the  shell,  into  the  thickness  of  the  shell  in 
the  neighbourhood  of  that  element,  is  called  "surface  den 
sity."  We  see  from  the  above  that,  if  the  surface  density 
be  expressed  by  the  series 


the  potentials  at  an  internal  and  an  external  point  will  seve 
rally  be 


, 

•  •  •   "T 


o          l      ~  3  •  •  •          f.  •          i  i+  • 

3     r         or  2t  +  1    r 

This  variation  in  surface  density  may  be  obtained  either 
by  combining  a  variable  volume  density  with  an  uniform 
thickness,  as  we  have  supposed,  or  by  combining  a  variable 
thickness  with  a  uniform  volume  density,  or  by  varying  both 
thickness  and  density. 

8.  We  have  seen,  in  Chap.  II.,  that  any  positive  integral 
power  of  p,  and  therefore  of  course  any  rational  integral 
function  of  JJL,  may  be  expressed  by  a  finite  series  of  zonal 
harmonics.  It  follows,  therefore,  that  we  can  determine  the 
potential  of  a  spherical  shell,  whose  density  is  any  rational 
integral  function  of  yu,. 

Suppose,  for  instance,  we  have  a  shell  whose  density 
varies  as  the  square  of  the  distance  from  a  diametral  plane. 
Taking  this  plane  as  that  of  xy,  the  density  may  be  ex- 

z2 
pressed  by  p//,2,  or  p  ^  .     We  have  seen  (Chap.  II.  Art.  20) 

that 


Hence,  by  the  result  of  the  last  Article,  the  potential 
will  be 

V  /I     2  P  T*\ 
p  -  -  (  j-  +  -  -^-  )  at  an  internal  point, 

O  \0         0      0     J 
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U  (I     2Ptf\ 
P  3  (r  +  5  ~p~J  at  an  external  point  ; 

or,  since  P/  =  *£l±  r2  =  *?^,  we  obtain 

Z7  /I      1  3s2  -  r2\ 
P  3~  i  6  +  5  —  63  —  y        *ke  Potential  at  an  internal  point, 

U(l     V  /    1      3s'\)  „ 
^3  1r      5"  V    p  "*"  "p"  J  )  *  at  an  external  point. 

9.  As  an  example  of  the  case  in  which  the  density  is  re 
presented  by  an  infinite  series  of  zonal  harmonics,  suppose  we 
wish  to  investigate  the  potential  of  a  spherical  shell,  whose 
density  varies  as  the  distance  from  a  diameter.  Taking  this 
diameter  as  the  axis  of  z,  the  density  will  be  represented 
by  p  sin  6,  or  p  (l-p?)*.  We  have  investigated  in  Chap.  u. 
Art.  21,  the  expansion  of  sin  6  in  an  infinite  series  of  zonal 
harmonics.  Employing  this  expansion,  we  shall  obtain  for 
the  potential 

rl_         1.3...fr-l)     1.3...  (1-8)  p  r*        \ 

'  •i~*"   ' 


_ 

2  b  (2    °     16 


or 


D__  ^ 

(2r     16    V       '    2.4...t(t+2)'2A..(i-2)i    •  Vi+1" 
according  as  the  attracted  point  is  internal  or  external  to  the 
spherical  shell,  i  being  any  even  integer.     All  these  expres 
sions  may  be  obtained  in  terms  of  surface  density,  by  writing, 
instead  of  p  U,  4nrc2cr. 

10.  We  may  next  proceed  to  shew  how  the  potential  of 
a  spherical  shell  of  finite  thickness,  whose  density  is  any  solid 
zonal  harmonic,  may  be  determined.  Suppose,  for  instance, 
that  we  have  a  shell  of  external  radius  a,  and  internal  radius 
a,  whose  density,  at  the  distance  c  from  the  centre,  is 

~  Ptc\  h  being  any  line  of  constant  length. 

Dividing  the  sphere  into  concentric  thin  spherical  shells, 
of  thickness  dc,  the  potential  of  any  one  of  these  shells,  of 
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radius  c,  at  an  internal  point  distant  r  from  the  centre  will 
be  obtained  by  writing  c  for  b,  ^  for  C,  4mc2dc  for  U,  in 
the  first  result  of  Art.  6.  This  gives 


p 

'  r 


47T 


To  obtain  the  potential  of  the  whole  shell,  we  must  inte 
grate  this  expression,  with  respect  to  c,  between  the  limits 
a  and  a.  This  gives 


Again,  the  potential  of  the  shell  of  radius  c,  at  an  external 
point,  will  be 


Integrating  as  before,  we  obtain  for  the  potential  of  the 
whole  shell, 

_  47T  p       (a*i+5-a^ 

t2t  +  l)(2M-3)tf    f         rm 

Suppose  now  that  we  wish  to  obtain  the  potential  of  the 
whole  shell  at  a  point  forming  a  part  of  its  mass,  distant  r 
from  the  centre.  We  shall  obtain  this  by  considering  sepa 
rately  the  two  shells  into  which  it  may  be  divided,  the 
external  radius  of  the  one,  and  the  internal  radius  of  the 
other,  being  each  r.  Writing  r  for  a,  in  the  first  of  the  fore 
going  results,  we  obtain 


And  writing  r  for  a  in  the  other  result,  we  obtain 


Adding  these,   we  get  for   the  potential  of  the  whole 
sphere 

47r    pP,(a^-^  {       r*»*  -  a™  \ 
+l  V  1     2      T  +(2»  +  3)rmJ' 
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It  is  hardly  necessary  to  observe  that  the  corresponding 
results  for  a  solid  sphere  may  be  obtained  from  the  foregoing, 
by  putting  a  =  0. 


If  the  density,  instead  of  being  j^c1,  be  j^P^,  similar 

reasoning  will  give  us,  for  the  potential  of  the  thin  shell  of 
radius  c  and  thickness  do  at  an  internal  and  external  point 
respectively, 

^-'P-Pr'c^dc      nd      47r     p   PC"m+*<1r 

2i  + 1  km ±il         dC}  a  j  2T+1 7?  Pi  7+r  ^ 

And,  integrating  as  before,  we  obtain  for  the  potential  of 
the  whole  shell, 

;p-a"'"it8-a"""t!)  '•' at  an  i 


__       '»>i+i+3 

~i+l —    -  at  an  external  point. 


And,  at  a  point  forming  a  part  of  the  mass, 

k-t+2          ft,m-if2  wi+i+3  '  m+i+3     -f 

~~  /  •          I  —  (L 


11.  Suppose,  for  example,  that  we  wish  to  determine,  in 
each  of  the  three  cases,  the  potential  of  a  spherical  shell 
whose  external  and  internal  radii  are  a,  a,  respectively,  and 
whose  density  varies  as  the  square  of  the  distance  from  a 
diametral  plane. 

Taking  this  plane  as  that  of  xy,  the  density  may  be  ex 
pressed  by  j^2s2,  or  jrjcy.  Now  fS=~  2J"  .  Hence  the 
density  of  this  sphere  may  be  expressed  as 


The  several  potentials  due  to  the  former  term  will  be, 

2 
writing  2  for  i  and  multiplying  by  -  , 

o 


TO   THE   THEORY   OF  ATTRACTION.  01 

af-a"V  —  P  pal^°L_ 
1.5    h*  ^     a'r'105/t2    8~p~ 

And  for  the  latter  term,  writing  0  for  i,  and  2  for  771,  and 
multiplying  by  ^ , 

47T  p   ,  4  _    ,4      4?r  p  a5  -  a'5     4?r  p   /a4  —  r4  .  r5  -  a'5 

1  9   7/2  \  ^     1  ^    7,^         *,'        *     ~o     7,^ 

±-t   Iv  ±O    ft  I  o     II 

And,  since  P2r2  =  1-:^ —  ,  we  get  for  the  potential  at  an 
internal  point 


at  an  external  point 

p  UTT  J-a'7  47ra5-a' 


at  a  point  forming  a  part  of  the  mass 


1  2.  We  may  now  prove  that  by  means  of  an  infinite  series 
of  zonal  harmonics  we  may  express  any  function  of  ^  what 
ever,  even  &  discontinuous  function.  Suppose,  for  instance, 
that  we  wish  to  express  a  function  which  shall  be  equal  to 
A  from  ^  =  1  to  //,  =  X,  and  to  B  from  ^=\  to  yu,  =  -l. 
Consider  what  will  be  the  potential  of  a  spherical  shell, 
radius  c,  of  uniform  thickness,  whose  density  'is  equal  to  A 
for  the  part  corresponding  to  values  of  p  between  1  and  X, 
and  to  B  for  the  part  corresponding  to  values  of  'u,  between  X 
and  -  1. 

Divide  the  shell,  as  before,  into  indefinitely  narrow  strips 
by  parallel  planes,  the  distance  between  any  two  successive 
planes  being 
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We  have  then,  for  the  potential  of  such  a  sphere  at  any 
point  of  the  axis,  distant  z  from  the  centre, 

for  the  first  part  of  the  sphere 


and  for  the  latter  part 


These  are  respectively  equal  to 


<  C 


at  an  internal  point ;  and  to 


>C 


at  an  external  point. 

Now  it  follows  from  Chap.  II.  (Art.  23)  that  if  i  be  any 
positive  integer, 


rl 

whence,  since      P^  =  0,  it  follows  that 
rf-i 


TO   THE  THEORY  OF  ATTRACTION.  63 

r1  r A 

Also  I  PQdp  =  l-\  /   P0dp=*l+\. 

Hence  the  above  expressions  severally  become  : 
For  the  potential  at  an  internal  point  on  the  axis 

rjp  [^  (i  -  x) + .B  (i  +  x)  - ^?  {p2 w  _PO(X) j  ^ 


and  for  the  potential  at  an  external  point  on  the  axis 


Hence  the  potentials  at  a  point  situated  anywhere  are 
respectively 


at  an  internal  point; 
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and 

27rc23c  l[A  (1  -  X)  +  B(I  +  X)}  -^ 
A-B 


P    (Mc 

~P<-l(x)}"m~" 


at  an  external  point. 

Now,  if  we  inquire  what  will  be  the  potential  for   the 
following  distribution  of  density, 

-X)  +  5(1  +  X)  -  (A-B)(P.(\]  -  P0(X)}P» 


we  see  by  Art.  6  that  it  will  be  exactly  the  same,  both^at 
an  internal  and  for  an  external  point,  as  that  above  in 
vestigated  for  the  shell  made  up  of  two  parts,  whose  densities 
are  A  and  B  respectively, 

But  it  is  known  that  there  is  one,  and  only  one,  dis 
tribution  of  attracting  matter  over  a  given  surface,  which 
will  produce  a  specified  potential  at  every  point,  both  ex 
ternal  and  internal.  Hence  the  above  expression  must 
represent  exactly  the  same  distribution  of  density.  That  is, 
writing  the  above  series  in  a  slightly  different  form,  the 
expression 
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is  equal  to  A,  for  all  values  of  ft  from  1  to  X,  and  to  B  for  all 
values  of  //,  from  X  to  —  1. 

13.  By  a  similar  process,  any  other  discontinuous  function, 
whose  values  are  given  for  all  values  of  //<  from  1  to  —  1,  may 
be  expressed.  Suppose,  for  instance,  we  wish  to  express  a 
function  which  is  equal  to  A  from  fju  =  1  -to  /JL  =  \,  to  B  from 
fjb  =  \  to  ft  —  X2,  and  to  G  from  p  =  X2  to  ./-t  =  —  1.  This  will 
be  obtained  by  adding  the  two  series 

-[\+  {P,W  -P. 


For  the  former  is  equal  to  A  —  B  from  JJL  =  1  to  /J>  =  \, 
and  to  0  from  //,  =  X4  to  //,  =  —  1  ;  and  the  latter  is  equal  to 
B  from  p  =  1  to  p  =  \,  and  to  C  from  ^  =  X2  to  /A  =—  1. 

By  supposing  ^  and  C  each  =  0,  and  B  =  1,  we  deduce  a 
series  which  is  equal  to  1  for  all  values  of  /z  from  p  —  \  to 
/u.  =  X2,  and  zero  for  all  other  values.  This  will  be 


-  PO(\)  - 


This  may  be  verified  by  direct  investigation  of  the 
potential  of  the  portion  of  a  homogeneous  spherical  shell, 
of  density  unity,  comprised  between  two  parallel  planes, 
distant  respectively  c\  aad  cX2  from  the  centre  of  the 
spherical  shell. 

14  In  the  case  in  which  \  and  X2  are  indefinitely  nearly 
equal  to  each  other,  let  X2  =  X,  and  Xj  ==  X  +  cZX.  We  then 
have,  ultimately, 


F.  H. 
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Hence  Pi+l(\)  -  Pi+l(\)  -  ^(XJ  -  P 


Hence  the  series 

£{i+3p;(x)p>). 


is  equal  to  1  when  //,  =  X  (or,  more  strictly,  when  JJL  has  any 
value  from  X  to  X  +  dX)  and  is  equal  to  0  for  all  other  values 
of  //,. 

We  hence  infer  that 


is  infinite  when  /i  =  X,  and  zero  for  all  other  values  of  fi. 
15.     Representing  the  series 


by  <£(X)  for  the  moment,  we  see  that  p<£(X)c?X  is  equal  to  p 
when  fj<  =  \  and  to  zero  for  all  other  values.  Hence  the 
expression 


is  equal  to  p^  when  //.  =  Xa,  to  p^  when  /^  =  X2...  Supposing 
now  that  Xlt  X2...  are  a  series  of  values  varying  continuously 
from  1  to  —  1,  we  see  that  this  expression  becomes 


f 
*  - 


p  being  any  function   of  X,  continuous   or  discontinuous. 
Hence,  writing  $(X)  at  length,  we  see  that 


is  equal,  for  all  yalues  of  p,  from  —  1  to  +  1,  to  the  same 
function  of  /LA  that  p  is  of  X. 
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16.     The  same  conclusion  may  be  arrived  at  as  follows  : 

The  potential  of  a  spherical  shell,  whose  density  is  p, 
and  volume  U,  at  any  point  on  the  axis  of  z,  is 

pd\ 


which  is  equal  to  ^-  j  I    pd\  +  -  1   pP1(\)d\+  ... 
^c  U  -i  cJ-i 


for  an  internal  point, 

and  to    f  {-  f  pd\  +  4  I'  pP^dX  +  .  .  . 
*  (z"  -i          z  J  _i 


for  an  external  point. 

It  hence  follows  that  the   potential,  at  a  point  situated 
anywhere,  is 


for  an  internal  point, 
andto 


for  an  external  point. 

And  these  expressions  are  respectively  equal  to  those 
for  the  potentials,  at  an  internal  and  external  point  re 
spectively,  for  matter  distributed  according  to  the  following 
law  of  density  : 


5—2 
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It  will  be  observed,  in  applying  this  formula,  that  if  p  be 
a  discontinuous  function  of  X,  each  of  the  expressions  of  the 

form  I   pP^dX  will  be  the  sum  of  the  results  of  a  series  of 
-i 


-i 

integrations,  each  integration  being  taken  through  a  series  of 
values  of  X,  for  which  p  varies  continuously. 


CHAPTER  IV. 

SPHEEICAL  HAKMONICS  IN  GENEEAL.  TESSEEAL  AND  SEC- 
TOEIAL  HAEMONICS.  ZONAL  HAEMONICS  WITH  THEIR 
AXIS  IN  ANY  POSITION.  POTENTIAL  OF  A  SOLID  NEAELY 
SPHEEICAL  IN  FOEM. 

1.  WE  have  hitherto  discussed  those  solutions  of  the 
equation  V2F=0  which  are  symmetrical  about  the  axis  of  z, 
or  in  other  words,  those  solutions  of  the  equivalent  equation  in 
polar  co-ordinates  which  are  independent  of  </>.  We  propose, 
in  the  present  Chapter,  to  consider  the  forms  of  spherical 
harmonics  in  general,  understanding  by  a  Solid  Spherical 
Harmonic  of  the  iih  degree  a  rational  integral  homogeneous 
function  of  x,  y,  z,  of  the  iih  degree  which  satisfies  the  equa 
tion  V2F=0,  and  by  a  Surface  Spherical  Harmonic  of  the 
ith  degree  the  quotient  obtained  by  dividing  a  Solid  Sphe 
rical  Harmonic  by  (#2  +y*+  z*)*.  Such  an  expression,  as  we 
see  by  writing  x  =  r  sin  Q  cos  <£,  y  —  r  sin  6  sin  <£,  z  —  r  cos  6, 
will  be  of  the  Ith  degree  in  sin  6  cos  <£,  sin  0  sin  </>,  cos  0;  and 
will  satisfy  the  differential  equation  in  Yiy 


sin  6 

or,  writing  p  for  cos  &, 
d 


It  will  be  convenient,  before  proceeding  to  investigate  the 
algebraical  forms  of  these  expressions,  to  discuss  some  of 
their  simpler  physical  properties. 

2.  We  will  then  proceed  to  shew  how  spherical  har 
monics  may  be  employed  to  determine  the  potential,  and 
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consequently  the  attraction,  of  a  spherical  shell  of  indefinitely 
small  thickness. 

We  will  first  establish  an  important  theorem,  connecting 
the  potential  of  such  a  shell  on  an  external  point  with  that 
on  a  corresponding  internal  point.  The  theorem  is  as  follows : 

If  0  be  the  centre  of  such  a  sliell,  c  its  radius,  P  any  in 
ternal  point,  P'  an  external  point,  so  situated  that  P'  lies  on 
OP  produced,  and  that  OP .  OP'  =  c2,  and  if  OP  =  r,  OP'  =  r', 
then  the  potential  of  the  shell  at  P  is  to  its  potential  at  P' 
as  c  to  r,  or  (which  is  the  same  thing)  as  r'  to  c. 

For,  let  A  be  the  point  where  OP'  meets  the  surface  of 
the  sphere,  Q  any  other  point  of  its  surface.  Then,  by  a 
known  geometrical  theorem, 

QP:  QF  ::AP:AP'  ::c-r:  r'-c. 

,  c  —  r       cr  —  ?'2      cr  —  r2     r      c 

And          -, =  — =  —, =  -  =  -7 , 

r  —  c     rr  —  cr     c  —  cr     c      r 

Again,  considering  the  element  of  the  shell  in  the  im 
mediate  neighbourhood  of  Q,  its  potential  at  P  is  to  its 
potential  at  P'  as  QP  is  to  QP,  that  is,  as  c  to  r,  or  (which 
is  the  same  thing)  as  r'  to  c,  which  ratio,  being  independent 
of  the  position  of  Q,  must  be  true  for  every  element  of  the 
spherical  shell,  and  therefore  for  the  whole  shell.  Hence 
the  proposition  is  proved. 

3.     Now,  suppose  the  law  of  density  of  the  shell  to  be 

r1 

such  that  its  potential  at  any  internal  point  is  F  (/JL,  <£)  —t . 

c 

r1 
Then  F(jjb,  </>)  —  must  be  a  solid  harmonic  of  the  degree  i. 

Hence  F  (JJL,  </>)  must  be  a  surface  harmonic  of  the  degree  i. 
Let  us  represent  it  by  Y{. 

By  the  proposition  just  proved,  the  potential  at  any 
external  point,  distant  r  from  the  centre,  must  be 
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Hence,  the  component  of  the  attraction  of  the  sphere  on 
the  internal  point  measured  in  the  direction  from  the  point 
inwards,  i.  e.  towards  the  centre  of  the  sphere,  is 


And  the  component  in  the  same  direction  of  the  attraction 
on  the  external  point,  measured  inwards,  is 


Now  suppose  the  two  points  to  lie  on  the  same  line 
passing  through  the  centre  of  the  sphere,  and  to  be  both 
indefinitely  close  to  the  surface  of  the  sphere,  so  that  r  and  r 
are  each  indefinitely  nearly  equal  to  c. 

And  the  attraction  on  the  external  point  exceeds  the 
attraction  on  the  internal  point  by  - 


Now,  supposing  the  shell  to  be  divided  into  two  parts, 
by  a  plane  passing  through  the  internal  point  perpendicular 
to  the  line  joining  it  with  the  centre,  we  see  that  the  at 
traction  of  the  larger  part  of  the  shell  on  the  two  points  will 
be  ultimately  the  same,  while  the  component  attractions  of 
the  smaller  portions,  in  the  direction  above  considered,  will 
be  equal  in  magnitude  and  opposite  in  direction.  Hence  the 

difference  between  these  components,  viz.  (Zi  +  1)  —  -  ,  will  be 

c 

equal  to  twice  the  component  attraction  of  the  smaller 
portion  in  the  direction  of  the  line  joining  the  two  points. 
But  if  p€  be  the  density  of  the  shell,  So  its  thickness,  this 
component  attraction  is  ^irp^c. 

Hence  Kl% 
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And,  if  ai  be  the  corresponding  surface  density, 

~v  "T"    •*•      "T7* 

O";  =  'A -if 

4f* 

It  hence  follows  that  if  the  potential  of  a  spherical  shell, 
of  indefinitely  small  thickness,  be  a  surface  harmonic,  its 
potential  at  any  internal  point  will  be  proportional  to  the 
corresponding  solid  harmonic  of  positive  degree,  and  its  po 
tential  at  any  external  point  will  be  proportional  to  the 
corresponding  solid  harmonic  of  negative  degree. 

That  is,  the  proposition  proved  for  zonal  harmonics  in 
Chap.  in.  Art.  6,  is  now  extended  to  spherical  harmonics  in 
general. 

4.  The  spherical  ha'rmonic  of  the  degree  i  will  involve 
2i  + 1  arbitrary  constants. 

For  the  solid  spherical  harmonic,  /F.,  being  a  rational 
integral  function  of  x,  y,  z  of  the  Ith  degree,  will  consist  of 

;he    expression   V2F,   being  a 


t 
rational  integral  function  of  x,  y,  z  of  the  degree  i  —  2,  will 

consist  of  ^  ~  '  ^  terms  ;  and  the  condition  that  it  must  be 
=  0  for  all  values  of  x,  y,  z,  will  give  rise  to  - — ^ —  relations 
among  the  ^+  -£ +  '  coefficients  of  these  terms,  leaving 

ff±j)fr'  +  2)  -fr'"1)1,  or  2i  + 1,  independent  coefficients. 
2  Z 

5.  We  proceed  to  shew  how  the  spherical  harmonic  of  the 
degree  i  may  be  arranged  in  a  series  of  terms,  each  of  which 
may  be  deduced  by  differentiation  from  the  Zonal  Harmonic 
symmetrical  about  the  axis  of  z.  The  solid  zonal  harmonic, 
which,  in  accordance  with  the  notation  already  employed,  is 
represented  by  r'P.  (M),  is  a  function  of  z  and  r  of  the  degree  i, 

d*V     d*V     d*V 
satisfying  the  equation  V2  V=  0,  or  -pr  +  ^— „-  4-  -p-  =  0. 

Now,  if  we  denote  this  expression  by  Pt  (z),  we  see  that 
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since  it  is  a  function  of  z  and  r,  it  is  a  function  of  the  dis 
tance  (z)  from  a  certain  plane  passing  through  the  origin,  and 
of  the  distance  (r)  from  the  origin.  Further,  if  we  write  for  z 
the  distance  from  any  other  plane  passing  through  the  origin, 


leaving  r  unaltered,  the  equation  ^-j 
continue  to  be  satisfied. 

Now  z+}a(x-{-J—~Ly\  a.  teing  any  quantity  whatever, 
represents  the  distance  from  a  certain  plane  passing  through 
the  origin,  since  in  this  expression,  the  sum  of  the  squares 
of  the  coefficients  of  z\  x,  y  is.  equal  to  unity.  Hence 
—  !  is  a  solid  zonal  harmonic  of  the 


degree  i,  its   axis  being  the  imaginary  line  -  =  —  j=  =  z. 

Therefore  the  equation 

d'V    d'V    d*V_ 
dx*  +  dy*  +  dz2  ~ 

is  satisfied  by  V=  Pt  {z  +  a.  (x  4-  V—  1  y)}>  that  is,  expanding 
by  Taylor's  Theorem,  it  is  satisfied  by' 


_L       .      **l       «     *    «         O  t-t  rV 

/or  a^  values  of  a. 

Hence,  since  the  equation  in  V  is  linear,  it  follows  that 
it  is  satisfied  by  each  term  separately,  or  that,  besides  Pi  (z) 
itself,  each  of  the  i  expressions, 


satisfies  the  equation 

By  similar  reasoning  we  may  shew  that  each  of  the  i  ex 
pressions, 


satisfies  the  same  equation. 
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Now  each  of  the  2i  solutions,  thus  obtained,  is  imaginary. 
But  the  sum  of  any  two  or  more  of  them,  or  the  result 
obtained  by  multiplying  any  two  or  more  by  any  arbitrary 
quantities,  and  adding  the  results  together,  will  also  be  a 
solution  of  the  equation.  Hence,  adding  each  term  of  the 
first  series  to  the  corresponding  term  of  the  second,  we  ob 
tain  a  series  of  i  real  solutions  of  the  equation.  Another 
such  series  may  be  obtained  by  subtracting  each  term  of  the 
second  series  from  the  corresponding  term  of  the  first,  and 
dividing  by  V—  1.  We  have  thus  obtained  (including  the 
original  term  Pt(z))  a  series  of  2i+l  independent  solutions 
of  the  given  equation,  which  will  be  the  2i  +  1  independent 
solid  harmonics  of  the  degree  i. 

6.  We  may  deduce  the  surface  harmonics  from  these  by 
writing  r  sin  0  cos  <£  for  x,  r  sin  0  sin  (/>  for  y,  r  cos  6  for  z, 
and  dividing  by  r*.  Then,  putting  cos  9  =  /-t,  and  observing 


that  Pt  (z}  =  W  (/,),       J-      =  rl  ...  we  obtain  the  fol- 

dz  a/A 

lowing  series  of  2i  +  1  solutions  : 

3  Mi 


dp 

n20- 

Expressions  of  the  form 


.    .    a-  .    OJL    .   .a  •     .,    •   ia 

sm  6  sin  6  —  —-^  ,  sin  26  sm2  6  —  7  2  -  ,  ...  sin  16  sm*  0 
dp  dyS 


or  tfsmo^sm-0^-^, 

dp* 

or  their  equivalents, 
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(0  and  S  denoting  any  quantities  independent  of  6  and  <j>) 
are  called  Tesseral  Surface  Harmonics  of  the  degree  i  and 
order  cr.  The  particular  forms  assumed  by  them  when 
<r  =  i  are  called  Sectorial  Surface  Harmonics  of  the  degree  i. 

It  will  be  observed  that,  since  —  -  j~^-  is  a  numerical  constant, 
Sectorial  Harmonics  only  involve  6  in  the  form 

sin'0,  or  (I-//)*. 

The  product  obtained  by  multiplying  a  Tesseral  or 
Sectorial  Surface  Harmonic  of  the  degree  i  by  r*  (that  is, 
the  expression  directly  obtained  in  Art.  5)  is  called  a  Tesseral 
or  Sectorial  Solid  Harmonic  of  the  degree  i. 

7.  We  shall  denote  the  factor  of  a  Tesseral  or  Sectorial 
Harmonic  which  is  a  function  of  0,  that  is  sm'd  —  i  ^  ,  or 

(1  —  A*-2)2  —   '^     ,  by  the  symbol  T^\  or,  when  it  is  necessary 

to  particularize  the  quantity  of  which  it  is  a  function,  by 
or  Z*>cos0. 


It  will  be  convenient,  for  the  purpose  of  comparison  with 
the  forms  of  Tesseral  Harmonics  given  in  the  Mecanique 
Celeste,  and  elsewhere,  to  obtain  TW  in  a  completely  de 
veloped  form. 

1  di(u?—'\.}i 

Now,  since  Pt(p)  =  £.  ^   ^   g  —  ^.  —  ^  —  i.  ,  we  see  that 

dvPt(fi)=  1  df+'df-iy 

d*        2M.2.  3,..*      di+0> 
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t-  4)  (2t-  5)...(i  -  <r  - 


2.  4.  (2*  -1)  (2*  -3) 

And  therefore 


2.4(2i-l)(2*-3) 


The  form  given  by  Laplace  for  a  Tesseral  Surface  Har 
monic  of  the  degree  i  and  order  cr  is  (see  Hecanique  Celeste, 
Liv.  3,  Chap.  2,  pp.  40—47) 


^4  being  a  quantity  independent  of  6  and  </>.  The  factor  of 
this,  involving  /*,  is  denoted  by  Thomson  and  Tait  (Natural 
Philosophy,  Vol.  1,  p.  149)  by  trie  symbol  ©/').  Thomson 
and  Tait  also  employ  a  symbol  $•/*>,  adopted  by  Maxwell  in 
his  Treatise  on  Electricity  and  Magnetism,  Vol.  1,  p.  164, 
which  is  equal  to 


ft        . 
V 


or  2" 
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Heine  represents  the  expression 

-'-  (^-g-)(^-°-1)     i-cr-2 


2.4.(2f-l)(2t-3) 

<r 

or  (—  I)2  S<r\  by  the  symbol  PO-'(/A),  and  calls  these  expres 
sions  by  the  name  Zugeordnete  Functionen  Erster  Art  (Hand- 
buch  der  Kugelfunctionen,  pp.  117,  118)  which  Todhunter 
translates  by  the  term  "Associated  Functions  of  the  First 
Kind,"  which  we  shall  adopt. 

Heine  also  represents  the  series 

•  0     (fr-q-)O'-q--i)  ^-a 

**'  2(2i-l) 

+  —  2  .V(2]-  1)  (2?-  3)'  ""  °"  ~      /A*""  4> 

by  the  symbol  ^(/*),  (p.  117). 

The  several  expressions,  Tf\  ®|^,  ^J^,  Pi,  3J&,  are  con 
nected  together  as  follows : 

2*.  1.2.  3.. .t 


8.     It  has  been  already  remarked  that  the  roots  of  the 
equation  Pt  =  0  are  all  real.    It  follows  also  that  those  of  the 

-IT)  /72  P 

equations  -ri=  0,  -j-^  =  0.,.  are  real  also.     Hence  we  may 

arrive  at  the  following  conclusions,  concerning  the  curves, 
traced  on  a  sphere,  which  result  from  our  putting  any  one 
of  these  series  of  spherical  harmonics  =  0. 

By  putting  a  zonal  harmonic =0,  we  obtain /small  circles, 
whose  planes  are  parallel  to  one  another,  perpendicular  to 
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the  axis  of  the  zonal  harmonic,  and  symmetrically  situated 
with  respect  to  the  diametral  plane,  perpendicular  to  this 
axis.  If  i  be  an  odd  number  this  diametral  plane  itself 
becomes  one  of  the  series. 

By  putting  the  tesseral  harmonic  of  the  order  <r=0,  we 
obtain  i  —  cr  small  circles,  situated  as  before,  and  cr  great 
circles,  determined  by  the  equation  cos  <r<j>  =  0,  or  sin  a<j)  =  0, 
as  the  case  may  be,  their  planes  all  intersecting  in  the  axis 
of  the  system  of  harmonics,  the  angle  between  the  planes  of 

7T 

any  two  consecutive  great  circles  being  -  . 


By  putting  the  sectorial  harmonic  =  0,  we  obtain  i 
great  circles,  whose  planes  all  intersect  in  the  axis  of  the 
system,  the  angle  between  any  two  consecutive  planes  being 

7T 


9.  The  tesseral  harmonic  may  be  regarded  from  another 
point  of  view.  Suppose  it  is  required  to  determine  a  solid 
harmonic  of  the  degree  i,  and  of  the  form  Yf*,  such,  that  Yi 
shall  be  the  product  of  a  function  of  p,  and  of  a  function  of  </>, 
which  functions  we  will  denote  by  the  symbols  -M.,  <£0  respec 
tively.  The  differential  equation,  to  which  this  will  lead,  is 


Now  this  will  be  satisfied,  if  we  make  Mt  and  <!>,  satisfy 
the  following  two  equations  : 


The  latter  equation  gives 

<£>,  =  C  cos  <?(}>  +  C'  sin  cr0. 

And,  taking  a-  as  an  integer,  positive  or  negative,  the 
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°Y  d  \l'+<r 
former  is  satisfied  by  M^  Tp\  i.e.  (l^W(ft)      (1  -/*% 

as  we  proceed  to  prove. 
We  know  that 


Differentiate  cr  times,  and  we  get 


whence,  by  Leibnitz's  Theorem, 


or 


and,  multiplying  by  (1  —  /*2 


=  0  -  (I)- 


Now,  putting 
we  get 
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And       i  (i  +  1)  r(<'>  =  i  (i  +  1)  (1  -  tf  . 


Hence  the  equation  above  given  for  Mt  is  satisfied  by 
Mi  =  T^\  and  the  equation  in  Yi  is  satisfied  by 

Yt  =  CTf*  cos  o-^>  +  (7'  Tf»  sin  o-0. 

10.  In  Chap.  n.  Art.  10  we  have  established  the  fundamental 
property  of  Zonal  Harmonics,  that  if  i  and  m  be  two  unequal 

positive  integers,  J     PtPmdfA  =  Q.     This  is  a  particular  case 

of  the  general  theorem  that  if  Yif   Ym  be  two  surface  har 
monics  of  the  degrees  i  and  m  respectively, 


P  [*"YY 
J-iJ,    '  ' 
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For,  let  Vt,  Vm  be  the  corresponding  solid  harmonics,  so 
that  F,=r*F,,  Vm  =  rmYm.  Then,  by  the  fundamental  pro 
perty  of  potential  functions,  we  have  at  every  point  at  which 
no  attracting  matter  is  situated, 


. 

da?      dtf       dz*       '      dx*  "      dy*  "  dz* 
and  therefore 


^         df         dzz  m\dx*       dy*       dzz     " 

or,  in  accordance  with  our  notation,  V^  Vm  —  FTOy2  V.  =  0. 


Now,  integrate  this  expression  throughout  the  whole 
space  comprised  within  a  sphere  whose  centre  is  the  origin 
and  radius  a,  a  being  so  chosen  that  this  sphere  contains  no 
attracting  matter.  We  then  have 

VmV*V,)  dx  Ay  dz  =  0. 

But  also,  when  the  integration  extends  over  all  space 
comprised  within  any  closed  surface,  we  have 


dS  denoting  an  element  of  the  bounding  surface,  and  -7- 


<MUO, 

-7- 


differentiation  in  the  direction  of  the  normal  at  any  point. 

Now,  in  the  present  case,  the  bounding  surface  being  a 
sphere  of  radius  a,  and  Vif  Vm  homogeneous  functions  of  the 
degrees  i,  m,  respectively, 


and,  the  integration  being  extended  all  over  the  surface  of 
the  sphere,  the  limits  of  p  are  —  1  and  1,  those  of  <f>,  0  and  STT. 
Hence 


F.H. 
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whence,  if  m  -  i  be  not  =  0, 

/-JT  '^*Mty  =  0. 

ri   r27r 

The  value  of  /     I     Y?dp<ty  will  be  investigated  here- 
•*  —V  o 

after. 

11.  We  may  hence  prove  that  if  a  function  of  /j,  and  <f> 
can  be  developed  in  a  series  of  surface  harmonics,  such  de 
velopment  is  possible  in  only  one  way. 

For  suppose,  if  possible,  that  there  are  two  such  develop 
ments,  so  that 


and  also 


Then  subtracting,  we  have 

0=  FO-  F;+  y;-  F;+  ...  +  F,-  Y:  +  ...  identically. 

^  Now,  each  of  the  expressions  Y9  -  F0',  Ft  -  T/...  F.  -  F/ 
being  the  difference  of  two  surface  harmonics  of  the  degree 
0,  1,  ...i...  is  itself  a  surface  harmonic  of  the  degree 
0,  1,  ...i  —  Denote  these  expressions  for  shortness  by 
ZQ,  Zl  .  .  .  Zi  .  .  .  so  that 

0==Z0-\-Zl+...+Z.  +  ...  identically. 

Then,  multiplying   by  Z4   and  integrating  all  over  the 
surface  of  the  sphere,  we  have 


-i   o 


That  is,  the  sum  of  an  infinite  number  of  essentially 
positive  quantities  is  =0.  This  can  only  take  place  when 
each  of  the  quantities  is  separately  =  0.  Hence  Z.  is  identi 
cally  =  0,  or  F/  =  F0  and  therefore  the  two  developments 
are  identical. 

We  have  not  assumed  here  that  such  a  development  is 
always  possible.  That  it  is  so,  will  be  shewn  hereafter. 
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12.  By  referring  to  the  expression  for  a  surface  har 
monic  given  in  Art.  4,  we  see  that  each  of  the  Tesseral  and 
Sectorial  Harmonics  involves  (1  —  //,*)%  or  some  power  of 
(1  —  //,2)2,  as  a  factor,  and  therefore  is  equal  to  0  when  //,  =  ±  1. 
From  this  it  follows  that  when  /*  =  ±  1,  the  value  of  the 
Surface  Harmonic  is  independent  of  </>,  or  that  if  Y  (/*,  <j>)  repre 
sent  a  general  surface  harmonic,  Y  (±  1,  <f>)  is  independent  of 
<£>,  and  may  therefore  be  written  as  Y  (±  1).  Or  F(l)  is  the 
value  of  F(X  <W  a^  the  pole  of  the  zonal  harmonic  P4(/<t), 
F(— 1)  at  the  other  extremity  of  the  axis  of  Pi  (/A). 

We  may  now  prove  that 

Jo 
For,  recurring  to  the  fundamental  equation, 

Now,  if  we  integrate  this  equation  with  respect  to  <£, 
between  the  limits  0  and  2?r,  we  see  that,  since 


and  the  value  of  Yt  only  involves  <£  under  the  form  of  cosines 
or  sines  of  d>  and  its  multiples,  and  therefore  the  values  of 

7  y 

-~T~  are  the  same  at  both  limits,  it  follows  that 


Hence 


rar 
Hence    I     y,^  is  a  function  of  /*  which  satisfies  the 

J  o 

fundamental  equation  for  a  zonal  harmonic,  and  we  therefore 
have 

6—2 
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f 


0 

C  being  a  constant,  as  yet  unknown. 

To  determine  C,  put  ^=1,  then  by  the  remark  just  made, 
Yi  becomes   Yt(T),  and  is  independent  of  <£.     Hence,  when 

A*  =  1,  f'j  Yid$  =  27r  F.  (1).     Also  Pt(fj)  =  1.     We  have  there 
fore 


r. 
It  follows  from  this  that 


13.     We  may  now  enquire  what  will  be  the  value  of 


F.,  Zi  being  two  general  surface  harmonics  of  the  degree  i. 
Suppose  each  to  be  arranged  in  a  series  consisting  of  the 
zonal  harmonic  P.  whose  axis  is  the  axis  of  z,  and  the  system 
of  tesseral  and  sectorial  harmonics  deduced  from  it.  Let  us 
represent  them  as  follows  : 


TW  cos  2<£  +  ...  + 

+ 
^W  sin  2^  +  ...  +  ft^W  sino-^  +  ... 


cos 


+  52T/2)  sin  20  +  .  .  .  +  ^  W  sin  o-^  +  ... 


Hence  the  product  Y^.  will  consist  of  a  series  of  terms, 
in  which  <£  will  enter  under  the  form  cos  <r<£  cos  o-'ijk,  or 
cos  a(f>  sin  cr'<£.  This  expression  when  integrated  between 
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the  limits  0  and  STT  vanishes  in  all  cases,  except  when 
cr  =  a-  and  the  expression  consequently  becomes  equal  to 
cos2  cr<£,  or  sin2tcr<£.  In  these  cases  we  know  that,  cr  being  any 
positive  integer, 


/27r  rEir 

cos2  CT<£  d(p  =       si 
j  Jo 


sn  c-<     <   =  TT. 


Hence  the  question  is  reduced  to  the  determination  of  the 
value  of 


a 

l 


But,  by  the  theorem  of  Rodrigues,  proved  in  Chap.  II. 
Art.  8,  we  know  that 


Hence  5^^  may  also  be  expressed  under  the  form 


whence  it  follows  that 


2\1.2.3...i7  \i-a- 


Now,  putting  Qu,2-l)'  =  Jf  for  the  moment,   and   inte 
grating  by  parts, 

^Mj    __di+'-lM  dt-'M 
*~  ^~  dp^*-1  ^dpF^ 

rM+^Mtt-^M 
J  "d^+*-i   dp*-***-    ^ 
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- 

The  factor  —3-3  —  vanishes  at  both  limits,  hence 
- 


*M  ,          f1  di+*-lM  d*-"*1  M., 

*  ^""-1  ~^i 


by  a  repetition  of  the  same  process. 

And  by  repeating  this  process  <r  times,  we  see  that 


=  (-  1)-  (21.  1.2.3  ...i')2 


and  therefore 


/i 
..( 


=  (l    r(TV  sin  o-< 
J-iJo 

<r     2?r 


It  will  be  observed  that  this  result  does  not  hold  when 
a-  =  0,  in  which  case  we  have 


P"   \^  Y^ 

J  -iJo 


Hence 

J  -iJo 

4t7T 


In  this  case   f  *"  cos2  a<bd(j>=  I      sin2 
Jo  •'O 
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I'll,, 


14.  We  have  hitherto  considered  the  Zonal  Harmonic 
under  its  simplest  form,  that  of  a  "  Legendre's  Coefficient "  in 
which  the  axis  of  z,  i.  e.  the  line  from  which  0  is  measured,  is 
the  axis  of  the  system.  We  shall  now  proceed  to  consider  it 
under  the  more  general  form  of  a  "Laplace's  Coefficient," 
in  which  the  axis  of  the  system  of  zonal  harmonics  is  in  any 
position  whatever,  and  shall  shew  how  this  general  form  may 
be  expressed  in  terms  of  Pi  (^  and  of  the  system  of  Tesseral 
and  Sectorial  Harmonics  deduced  from  it. 

Suppose  that  &y  <£>'  are  the  angular  co-ordinates  of  the 
axis  of  the  Zonal  Harmonic,  i.e.  that  the  angle  between  this 
axis  and  the  axis  of  z  is  &,  and  that  the  plane  containing 
these  two  axes  is  inclined  to  a  fixed  plane  through  the  axis 
of  z  which  we  may  consider  as  that  of  zx,  at  the  angle  <£>'. 
In  accordance  with  the  notation  already  employed,  we  shall 
represent  cos  &  by  fjf. 

The  rectangular  equations  of  the  axis  of  this  system 
will  be 

«      =      y  z 

sin  6'  cos  <£'      sin  &  sin  <£'      cos  &  ' 

Hence  the  Solid  Zonal  Harmonic  of  which  this  is  the  axis 
is  deduced  from  the  ordinary  form  of  the  solid  zonal  har 
monic  expressed  as  a  function  of  z  and  r  by  writing,  in  place 
of  z,  x  sin  &  cos  <£'  +  y  sin  &  sin  </>'  +  z  cos  6'. 

To  deduce  the  Surface  Zonal  Harmonic,  transform  the  solid 
zonal  harmonic  to  polar  co-ordinates,  by  writing  r  sin  6  cos  <f> 
for  x,  r  sin  6  sin  <£  for  y,  r  cos  6  for  z,  and  divide  by  r*. 

The  transformation  from  the  special  to  the  general 
form  of  surface  zonal  harmonic  may  be  at  once  effected, 
by  substituting  for  /u-,  or  cos  6,  cos0cos0'+sm0sin0'cos(0— $'). 

Now,  in  order  to  develope 

Pt  {cos  0  cos  &  +  sin  6  sin  &  cos  (<£  — 
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in  the  manner  already  pointed  out,  assume 

Pt  {cos  6  cos  &  +  sin  6  sin  &  cos  (<£  —  <£')] 

cos  </>  +  £<«  sin  0)  2J 

cos  2</>  +  £<2>  sin  20)  TV  +  ... 


cos  (T(>  +     *  sn  0-        .    +  .  .  . 
+  (  <?<*>  cos  ty  +  flW  sin  *0)  ZJO 

the  letters  A,  ...  C^\  S^  ...  denoting  functions  of  pf  and 
0',  to  be  determined. 


To  determine  C^\  multiply  both  sides  of  this  equation 
by  cos  o-(f>T^  and  integrate  all  over  the  surface  of  the  sphere, 
i.  e.  between  the  limits  —  1  and  1  of  fju,  and  0  and  2?r  of  0. 
We  then  get 

I  l    [  "  Pi  {cos  e  cos  &  +  sin  0  sin  &  cos  (0  -  0')]  cos  <r$  ZJ 
J  -i./  o 


=  (7  W  (cos  <76T 

-IJO 


It  remains  to  find  the  value  of  the  left-hand  member  of 
this  equation. 

Now  cos  <7<£  2y<7)  is  a  surface  harmonic  of  the  degree  i,  and 
therefore  a  function  of  the  kind  denoted  by  Yi  in  Art.  12. 

And  we  have  shewn,  in  that  Article,  that 


that  is,  that  i/  a?!?/  surface  harmonic  of  the  degree  i  6e  multi 
plied  by  the  zonal  harmonic  of  the  same  degree,  and  the  product 
integrated  all  over  the  surface  of  the  sphere,  the  integral  is 

equal  to  ^  —  =-  into   the  value  which  the  surface  harmonic 
2i+  1 

assumes  at  the  pole  of  the  zonal  harmonic. 
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Hence 

P    I*"  P6  (cos  9  cos  &  +  sin  0  sin  &  cos  (0  -  0')}  Y,  (/*,  <£)  cfy&ety 
J  -iJo 


and  therefore 


P  [^  P4  {cos  0  cos  0'  +  sin  6  sin  0'  cos  (</>  -  <£')}  cos  v^T 


Hence 


\i  —  cr 
or  GW  =  2  ,  .         cos  vd>  Tr>  (ft ). 

\%  +  <T 


Similarly  jSfW  =  2  sin  erf  T^  (^)- 

t  ~T~  O" 

And  to  determine  A,  we  have 

I  ^P,  {cos  0  cos  0'  +  sin  0  sin  0'  cos  (<£  -  £' 

J  o 


—i  o 


or        - 
Hence,  Pf  {cos  0  cos  6'  +  sin  0  sin  0'  cos  (</>  -  <£' 


=  P<  00  P<  M  +  2        cos  (^  -  «#,') 

I  t  -r  J. 


+  2  [=  cos  2  (^  -  ^  r,»  0*')  2J»  W  + 
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+  2        cos  »•  ft  -  f  )  T/0  (/)  ZJM 


15.  We  have  already  seen  (Chap.  n.  Art.  20)  how  any 
rational  integral  function  of  p  can  be  expressed  by  a  finite 
series  of  zonal  harmonics.  We  shall  now  shew  how  any 
rational  integral  function  of  cos  6,  sin  6  cos  <£,  sin  0  sin  <£, 
can  be  expressed  by  a  finite  series  of  zonal,  tesseral,  and 
sectorial  harmonics. 

For  any  power  of  cos  </>  or  sin  $,  or  any  product  of  such 
powers,  may  be  expressed  as  the  sum  of  a  series  of  terms  of 
the  form  cos  <r(f>,  or  sin  cr^>,  the  greatest  value  of  cr  being  the 
sum  of  the  indices  of  cos  </>  and  sin  <£,  and  the  other  values 
diminishing  by  2  in  each  successive  term.  Hence  any 
rational  integral  function  of  cos  0,  sin  0  cos  0,  sin  6  sin  cf,  will 
consist  of  a  series  of  terms  of  the  form 

cosm  0  sin"  0  cos  a(j>  or  cosm  0  sin*  0  sin  o-<£, 
where  n  is  not  less  than  <r. 

If  n  be  greater  than  cr,  n  —  o-  must  be  an  even  integer.  Let 
n  —  a  =  2s,  then  writing  sin"#  under  the  form  (1  —  cos*0)*  sin0"^, 
we  reduce  cosm  0  sinn  0  cos  a(f>  to  the  sum  of  a  series  of  terms 
of  the  form  cosp  0  sin0"  0  cos  cr^>,  or,  writing  cos  0  =  p,  of  the 


<r 


form  /*p  (1  —  /A2)  *  cos  o-0. 

Similarly    cos™  0  sin"  0  sin  <n£  is  reduced  to  a  series   of 

cr 

terms  of  the  form  fjf  (1  -  ^2)  2  sin  o-(£. 


and  /AP+O>  can  be  developed  in  a  series  of  terms  of  the  form 
of  multiples  of  Pp+<r,  Pp+<r-z (Chap.  II.  Art.  17.) 

Hence  yfl  can  be  expressed  in  a  series  of  the  form 
~(A0Pp 
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AQ)  At  representing  known  numerical  constants,  and  therefore 

cr 

IJP  (1  —  //,2)a  assumes  the  form 


( A    T(<T) 

V"-0  •*•#+< 

consequently  multiplying  these  series  by  cos  a$>  or  sin  a$>,  we 
obtain  the  developments  of 

E.  — 

uP  (1  —  a2) 2  cos  crd)  and  fjP  (1  —  tt2) *  sin  crd> 

in  series  of  tesseral  harmonics. 

16.    We  will  give  two  illustrations  of  this  transformation. 

First,  suppose  it  is  required  to  express  cos2  6  sin20  sin  cf>  cos<£ 
in  a  series  of  Spherical  Harmonics. 

Here  we  have    sin  6  cos  <f>  =  -  sin  2$. 

Z 

Hence  cos2  6  sin2  6  sin  <£  cos  <£  =  -  cos2  ^  sin2  6  sin  2(£. 

& 

Comparing  this  with  cos™  0  sinn  6  sin  cr0,  we  see  that  n  is 
not  greater  than  cr. 

Hence  cos2  0  sin2  0  sin  <£  cos  cj>  —  -  $  (1  —  /^2)  sin  2(/>. 

.z 

2_      1        6?2       4 

and  /u,4  = ; 

4 

*'*  ^  =  12  V35 


105  d/42      21 
.'.  cos8  6  sin2  5  sin  <   cos  < 


92  SPHERICAL  HARMONICS  IN  GENERAL. 

Next,  let  it  be  required  to  transform  cos3  6  sin3  6  sin  0  cos20 
into  a  series  of  Spherical  Harmonics. 

Here  sin  0  cos2  0  =  ^  sin  20  cos  0  =  -  (sin  30  4-  sin  0). 
Now  cos3  6  sin3  0  sin  30  =  fj?  (1  —  /^  sin  30 


Also  cos3  6  sin3  #  sin  $  =  ^  (1  —  /A")  (1  —  /i2)4  si 


sn 


Also  (Chap.  II.  Art.  17) 


231    6     77    •     21 
Hence  cos3  6  sin3  6  sin  3</> 


1    /16<?P6  ,  ,.        ... 

120  (iTl  d     +  77  ~        sm 


p      77  c?^2      63 
2  ciP 


4 


. 
(      ^ 


f  8   dPR       2    dP4      4 
"(693  ^--- 


693    6        385    4        63    2 
^  sin0  cos20  =   JL 


693    a        770    4       63 
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17.  The  process  above  investigated  is  probably  the  most 
convenient  one  when  the  object  is  to  transform  any  finite 
algebraical  function  of  cos  0,  sin  0  cos  <£,  and  sin  0  sin  <£,  into 
a  series  of  spherical  harmonics.  For  general  forms  of  a 
function  of  p  and  <£,  however,  this  method  is  inapplicable, 
and  we  proceed  to  investigate  a  process  which  will  apply 
universally,  even  if  the  function  to  be  transformed  be  discon 
tinuous. 

We  must  first  discuss  the  following  problem. 

To  determine  the  potential  of  a  spherical  shell  whose 
surface  density  is  F '(//,,  <£),  F  denoting  any  function  whatever 
of  finite  magnitude,  at  an  external  or  internal  point. 

Let  c  be  the  radius  of  the  sphere,  r'  the  distance  of  the 
point  from  its  centre,  0',  <j>'  its  angular  co-ordinates,  V  the 
potential.  Then  fju  being  equal  to  cos  0 

y  =  I      I      • — ____________ — _ 

J  - iJ  o    \r'z—  2cr  (cos  0  cos  0'  +  sin  0  sin  &  cos  (<£  -  <£'))  +  c2] * 


The  denominator,  when  expanded  in  a  series  of  general 
zonal  harmonics,  or  Laplace's  coefficients,  becomes 

(  (ft  fl     +  . 


for  an  internal  and  an  external  point  respectively,  Pt  (p, 
being  written  for 

Pt  (cos  6  cos  0'  +  sin  0  sin  0'  cos  ($  -  0')}. 


Hence,  Fx  denoting  the  potential  at  an  internal,  F2  at  an 
external,  point, 


-1J  0 


tj  r  f  fcpt  (M> 

C  J^jJ  o 

...+J  r  rpi(/U)  ^^ 

c  ^_^o 
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=  r'  {f  JT* 


It  will  be  observed  that  the  expression  Pi  (p,  </>)  involves 
/j,  and  //  symmetrically,  and  also  <j>  and  <j>'.  Hence  it  satisfies 
the  equation 


And,  since  ^u.  and  <£  are  independent  of  fjf  and  (/>',  this 
differential  equation  will  continue  to  be  satisfied  after  Pt  has 
been  multiplied  by  any  function  of  p,  and  </>,  and  integrated 
with  respect  to  p  and  <£.  That  is,  every  expression  of  the 
form 


T 

J  -iJ  0 


is  a  Spherical  Surface  Harmonic,  or  "Laplace's  Function" 
with  respect  to  p  and  </>'  of  the  degree  i.  And  the  several 
terms  of  the  developments  of  Vt  are  solid  harmonics  of  the 
degree  0,  1,  2...i...  while  those  of  F2  are  the  corresponding 
functions  of  the  degrees  —  1,  —2,  —  3...  —  (i+  1),  ...  And 
these  are  the  expressions  for  the  potential  at  a  point  (/,  //,  <£') 
of  the  distribution  of  density  F(p!,  <£')  at  a  point  (c,  p,  <£')• 

Now,  the  expressions  for  the  potentials,  both  external 
and  internal,  given  in  the  last  Article,  are  precisely  the  same 
as  those  for  the  distribution  of  matter  whose  surface  density  is 


ff  p 

7T  ^/  _jj  Q 


or,  as  it  may  now  be  better  expressed, 
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+  3  f  f**P1  {cos0  cos0f +  sin0  sinfl'  cos  (<£-<£')  F(fit 
J  ~iJ  o 

r1  [2*  *  *  1 

J-tJo 


And,  since  there  is  only  one  distribution  of  density  which 
will  produce  a  given  potential  at  every  point  both  external 
and  internal,  it  follows  that  this  series  must  be  identical 
with  F(p!t  <£').  We  have  thus,  therefore,  investigated  the 
development  of  F(fi9  <j>)  in  a  series  of  spherical  surface 
harmonics*. 

The  only  limitation  on  the  generality  of  the  function 
F(fA,  (/>')  is  that  it  should  not  become  infinite  for  any  pair  of 
values  comprised  between  the  limits  —  1  and  1  of  //.,  and  0 
and  2-Tr  of  <p. 

18.  Ex.  To  express  cos  20'  in  a  series  of  spherical  har 
monics. 

For  this  purpose,  it  is  necessary  to  determine  the  value  of 


Now  Pt  {cos  6  cos  &  -f  sin  6  sin  &  cos  (<£  - 

=  P<  (cos  flJP^cosfl1) 

2        .      dP^cosfl)   .    /K<ZP1(cos0') 

-^     "-" 


rz* 

Now  I     cos  o-  (A  -  A')  cos  26  ^  =  0, 

^o 

for  all  values  of  cr  except  2. 

*  In  connection  with  the  subject  of  this  Article,  see  a  paper  by  Mr  G.  H. 
Darwin  in  the  Messenger  of  Mathematics  for  March,  1877. 
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fZir 

And         j    cos  2  (£  -  £')  cos  2<£  c?6  =  TT  cos  26'. 

JO 

Also 


And 


Now  when  yu-  =  1, 


And  when  /x  =  —  1, 


dfj+l  d/j? 

Hence 


n 


=  4  or  0,  as  i  is  even  or  odd ; 

J_JQ   '  djj?  ^     ™' 

=  4t7r  cos  2<£'  or  0,  as  i  is  even  or  odd ; 
.*.  cos  2^)' 

2         .    .  2  fi,  tfP^  (cos  ^) 
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1 1 

d  /-»    t  //"•-*-  o  .          •  /.-«-,  JL  O  . 

=  2  cos 


0.,  /  5.rtp 

2*(i±t: 


4      3.4.56      5.6.778 


Hence  the  potential  of  a  spherical  shell,  of  radius  c  and 
surface  density  cos  2<£',  will  be 

/       rn(2)  ,2  m(2J  '4  m(2)  '6 

8.  cos  2^'  (rfe  ?  +  ^^  ^-  +  .--{V:8  -+  . 
and 


at  an  internal  and  external  point  respectively. 

19.  We  will  now  explain  the  application  of  Spherical 
Harmonics  to  the  determination  of  the  potential  of  a  homo 
geneous  solid,  nearly  spherical  in  form.  The  following 
investigation  is  taken  from  the  Mecanique  Celeste,  Liv.  in. 
Chap.  II. 

Let  *•  be  the  radius  vector  of  such  a  solid,  and  let 

r  =  a  -f-  a  (ax  Ft  4-  a2  Y2  +  ...  +  a,  Y,  +  ...), 

a  being  a  small  quantity,  whose  square  and  higher  powers 
may  be  neglected,  at,  a.2,...«£...  lines  of  arbitrary  length,  and 
Ft,  F2,..,Y....  surface  harmonics  of  the  order  1,  2,...i...  re 
spectively. 

4 

The  volume  of  the  solid  will  be  _  ira*. 

o 


For  it  is  equal  to 


rr  rl 

J  Q  J  -iJ  o 


r1  r27r 

,  since   |     |      YidjJid(j>  =  Q, 


-i  .'  0 

for  all  values  of  i. 
F.  H. 
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Again,  if  the  centre  of  gravity  of  the  solid  be  taken  as 
origin,  ^  =  0. 

For  if  ~z  be  the  distance  of  the  centre  of  gravity  from  the 
plane  of  xy, 


4  r  r1  r27r 

,37ra32=       I          r3u,drdad<t> 
J-iJo 


a1        I     pY 

J  _  J  o 


„  4a3  a  .  a,  f    f  *(1  -  /^  cos 

o  J  -i  J  o 

4  r1    ran- 

gwa8y  =  4a8a.ai  I     I     (l-^sin 


) 


Similarly 


Now  Fj  is  an  expression  of  the  form 

AJJ,  +  B  (1  —  /j?j*  cos  </>  +  (7(1  —  yu,2)^  sin  <£, 

and  therefore  all  the  expressions  x,  y,  z  cannot-  be  equal  to  0, 
unless  at  =  0. 

We  may  therefore,  taking  the  centre  of  gravity  as  origin, 
write 

as  the  equation  of  the  bounding  surface  of  thfe  solid. 

Now  this  solid  may  be  considered  as  made  up  of  a  homo 
geneous  sphere,  radius  a,  and  of  a  shell,  whose  thickness  is 

a  (a,F,+  ...  +  a,r,  +  ...)- 

The  potential  of  this  shell,  at  least  at  points  whose  least 
distance  from  it  is  considerable  compared  with  its  thickness, 
will  be  the  same  as  that  of  a  shell  whose  thickness  is  aa,  and 
density 

ai  i"          \ 


..+  . 
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p0  being  the  density  of  the  solid.     Therefore  the  potential, 
for  any  external  point,  distant  R  from  the  centre,  will  be 


The  potential  at  any  internal  point,  distant  E  from  the 
centre,  will  be  made  up  of  the  two  portions 


vpJP  +  27TPo  (a2  -  £2)  or  2^    a«  -     - 

for  the  homogeneous  sphere, 

Ytf  aY 


for  the  shell,  and  will  therefore  be  equal  to 


20.  If  the  solid,  instead  of  being  homogeneous,  be  made 
up  of  strata  of  different  densities,  the  strata  being;  concentric, 
and  similar  to  the  bounding  surface  of  the  solid,  we  may 

x» 

deduce  an  expression  for  its  potential  as  follows.     Let  -  r  be 

a 

the  radius  vector  of  any  stratum,  p  its  density,  r  having  the 
same  value  as  in  the  last  Article,  and  p  being  a  function 
of  c  only.  Then,  Sc  being  the  mean  thickness  of  the  stratum, 
that  is  the  difference  between  the  values  of  c  for  its  inner 
and  outer  surfaces,  the  potential  of  the  stratum  at  an  ex 
ternal  point  will  be 


I          *      »    I 

^St  +  ljB1^ 

To  obtain  the  potential  of  the  whole  solid  at  an  external 
point  we  must  integrate  this  expression  with  respect  to  c, 
between  the  limits  0  and  a,  remembering  that  p  is  a  func 
tion  of  c. 

7-2 
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Again,  the  potential  of  the  stratum,  above  considered, 
at  an  internal  point  will  be 

*  '     - 


To  obtain  the  potential  of  the  whole  solid  at  an  internal 
point  we  must  integrate  the  expression  (1)  with  respect  to  c 
between  the  limits  0  and  R,  and  the  expression  (2)  with 
respect  to  c  between  the  limits  R  and  a,  remembering  in 
both  cases  that  p  is  a  function  of  c,  and  add  the  results 


together. 


|Z^v    cu  ^C 


CHAPTER  V. 

SPHERICAL  HARMONICS  OF  THE  SECOND  KIND. 

1.  WE  have  already  seen  (Chap.  n.  Art.  2)  that  the 
differential  equation  of  which  Pi  is  one  solution,  being  of 
the  second  order,  admits  of  another  solution,  viz. 


Now  if  //-  between  the  limits  of  integration  be  equal 
to  ±  1,  or  to  any  roots  of  the  equation  Pt  =  0  (all  of  which 
roots  lie  between  1  and  —  1),  the  expression  under  the 
integral  sign  becomes  infinite  between  the  limits  of  inte 
gration.  We  can  therefore  only  assign  an  intelligible 
meaning  to  this  integral,  by  supposing  /*  to  be  always  be 
tween  1  and  oo  ,  or  between  —  1  and  —  GO  .  We  will  adopt 
the  former  supposition,  and  if  we  then  put  (7=  —  1,  the 


.  -11  v   i 

expression    ^^-  --  §s(i.e,«r7~i  —  it    Wl11  be  always  P°si- 

•t  i  (L  ~  P  )  \      *i  (P  ~  L)/ 
tive.     We  may  therefore  define  the  expression 


as  the  zonal  harmonic  of  the  second  kind,  which  we  shall 
denote  by  Qiy  or  Q.  (/*),  when  it  is  necessary  to  specify  the 
variables  of  which  it  is  a  function. 

It  will  be  observed  that,  if  fj,  be  greater  than  1,  Pt  is 
always  positive.  Hence,  on  the  same  supposition,  Qt  is 
always  positive. 

Weseethat  Q.  =  j"  ^  =  \  log J±| , 
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And,  in  a  similar  manner,  the  values  of  Q2,  Q3,...  may 
be  calculated. 

2.  But  there  is  another  manner  of  arriving  at  these 
functions,  which  will  enable  us  to  express  them,  when  the 
variable  is  greater  than  unity,  in  a  converging  series,  with 
out  the  necessity  of  integration. 

This  we  shall  do  in  the  following  manner. 

Let  Z7=  -  — ,  v  being  not  less,  and  fi  not  greater,  than 
unity. 

Then  ^=-   _!_  £?=     J_ 


''-& 
d_ 
~dv 


am    d  (         au 


Now,  let  -  be  expanded  in  a  series  of  zonal  harmonics 

P»,  PM...PM,  so  that 

^^~  =*.(")  -P.M  +  *,M  ^,00  +  ...  +  *. 

Then  |^  {(1-^)  °^J  =  ...-i  (.'+!)  *,(••)  P,0*)+  ... 
by  the  definition  of  P  (p). 
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d  *dU  d 


And  these  two  expressions  are  equal.     Hence,  equating 
the  coefficients  of  Pf  (//.), 


Hence  ^(v)  satisfies  the  same  differential  equation  as  P£ 
and  Qt.  But  since  £7=  0  when  v  =  oo  ,  it  follows  that  fa  (v)  =  0 
when  i/=  oo  .  Hence  <£.(*>)  is  some  multiple  of  ftW^^QtW 
suppose.  It  remains  to  determine  A. 

Now,  <f>t(v)  may  be  developed  in  a  series  proceeding  by 
ascending  powers  of  —  ,  as  follows. 


We  have      -J—  =  -  +  ^  +  ...  +-^+  .... 
v  —  p      v      v  v 

and  also        =  </>»  P»  +  £»  P,  (/*)  +..-  +  *, 

Now,   by  Chap.  n.  Art.  17,  we  see  that,  if  m  be  any 
integer  greater  than  i,  the  coefficient  of  Pi  in  yu,TO  is 


m  —  i  being  always  even. 

Hence,  writing  for  m  successively  i,  i  +  2,  /  +  4,  ...  we  get 


4.  6. ../ 
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and    =^+ 


4. 6.. .(1  +  2)  J_ 

;(2t+8)(2t  +  l)...(t  +  3)  vl™ 

6.8...fr'  +  4)  1  \.      • 

+  jy-. — ~x   g  . — ON — p — ^r  T«  +••••(•«*  w  even. 

Now,  recurring  to  the  equation 


we  see  that,  if  Qt(y)  be  developed  in  a  series  of  ascending 
powers  of  -,  the  first  term  will  be  -/T(9~""W~™  '  wliere  ^ 
is  the  coefficient  of  p?  in  the  development  of  P. 


and  . 

L  ,  •*  .  0.  .  .£ 

Hence  the  first  term  in  the  development  of  Qt  (i/)  is 
2.  4.  6..  .(t-1) 

1f^ 


and     =  (f+iyF+^     if  *  b-e  even> 

which  is  the  same  as  the  first  term  of  the  development  of 
P,  (v],  divided  by  -r-.  —  =-  . 

Hence  A  =  2i+l,  and  we  have 

-~  =  <?»  W  P»  +  3  <?,  M  P,  M  +  5Q2  („)  Ps  M  +  .  .. 

3.     The  expression  for  Qt  may  be  thrown  into   a    more 
convenient  form,  by  introducing  into  the  numerator  and  de- 
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nominator  of  the  coefficient  of  each  term,  the  factor  neces 
sary  to  make  the  numerator  the  product  of  i  consecutive 
integers.  We  shall  thus  make  the  denominator  the  product 
of  i  consecutive  odd  integers,  and  may  write 


1.3.  5...(2/  +  1)  i;m     3  .  5  .  7...\2i  +  3) 


5.  7.  9..  . 


(2/j-fl)  (2&  +  3)...(2£+2&  + 
whether  i  be  odd  or  even. 

4.  We  shall  not  enter  into  a  full  discussion  of  the  pro 
perties  of  Zonal  Harmonics  of  the  Second  Kind.  They  will  be 
found  very  completely  treated  by  Heine,  in  his  Handbuch  der 
Kugelfunctionen.  We  will  however,  as  an  example,  investi 
gate  the  expression  for  -^  in  terms  of  Qi+l,  Q. 

Recurring  to  the  equation 

+  (2i  + 1 
we  see  that 


Now  we  have  seen  (Chap.  n.  Art.  22)  that 


Hence 


106        SPHERICAL  HARMONICS   OF  THE   SECOND  KIND. 


2,'+  9)  Pi+»  +  (2*  +  5) 


And  therefore  the  coefficient  of  P.(»  in  the  expansion 
f  d      1 

of  -j-        -  is 

Cl//,  J/  —  fj, 

(2i+l)  {(2 
Again, 


A     ,  J       1  d       1 

And  -j-        -  +  j-        -  =  0. 

civ  v  —  p      a/j,  v  —  /j, 

Hence,  comparing  coefficients  of  P.  (/^), 

W  =  -  (2i  +  3)  Qi+1W  -  (2»  +  7)  C«.M 


Hence  it  follows  that 


and  therefore  that 


5.  By  similar  reasoning  to  that  by  which  the»existence  of 
Tesseral  Harmonics  was  established,  we  may  prove  that  there 
is  a  system  of  functions,  which  may  be  called  Tesseral  Har 
monics  of  the  Second  Kind,  derived  from  T^  in  the  same 
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manner  as  Qt  is  derived  from  Pf.     The  general  type  of  such 
expressions  will  be 


*I?«Ml' (/*'-!)' 

and  this  when  multiplied  by  cos  a$  or  sin  a(f>,  will  give  an 
expression  satisfying  the  differential  equation 

(1  -  ft  Al*  U  +  {i  (i  +  !}(!-  ^)  -  a8]  V=  0, 


and  which  may  be   called   the  Tesseral  Harmonic  of  the 
second  kind,  of  the  degree  i  and  order  er. 


CHAPTER  VI. 

ELLIPSOIDAL  AND   SPHEROIDAL   HARMONICS. 


1.  THE  characteristic  property  of  Spherical  Harmonics 
is  thus  stated  by  Thomson  and  Tait  (p.  400,  Art.  537). 

"A  spherical  harmonic  distribution  of  density  on  a  spheri 
cal  surface  produces  a  similar  and  similarly  placed  spherical 
harmonic  distribution  of  potential  over  every  concentric 
spherical  surface  through  space,  external  and  internal." 

The  object  of  the  present  chapter  is  to  establish  the  ex 
istence  of  certain  functions  which  possess  an  analogous  pro 
perty  for  an  ellipsoid.  They  have  been  treated  of  by  Lame', 
in  his  Lemons  sur  les  fonctions  inverses  des  transcendantes  et 
les  fonctions  isothermes,  and  were  virtually  introduced  by 
Green,  in  his  memoir  On  the  Determination  of  the  Exterior 
and  Interior  Attractions  of  Ellipsoids  of  Variable  Densities, 
(Transactions  of  the  Cambridge  Philosophical  Society,  1835). 
We  shall  consider  them  both  as  functions  of  the  elliptic  co 
ordinates  (as  Lamd  has  done)  and  also  as  functions  of  the 
ordinary  rectangular  co-ordinates  ;  and  after  investigating 
some  of  their  more  important  general  properties,  shall  pro 
ceed  to  a  more  detailed  discussion  of  the  forms  which  they 
assume,  when  the  ellipsoid  is  a  surface  of  revolution. 

2.  For  this  purpose,  it  will  be  necessary  to  transform 
the  equation 


.  2T7     . 

j-s  +  -j-i  +  ~T^  =  °>  or  V  V=  °> 
dx        ay        dz 

into  its  equivalent,  when  the  elliptic  co-ordinates  e,  u,  i/  are 
taken  as  independent  variables.  If  a,  b,  c  be  the  semiaxes 
of  the  ellipsoid,  the  two  sets  of  independent  variables  are 
connected  by  the  relations 
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-±j£-  _fL=i   *2_+  £__    **  _-, 

6+6*  +  6      c"  +  e       >J  aa  +  u  +  52  +  i;  +?T^~    ' 
-^2_         3/2  g2 

+"l 


Thus  a2  +  e,  &2  +  e,  c24-  e  are  the  squares  on  the  semiaxes 
of  the  confocal  ellipsoid  passing  through  the  point  x,  y,  z. 

a?  +  v,  62  +  v  ,  c2  +  v,  the  squares  on  the  semiaxes  of  the 
confocal  hyperboloid  of  one  sheet. 

a2  +  v  ',  b2  +  v  ',  c2  +  v,  the  squares  on  the  semiaxes  of  the 
confocal  hyperboloid  of  two  sheets. 

Thus,  e  is  positive  if  the  point  #,  y,  z  be  external  to  the 
given  ellipsoid,  negative  if  it  be  internal. 

And,  if  a2  be  the  greatest,  c2  the  least,  of  the  quantities 
a',  &2,  <?, 

e  will  'lie  between  —  c2  and  oo  , 
v         »         ,>         -F    „    -c2, 


AT 

JN  ow   ^  >a  +  -^r  2  4-  ^-j  =  0  is  the  condition  that 


,dx)  TUW  Tvaf 

taken  throughout  a  certain  region  of  space,  should  be  a  mini 
mum.  In  the  memoir  by  Green,  above  referred  to,  this 
expression  is  transformed  into  its  equivalent  in  terms  of  a 
new  system  of  independent  variables,  and  the  methods  of  the 
Calculus  of  Variations  are  then  applied  to  make  the  resulting 
expression  a  minimum.  We  shall  adopt  a  direct  mode  of 
transformation,  as  follows : 

Suppose  a,  /?,  7  to  be  three  functions  of  x,  y,  z,  such  that 

such  also  that  the  three  families  of  surfaces  represented  by 
the  equations  a  =  constant,  ft  =  constant,  7  =  constant,  inter 
sect  each  other  everywhere  at  right  angles,  i.e.  such  that 
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dfidy  ,^^7  ,d/3dy^         dydz     dych     dydz_ 
dx  dx      dy  dy      dz  dz        '     dx  dx     dy  dy      dzdz~    ' 


++=0 

x  dx     dy  dy     dz  dz 
Then 


=         + 

dx       dz  dx      dfi  dx      dy  dx  ' 

<fF_  tfVfdtf?     <FZfd&\ 
dx2  ~  ~dl  (dx)  +  dP  (dx) 

[  2  d*Vd/3dy  |  ^V 

dx  dx        dydz  dx  dx        dadft  dx  dx 


do.  dx*  +  dp  da?  +  dy  drf' 


—,-?  and  -y-^  being  similarly  formed,  we  see  that,  when  the 

three  expressions  are  added  together,  the  terms  involving 
dV  dV   dV    ...  ..  ,.A.        „.       ,A 

~fT  '  ~dti'  ~J~  disappear  by  the  conditions  (1),  and  thoso 


involving  -j^-r  >  i  —  r  >   T~TO  by  the  conditions  (2).     Hence 
&  dydx  J 


4.     Now,  let 
a  = 
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All  these  expressions  satisfy  the  conditions  (1),  for  a.  is 
the  potential  of  a  homogeneous  ellipsoidal  shell,  of  proper 
density,  at  an  external  point,  and  ft  and  7  possess  the  same 
analytical  properties. 

Again,  a  is  independent  of  v  and  v,  and  is  therefore  con 
stant  when  e  is  constant.  Similarly  ft  is  constant  when  v  is 
constant,  and  7  is  constant  when  v  is  constant.  Hence  a,  ft 
7  satisfy  the  conditions  (2). 

Now 


And 


O  O  Q  \          7  r* 

ar  2/  £       1  de  _     2x 

•  \(^+W' +  (F+T2    (?+*?} £> " ^T* • 

with  similar  expressions  for  -,-  and  -v- .      Hence,  squaring 
and  adding, 


(a2  +  6)s+(6^e)2' 
But  from  the  equations 


v      62  -f  i;      c2  +  i; 


*  /3  is  a  purely  imaginary  quantity.    We  may,  if  we  please,  •write  *J—  1/3' 
forjS. 
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-__         - 

o'+i;     £*  +  u 


we  deduce 


.Vs         __^_    .   (a>~e)  (a  -  u)  (w  -  i/) 


___ 

a2  +  ft)          +  a)      c2  +  co  ~  (ftHh  a2)  (a>  "+  &2)  (aT+  c2)  ' 

w  being  any  quantity  whatever.  For  this  expression  is  of 
0  dimensions  in  CD,  e,  v,  v,  it  vanishes  when  CD  =e,  v,  or  i/, 
and  for  those  values  of  co  only,  it  becomes  infinite  when 
&  =  -  «2,  -  &",  or  —  c2,  and  for  those  values  of  co  only,  and  it  is 
=  1  when  &)  =  oo  . 

From    this,    multiplying   by   a2  +  co,    and   then    putting 
o>  =  —  a2,  we  deduce 

2  =  (e 


a  result  which  will  be  useful  hereafter. 

Again,  differentiating  with  respect  to  co,  and  then  puttin 


__ 
+  e)2      (c"  +  e/      (e  +  a<)  (e  +  6^  (e  -f- 


;«-«0(«'  -«)  («-«) ""  v)  -•-""> 


The  equation  V*F=  0  is  thus  transformed  into 


or 
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(i/  -  v)  [{(e  +  a2)  (e  +  62)  (e  +  c2)}^]  V 


+  („  _  e) 

5.  A  class  of  integrals  of  this  equation,  presenting  a  close 
analogy  to  spherical  harmonic  functions,  may  be  investigated 
in  the  following  manner.  Suppose  E  to  be  a  function  of  6, 
satisfying  the  equation 


{(€  +  a2)  (e  +  V)  (e  +  c2)}  =  (me  +  r)  E, 

m  and  r  being  any  constants. 

Then,  if  H  and  II'  be  the  forms  which  this  function 
assumes  when  v  and  v  are  respectively  substituted  for  e, 
the  equation  V2F=0  will  be  satisfied  by  V=EHH'. 

6.  We  will  first  investigate  the  form  of  the  function 
denoted  by  E,  on  the  supposition  that  E  is  a  rational  integral 
function  of  e  of  the  degree  n,  represented  by 


We  see  that 

'2, 
\en 


fe+C«{f  +  «fr+«fll*-l 

Ll<?      a)   €        )(6     C"   de\ 

+ 
=  n  F (n  -  1)  (e  +  a2)  (e  +  62)  (e  +  c2)  je""2  +  (n  -  2)  jpl€n 

+  (n"^(g"3)j?,e"-4  + 
-a)  +  (e  +  a)(6+b)  L*,i+^ _ 


•^*-ti 


F.  H. 
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Hence  writing 


(e  +  a2)  (e  +  62)  (e  +  c")  =,  e3  +  3//  +  S/2e  +/„ 
we  see  that 


(n  -  1)  (e3  +  3//  +  3/2e  +/,)    «"  +  (n  -  2)  p,e 


Hence,  equating  coefficients  of  like  powers  of  e,  we  get 
»  f  »  +  2)  =  TO- 


n  fn  —  1) 


f  3 

n  |(n  -  l)/3^fl_2  +  2/tP-i 

or,  as  they  may  be  more  simply  written, 
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•  =  nmpi  +  r, 


^{^^(-- 


It  thus  appears  that  pl  is  a  rational  function  of  r  of  the 
first  degree,  p2  of  the  second,  pn  of  the  nth,  and  when  the 
letters  pl3p2...pn  have  been  eliminated,  the  resulting  equa 
tion  for  the  determination  of  r  will  be  of  the  (n  +  l)th  degree. 
Each  of  the  letters  plt  pr..pn  will  have  one  determinate 
value  corresponding  to  each  of  these  values  of  r\  and  we 

have  seen  that  m  =  n  (n  +  - J .    There  will  therefore  be  0  + 1) 

values  of  E,  each  of  which  is  a  rational  integral  expression 
of  the  nth  degree,  n  being  any  positive  integer. 

7.  But  there  will  also  be  values  of  E,  of  the  nih  degree, 
of  the  form 


We  thus  obtain 


=  (e  +  a")i  (<•  +  V)  (6  +  c")  (n  -  1)     -2  +  (n  -  2) 


—2 
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•)  (e  +  Z,2)  (e 


+  (e  +  a2)  *  (e  +  t2)  j  (n  -  1)  je""  +  (n  -  2)  ?1e" 


)  («-l)  (ft  -2) 


Hence 


J  (e  +  62)  (e  +  c2)  +  (e  +  a2)  (e  +  c2)  +  (e  +  a2)  (e  +  Z-2)l 


+  (e  +  a2)  (e  +  &2)  (e  +  c2)  (»  - 1)  (n  -  2)   e"'3  +  («  -  3)  j 


+  (n  -  1)  (n  -  2)  [a2  + 1?  +  c2  +  (n  -  3)  ^}  =  (n  -  1)  wyt  -f  r, 


f  ~n    J  (  _l_       2   2     I        27  2  \  •     /  ON       2 

v          /  )  \  2  '  -^n~2      \          / 
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By  a  similar  process  to  that  applied  above,  we  shall  find 
that  r  is  determined  by  an  equation  of  the  nth  degree,  and 

that  m=(^-l)  f  W--    ,    and  that  each  of  the  letters       t 


<?2-  "?n-i  ig  a  rational  function  of  r.     Thus,  there  will  be  n 
solutions  of  the  form 


(6  +  b      (e  +  c2)*  {e-'  +  (n  -  1)  qlt/»  +  .  .  .  +  £_,). 

There  will  also  be  n  solutions  of  a  similar  form,  in  which 
the  factors  (e  +  c'2)*  (e  +  a2)*,  (e  +  a2)^  (e  +  62)^  are  respectively 
involved.  Hence,  the  total  number  of  solutions  of  the  n*h 
degree  will  be  4w  +  1. 

8.  We  may  now  investigate  the  number  of  solutions  of 
the  degree  71  -J-  -  ,  w  being  any  positive  integer.  These  will 

be  of  the  following  forms  :  three  obtained  by  multiplying  <i 
rational  integral  function  of  e  of  the  degree  n  by  (e  +  a2)*, 
(e  +  &*)»,  (e  +  c2)2,  respectively,  and  one  by  multiplying  a 
rational  integral  function  of  e  of  the  degree  n  —  1  by  the 
product 


An  exactly  similar  process  to  that  applied  above  will 
shew  us  that  there  will  be  n  +  1  solutions  of  each  of  the 
first  three  kinds,  and  n  of  the  fourth.  Hence  the  total  number 
of  such  solutions  -  will  be  3  (n  +  1)  +  n,  or  4w  +  3,  that  is 


To  sum  up  these  results,  we  may  say  that  the  total 
number  of  solutions  of  the  nih  degree  is  4r&  + 1,  n  denoting- 
cither  a  positive  integer,  or  a  fraction  with  an  odd  numerator, 
and  denominator  2. 

Similar  forms  being  obtained  for  H,  H' ,  we  may  proceed 
to  transform  the  expression  EHH.'  into  a  function  of  x,  yt  z. 

9.     Consider  first  the  case  in  which 

E=  f+!l9jf*+*t*=l.fj*  +  ...  +p<e 
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Write  this  under  the  form 

Then          H=(v-  ft),)  (v  -  o>8)  . . ,  (v  -  G>B), 

TJ-t          I     I  \    f     '  \  f     '  \ 

0.    =  (y-o>J  (y-ft)J  ...(V  -G)J. 

Hence 
EHH'  =(e-  ft),)  (v  -  ft),)  (v  -  ft),) . . .  (e  -  wj  (y  -  o)n)  (y'  -  WM). 

Now  we  have  shewn  (see  Art.  4  of  the  present  Chapter) 
that  (e  -  «yj  (y  -  &)J  (y'  -  ft)J 


Each  of  the  factors  of  ^ZTZ?'  being  similarly  transformed, 
we  see  that  EHH'  is  equal  to  the  continued  product  of  all 
expressions  of  the  form 


+  cw 
the  several  values  of  co  being  the  roots  of  the  equation 

o. 


As  "this  equation  has  been  already  shewn  to  have  (n  +  1) 
distinct  forms,  we  obtain  (n  +  1)  distinct  solutions  of  the 
equation  V2F=0,  each  solution  being  the  product  of  n 
expressions  of  the  form 


That  is,  there  "will  be  n  +  1  independent  solutions  of  the 
degree  2n  in  x,  y,  z,  each  involving  only  even  powers  of  the 
variables. 

10.  To  complete  the  investigation  of  the  number  of  solu 
tions  of  the  degree  2n,  let  us  next  consider  the  case  in  which  E 
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The  object  here  will  be  to  transform  the  product 

(e  +  V)*  (v  +  V}*  (V  +  Z>2)*  (c  +  c2)'  (i;  +  c2)*  (v  +  c2)*, 

since  the  other  factors  will,  as  already  shewn,  give  rise  to  the 
product  of  7i  —  1  expressions  of  the  form 


a2  +  ft)     b2  +  ft)     c2  +  &) 

Now,  by  comparison  of  the  value  of  x2  given  in  Art.  4, 
we  see  that 

(e  +  62)  (v  +  b2)  (v  +  b2)  (e  +c2)  (i;  +  c2)  (i/  +  c2) 

=  (b2  -  c2)  (b2  -  a2)  (c2  -  a2)  (c2  -  b2)  fz\ 

Hence,  we  obtain  a  system  of  solutions  of  the  form  of 
the  product  of  (n  —  1)  expressions  of  the  form 

x*  2/2  z2 

multiplied  by  yz.  Of  these  there  will  be  n,  and  an  equal 
number  of  solutions  in  which  zx,  xy}  respectively,  take  the 
place  of  yz. 

Thus,  there  will  be  4tn  + 1  solutions  of  the  degree  2??,  in 
the  variables  of  which  n  +  I  are  each  the  product  of  n 
expressions  of  the  form 

a2  +  w      b2  -F&>     c2  +  G)  ~    J 

n  are  each  the  product  of  (n—l)  such  expressions,  multiplied 

by  yz, 
n  ...  ...  ...  zx, 


11.  We  may  next  proceed  to  consider  the  solutions  of  the 
degree  2n  + 1  in  the  variables  x,  y,  z. 

Consider  first  the  case  in  which 
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Here  the  product  (e  +  a2)*  (v  +  a2)*  (v  +  a2)*  will,  as  just 
shewn,  give  rise  to  a  factor  x  in  the  product  EHH'. 

Hence  we  obtain  a  system  of  solutions  each  of  which  is 
the  product  of  n  expressions  of  the  form 

**      ,      f      .      *2         ! 
2  * 


multiplied  by  x.  Of  these  there  will  be  n  +  1,  and  an  equal 
number  of  solutions  in  which  y,  z,  respectively  take  the 
place  of  the  factor  x. 

Lastly,  in  the  case  in  which 

rf  (e  +  irf  (e  +  c2 


we  see  that  in  EIIH'  the  product 

^  (v+a2)*  (i/+a2)*  (e-f  b^  (v+tf)1*  (v'  +  tf)* 


will  give  rise  to  a  factor  xyz. 

Hence  we  obtain  a  system  of  solutions  each  of  which  is 
the  product  of  (n  —  1)  expressions  of  the  form 


n. 


multiplied  by  xyz.     Of  these  there  will  be 

Thus  there  will  be  4<n  +  3  solutions  of  the  degree  2n  +  l 
in  the  variables,  of  which 

(n  +  1)  are  each  the  product  of  n  expressions  of  the  form 

y?  if  z* 

-r-  -  +  TT^  --  h  -3—  --  1  multiplied  by  x, 
a  +  a)      b  +  co      c2  -f  co 

(n  +  1)  are  each  the  product  of  n  such  expressions,  multiplied 
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n  are  each  the  product  of  (n  -  1)  such  expressions,  multi 
plied  by  xyz. 

12.     Now  an  expression  of  the  form  G  .  EH  IT,  C  being 
any  arbitrary  constant,  is  an  admissible  value  of  the  potential 

at  any  point  within  the  shell   ^  +  £2+c*  =  1'      ^ut   ^   *s 

not  admissible  for  the  space  without  the  shell,  since  ^  it 
becomes  infinite  at  an  infinite  distance.  The  factor  which 
becomes  infinite  is  clearly  E,  and  we  have  therefore  to 
enquire  whether  any  form,  free  from  this  objection,  can  be 
found  for  this  factor.  We  shall  find  that  forms  exist,  bearing 
the  same  relation  to  E  that  zonal  harmonics  of  the  second 
kind  bear  to  those  of  the  first. 

Now  considering  the  equation 

=me+r)  U, 


which  we  suppose  to  be  satisfied  by  putting  U=  E,  we  see 
that,  since  it  is  of  the  second  order,  it  must  admit  of  another 

particular  integral.     To  find  this,  substitute  for   U,  E  Ivde, 
we  then  have 


=  [{(6+  a2)  (e  +  J2)  (e  +  c2)}*  I]  E.jvde 


=  [{(*  +  a2)  (e  +  V)  (e  +  c2))*  |J  E  .  jvde 
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Now,  since  by  supposition,  the  equation  for  the  determi 
nation  of  U  is  satisfied  by  putting   U=E,  it  follows  that 

when  E  Ivde  is  substituted  for  U,  the  terms  involving  \vcle 

will  cancel  each  other,  and  the  equation  for  the  determina 
tion  of  v  will  be  reduced  to 


ldv,2_dE     1 
vde+E'de  +2 

whence    log  v  +  2  log  E  +  log  {(e  +  a2)  (e  +  62)  (e  +  c2))* 

=  log  u0+  2  log  .#0+  log  a&c, 
v0  and  .Z?0  being  the  values  of  v  and  ^,  corresponding  to  e  =  0. 

Hen  .*' ^ 


We  may  therefore  take,  as  a  value  of  the  potential  at 
any  external  point, 

de 


'0       •  TT>     f/  o. 

/  e  ^2  {(e  +  a  ) 

For  this  obviously  vanishes  when  e  =  op .  It  remains  so 
to  determine  v0  that  this  value  shall,  at  the  surface  of  the 
ellipsoid,  be  equal  to  the  value  C.EHH',  already  assumed 
for  an  internal  point.  This  gives 


C-  v  Efabc 


Hence,  putting  v0 .  E* .  abc  =  F0,  we  see  that  to  the  value 
of  the  potential 

V.EHH'T d€ 

Jo  E*  l\ 
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for  any  internal  point,  corresponds  the  value 

cle 


V.EHH' 


for  any  external  point. 

13.  We  proceed  to  investigate  the  law  of  distribution  of 
density  of  attracting  matter  over  the  surface  of  the  ellipsoid, 
corresponding  to  such  a  distribution  of  potential. 

Now,  generally,  if  §n  be  the  thickness  of  a  shell,  p  its 
volume  density,  the  difference  between  the  normal  compo 
nents  of  the  attraction  of  the  shell  on  two  particles,  situated 
close  to  the  shell,  on  the  same  normal,  one  within  and  the 
other  without  will  be  4<7rp$n.  This  is  the  attraction  of  the 
shell  on  the  outer  particle,  minus  the  attraction  on  the  inner 
particle. 

But  the  normal  component  of  the  attraction  on  the  outer 

,   .      dV 
particle  estimated  inwards  is  —  -r-  . 

And,  if  V  denote  the  potential  of  the  shell  on  an  in 
ternal  particle,  the  normal  component  of  the  attraction  on 

-   .       dV 

it  estimated  inwards  is  —  ,  -  . 

an 

TT  *  dV'        dV 

Hence  4^  =  _  -  _  . 


dV_dV.dx     dVJy     dVdz_ 
dn  ~~  dx  dn      dy  dn      dz  dn  ' 

And  -7-  is  the  cosine  of  the  inclination  of  the  normal  at 
dn 

the  point  #,  y,  z  to  the  axis  of  x,  and  is  therefore  generally 

/v» 

equal  to  e  —^  —  -  ,  e  denoting   the   perpendicular   from   the 

Cb     ~T~  € 

centre  on  the  tangent  plane  to  the  surface 
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And  we  have  shewn  that 


flf  = 


(a"  -  62)  (aa  -  c2) 


whence  -  —  =  -  — 

x  de      a?  +  e  ' 


or 


•  •  -r-  =  e  -y-  . 

a/i  de 

Similarly            ^=2e^  ^- 

'  " 


de 


dx  +  dVdy     dVdz^v 

de       dy  de      dz  de)  ~    *  de  ' 

Similarly  . =  2e  — 

dn  de  ' 

Now     V  =  F. .  7777^'  r     * 


And       V'=V0.EHH' 

therefore,  generally, 

^F      lr          „  dE  r  de 


•  ^•{(t+'rfJCf+iPXf^^ 

efe 


-  V^.EHH' 


But,  when  the  attracted  particle  is  in   the   immediate 
neighbourhood  of  the  surface,   e  =  0.     Hence,  the  first  line 


ELLIPSOIDAL  AND   SPHEROIDAL   HARMONICS.  125 

dV  dV 

of  the  value  of  —  becomes  identical  with  the  value  of  —  s 
de  de 

and  we  have 


de        de  "      °   E    abc  ' 


E0  denoting  the  value  which  E  assumes,  when  e  =  0. 

' 


Hence,  ivpSn  =  2e  F0  -- 


But  Sn,  being  the  thickness  of  the  shell,  is  proportional  to 

/?          r/ 

e,  and  we  may  therefore  write  ~-  =  K-  ,  £a  being  the  thick 
ness  of  the  shell  at  the  extremity  of  the  greatest  axis  ; 


*         ZirSaabc    E0    ' 

and  this  is  proportional  to  the  value  of  V  corresponding  to 
any  specified  value  of  e,  since  HH'  is  the  only  variable 
factor  in  either. 

Hence  functions  of  the  kind  which  we  are  now  considering 
possess  a  property  analogous  to  that  of  Spherical  Harmonics 
quoted  at  the  beginning  of  this  Chapter.  On  account  of 
this  property,  we  propose  to  call  them  Ellipsoidal  Harmonics, 
and  shall  distinguish  them,  when  necessary,  into  surface  and 
solid  harmonics,  in  the  same  manner  as  spherical  harmonics 
are  distinguished.  They  are  commonly  known  as  Lame's 
Functions,  having  been  fully  discussed  by  him  in  his  Lemons. 
The  equivalent  expressions  in  terms  of  x,  y,  z  have  been  con 
sidered  by  Green  in  his  Memoir  mentioned  at  the  beginning 
of  this  chapter,  and  for  this  reason  Professor  Cayley  in  his 
"  Memoir  on  Prepotentials,"  read  before  the  Royal  Society 
on  June  10,  1875,  calls  them  "  Greenians." 

We  may  observe  that  the  factor 

1     a     1 
4-7T  Ba  abc 

is  equal  to  -r  V-K-  ,  and  therefore  also  to  7   —  ^r  or  -r 
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Hence,  it  is  equal  to 


-Q-  (bc$a  +  ca&b 
o 

1 
or  to 


volume  of  shell ' 
and  the  potential  at  any  internal  point 

=  i  volume  of  shell  x  EE0 .  p  [°°  — - — 

and  the  potential  at  any  external  point 

=  i  volume  of  shell  x  EEn .  p  (   ^L 


where  for  p  must  be  substituted  its  value  in  terms  of  v  and  v. 

^  14.    We  will  next  prove  that  if  Vlt  F2  be  two  different 
ellipsoidal  harmonics,  dS  an  element  of  the  surface  of  the 

ellipsoid,   lie  F1F2c?>S  =  0,  the  integration  being  extended  all 
over  the  surface. 

We  have  generally 

jj[(  FXV2  F2  -  F2V2  FJ  dxdydz  =  jf(  Vl  d^  -  F2  ^ 


And  throughout  the  space  comprised  within  the  limits  of 
integration,  V2  Vl  =  0,  V2  F2  =  0.     Hence 


Now  it  has  been  shewn  already  that  Vlt  F2  are  each  of 
the  form  EHH',  where  E  is  a  function  of  e  only,  H  the  same 
function  of  v,  Hr  of  v.  We  may  therefore  write 


and  similarly  F,  =/a  (e)/,  («)/„  (u). 
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Hence  v*L-wfLQ 

KI  de  -  *  v;c<o  ' 
y4£-rr£®- 

2  de  ~     l    2e   ' 


*de          *  de   ~      l 
Now,  all  over  the  surface,  e  =  0.     Hence 


Hence,  unless  -  2.        —  '\..1  .      =  0,  which  cannot  happen 
/>  W)     /i  v^J 

unless  the  functions  denoted  by  f^  and  ft  are  identical*,  or 
only  differ  by  a  numerical  factor,  we  must  have 


Now  e  is  proportional  to  the  thickness  of  the  shell  at 
any  point.     Calling  this  thickness  80,  we  have  therefore 


Hence,  adding  together  the  results  obtained  by  integrating 
successively  over  a  continuous  series  of  such  surfaces,  we  get 


Fj ,  F2  now  representing  solid  ellipsoidal  harmonics,  and  the 
integration  extending  throughout  the  whole  space  comprised 
within  the  ellipsoid. 

*  This  may  be  shewn  more  rigorously  by  integrating  through  the 
space  bounded  by  two  confocal  ellipsoids,  denned  by  the  values  \  and  /x  of  e. 
We  then  get,  as  in  the  text, 


Now  the  factor  within  { }  cannot  vanish  for  all  values  of  X  and  p,  unless  the 
functions  devoted  by  /j,  and  /2  be  identical,  or  only  differ  by  a  numerical 
factor. 
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15.     It  will  be  well  to  transform  the  expression 

jj'WdS         *  '     . 

to  its  equivalent,  in  terms  of  v,  v. 

For  this  purpose  we  observe  that  if  ds,  ds  be  elements  of 
the  two  lines  of  curvature  through  any  point  of  the  ellipsoid, 
dS  =  ds  ds. 

Now, 

c?s2  is  the  value  of  dx*  +  cZ?/2  +  dz*  when  e  and  v  are  constant, 
ds'2  ...  ...  ...  e  and  v 


2  _  (e  +  ag)  (i;  +  a2)  (v  +  a2) 
*T*a 


therefore  if  e  and  v  '  do  not  vary, 

dv 


-,^7 

.'.  a^c  =  -  —  -  au. 

2  a  -1-  v 

Similarly 


- 


(6'+  uj'      (c*  •+  v)  ) 

Again,  differentiating  with  respect  to  co  the  expression 

x2  tf  zz 

obtained  for  -=  --  \-  -^  —  -f  -  ---  1,  we  get 
a  4-  a)      b~  +  co     c  +  co 

y?  f  z*  (v-a)}  (u'-ft)) 

^+~)2  +  -^+^Y  +  (<r  +  a))2  ~  ^H-  o>)^"+  "wKcN-  o) 


i>'  —  G>)  (e  —  o))  (e  —  6?)  (v  —  ft))  (i/  —  w) 

'  " 
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therefore,  putting  &>  =  v, 

(a2  4-  v)'2      (62  -f  v)*      (ca  +  vf     (a2  +  v)  (52  +  v)  (c*  +  v)  ' 
^2_1          (i/-v)(e-v)          ,, 


A  similar  expression  holding  for  ds'*  we  get 

2==      1 (v'-v)*(e-v)(e-v')   

"  16  (a*+v)  (6*  +  i;)  (c2+v)  («2+'      '  2' 

1  ^2  v2 

Again,  -y  =   7-2-— vi  +  /72     , 


(c2+e) 


writing  e  for  o>  in  the  expression  above  ; 


16  (aa+L>)  (6J+f)  (c2+v)  (a2+u)  (62+fr)  (c2+n 

It  has  been  shewn  that,  integrating  all  over  the  surface, 
the  limits  of  v  are  -  c2  and  —bz,  those  of  i/,  —  b2  and  —  a2. 

Hence,    FA,   F2,   denoting  two  different  ellipsoidal  har 
monics 


The  value  of  the  expression  1  1  1  V*dxdydz,  or  its  equiva 
lent 

"'8  f'6'  V*(v-v}dvdv' 

-b-2  J  -a?    a'+  v    ^ 


ale 
}C 


in   any  particular   case,  is   most  conveniently  obtained  by 
expressing  Fas  a  function  of  x,  y,  z. 

F.  H.  9 
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16.  Before  proceeding  further  with  the  discussion  of  ellip 
soidal  harmonics  in  general,  we  will  consider  the  special  case 
in  which  the  ellipsoid  is  one  of  revolution.  We  must  enquire 
what  modification  this  will  introduce  in  the  quantities  which 
we  have  denoted  by  a,  ft,  7,  viz. 


a=f°- 

J  .    a2 


d 


+  *  +          c  + 


and  in  the  differential  equation 


__          ._         ^v_=st 

We  will  first  suppose  the  axis  of  revolution  to  be  the 
greatest  axis  of  the  ellipsoid,  which  is  equivalent  to  supposing 
tf  =  c2.  To  transform  a  and  7,  put  a2  +  ^  =  02,  a2  +  e  =  17*, 
a2  +  i/  =  ft)2  ;  we  then  obtain 


= 
" 


= 


To  transform  /3,  we  must  proceed  as  follows. 

Put    i/r  =  _  c2  cos2  -or  -  62  sin2  w,    u  =  -  c2  cos2  </>  -  62  sin2  $, 
we  then  get  generally 

62  +  ^  =  (62  -  c2)  cosV,    c2  +  ^  =  (c2  -  6')  siriV  ; 

C?-x|r  =  2  (c2  —  62)  COS  «r  sin  -or  C^tzr  ; 
.       Q_ 


V  ^1   *  (a8  -  62 
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Hence,  ** 


.     2  2 

-j-  =  ~  (w  -  cr  4-  62)  -r-, 
ay     2  '  6/a> 


Also,  e  =  ?;2  -  a2,  u  =  a)'2  -  a2,  u  =  -  b\  and  our  differential 
equation  becomes 


-  a2  4-  J2)     a)2  -  a2  +  62     -      7 


or  a>2  -  a2  +  J2          -az 


This  equation  may  be  satisfied  in  the  following  ways. 

First,  in  a  manner  altogether  independent  of  0,  by  sup 
posing  F  to  be  the  product  of  a  function  of  77  and  the  same 
function  of  o>,  this  function,  which  we  will  for  the  present 
denote  by/  (17)  or/(o>),  being  determined  by  the  equation 
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CL   V 

Secondly,  by  supposing  ->  ,2   a  constant  multiple  of    F, 
—  o-2  F,  suppose. 

Our  equation  may  then  be  written 


-  <r2  (a2  -  62)  {(a>2  -  a2  +  62)  -  (^  -  a2  +  62))  F=  0, 

which  may  be  satisfied  by  supposing  the  factor  of  V  inde 
pendent  of  <f>  to  be  of  the  form  F  (rj)  F  (o>),  where 


(f-  <f  +  V]  ~F(rf>  -  °*  0"  -  62 
I  a^j 


The  factor  involving  (/>  will  be  of  the  form 

A  cos  a-(f)+£  sin  cr<£. 
Now,  returning  to  the  equation 


we  see  that,  supposing  the  index  of  the  highest  power  of  77 
involved  mffy)  to  be  it  we  must  have  m  =  i  (i  +  1). 


Now,  it  will  be  observed  that  rj  may  have  any  value 
however  great,  but  that  o>2,  which  is  equal  to  a2  +  v,  must 
lie  between  a2  —  62  and  0.  Hence,  putting  o>2  =  (a2  —  62)  yu-2, 
where  //  must  lie  between  0  and  1,  we  get 

=0. 
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Hence  this  equation  is  satisfied  by  f{(a*  —  &9)^}  =  CPt, 
C  being  a  constant ;  and  supposing  0=1  we  obtain  the 
following  series  of  values  for  /  (o>), 

t  =  0,    /(o))  =  l, 

••=i.  /M-rs-2*!. 


o      j?/   \      o®  —  ^cr  — q  J 

l  =  2'  /(m)=     2(J-6-)     ' 

t  =  3,  /(»)  =  ' 


2  (a* - 


Exactly  similar  expressions  may  be  obtained  for  /(??),  and 
these,  when  the  attraction  of  ellipsoids  is  considered,  will 
apply  to  all  points  within  the  ellipsoid.  But  they  will  be 
inadmissible  for  external  points,  since  rj  is  susceptible  of  in 
definite  increase. 

The  form  of  integral  to  be  adopted  in  this  case  will  be 
obtained  by  taking  the  other  solution  of  the  differential 
equation  for  the  determination  of  f(if),  i.e,  the  zonal  har 

monic  of  the  second  kind,  which  is  of  the  form  Q. 
where 

de 


=  P. 

-  J 


Or,  putting  r/2  =  (a2  -  £'2)  i/2,  ^  =  (a2  —  62)  A,2,  we  may  write 

rT\ 
p,x,±      ir 
•s  Wj  (A  —  -1) 

17.  We  may  now  consider  what  is  the  meaning  of  the 
quantities  denoted  by  TJ  and  a).  They  are  the  values  of  Sr 
which  satisfy  the  equation 
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and  are  therefore  the  semi-axes  of  revolution  of  the  surfaces 
confocal  with  the  given  ellipsoid,  which  pass  through  the 
point  x,  y,  z.  One  of  these  surfaces  is  an  ellipsoid,  and 
its  semi-axis  is  TJ.  The  other  is  an  hyperboloid  of  two  sheets 
whose  semi-axis  is  o>. 

Now,  if  6  be  the  eccentric  angle  of  the  point  x,  y,  z, 
measured  from  the  axis  of  revolution,  we  shall  have 

a?  =  7f  cos2  6. 

But  also,  since  rf,   co2,  are  the  two  values  of  ^2  which 
satisfy  the  equation  of  the  surface, 

i;V=±  (a*  -&*)*> 

Hence  a>2  =  (a2  -  62)  cos2  0, 

and  we  have  already  put 


whence  the  quantity  which  we  have  already  denoted  by  //, 
is  found  to  be  the  cosine  of  the  eccentric  angle  of  the  point 
x,  y,  z  considered  with  reference  to  the  ellipsoid  confocal 
with  the  given  one,  passing  through  the  point  xt  y,  z.  We 
have  thus  a  method  of  completely  representing  the  potential 
of  an  ellipsoid  of  revolution  for  any  distribution  of  density 
symmetrical  about  its  axis,  by  means  of  the  axis  of  revo 
lution  of  the  confocal  ellipsoid  passing  through  the  point 
at  which  the  potential  is  required,  and  the  eccentric  angle 
of  the  point  with  reference  to  the  confocal  'ellipsoid.  For 
any  such  distribution  can  be  expressed,  precisely  as  in  the 
case  of  a  sphere,  by  a  series  of  zonal  harmonic  functions  of 
the  eccentric  angle. 

18.  When  the  distribution  is  not  symmetrical,  we  must 
have  recourse  to  the  form  of  solution  which  involves  the  factor 
A  cos  o-<f)  +  B  sin  ox£.  It  will  be  seen  that,  supposing  F  to 
represent  a  function  of  the  degree  it  and  putting  m  =  i  (t+1), 
the  equation  which  determines  F(o>)  is  of  exactly  the  same 
form  as  that  for  a  tesseral  spherical  harmonic.  For  F(ij)t  if 
the  point  be  within  the  ellipsoid,  we  adopt  the  same  form, 
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if  without,  representing  the  tesseral  spherical  harmonic  by 

_  !?  —  .1  1  or  T&  (v),  we  adopt  the  form 

a  -b) 


19.  It  may  be  interesting  to  trace  the  connexion  of  sphe 
rical  harmonics  with  the  functions  just  considered.  This  may 
be  effected  by  putting  tf  =  a2.  We  see  then  that  77  will  become 
equal  to  the  radius  of  the  concentric  sphere  passing  through 
the  point,  and  jf  -  a?  4-  62  will  become  equal  to  rf.  Hence 
the  equation  for  the  determination  of/  (77)  will  become 


which  is  satisfied  by  putting  f(vj)  =  V,  or  T?' (m).  The  former 
solution  is  adapted  to  the  case  of  an  internal,  the  latter  to 
that  of  an  external  point. 

With  regard  to  /(&>),  it  will  be  seen  that  the  confocal 
hyperboloid  becomes  a  cone,  and  therefore  co  becomes  inde 
finitely  small.  But  /*,  which  is  equal  to  2^  i,  remains 

finite,  being  in  fact  equal  to  -  or  cos  0.  Hence /(/A)  becomes 
the  zonal  spherical  harmonic. 

Again,  the  tesseral  equations,  for  the  determination  of 
F(rj),  F(ai)t  become 


which  are  satisfied  by  F(rf)  =rf  or  ^  {l+l). 

And,  writing  for  o>2,  (a2- 62)  ju,2,  we  have,  putting  F(w)  = 

rhich 
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20.  ^  We  will  next  consider  the  case  in  which  the  axis  of 
revolution  is  the  least  axis  of  the  ellipsoid,  which  is  equi 
valent  to  supposing  a2=62.  To  transform  a  and  /3,  put 
c2  +  T/T  =  #2,  c2  +  e  =  rfy  <?  +  v  =  o>2,  we  thus  obtain 


d9 


CO 


To  transform  7,  we  must  proceed  as  follows : 

Put   yjr  =  —  a2  sin2  ^r  —  &2  cos2  •cr,     v  =  —  a?  sin2  <p  —  l*  cos 
we  then  get,  generally, 

a2  +  ^  =  (a2  —  i2)  cos2  OT,     &2  +  'vjr  =  —  (a?  —  62)  sin2  CT, 
c  -f  i/r=c2— a>2sin2<£— ^2cos20,  67^  =  — 2  (a2— 62)  sin'zxrcos'GT  c? 
Hence 
=  2r° ^OT 2j>       .     2_ 

*  /  /        9          •         O  .70  «  ->v  1  "~"  TTT     II      ft       ==     0     . 


Hence, 


also, 


u  =-a  , 
and  our  differential  equation  becomes 


'« 

-(a2-c2  +  ^)^a2-c2  +  to2)4iV 

aw 
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We  will  first  consider  how  this  equation  may  be  satisfied 
by  values  of  V  independent  of  </>. 

We  may  then  suppose  Fto  be  the  product  of  a  function 
of  TI,  and  the  same  function  of  &>,  this  function,  which  we  will 
suppose  to  be  of  the  degree  it  being  determined  by  the 
equation 


ft) . 


On  comparing  this  with  the  ordinary  differential  equa 
tion  for  a  zonal  harmonic,  it  will  be  seen  that,  on  account 
of  a2  being  greater  than  c2,  the  signs  of  the  several  terms  in 
the  series  for  f(tj)  will  be  all  the  same,  instead  of  being 
alternately  positive  and  negative.  We  shall  thus  have 


and  generally 


We  will  denote  the  general  value  of/  (i;)  by  p4  \  —  —  —i 

\  \G/  —  C  J  " 

or,  writing  TJ  =  (a2  -  C2)^i>,  by  pt(v). 
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For  external  points,  we  must  adopt  for  /  (77)  a  function 

which  we  will  represent  by  qt\ ^— if ,  or  g. (v),  which  will 

l(a2  — C2)*J 
be  equal  to 

ae 


It  is  clear  that  /(o>)  may  be  expressed  in  exactly  the 
same  way.  But  it  will  be  remembered  that  rf  and  co*  are 
the  two  values  of  S-2  which  satisfy  the  equation 


Hence  77,  as  before,  is  the  semi-axis  of  revolution  of  the 
confocal  ellipsoid  passing  through  the  point  (x,  y,  z).  But 
vfu**  =  —  (a2  —  c2)  32,  an  essentially  negative  quantity,  since 
a2  is  greater  than  c2.  Hence  <w2  is  essentially  negative.  Now, 
if  0  be  the  eccentric  angle  of  the  point  (x,  y,  z)  measured 
from  the  axis  of  revolution,  we  have  z2  =  rf  cos2#.  Hence 


and  therefore  w2  =  —  (a2  —  c2)  cos2  9 

=  —  (a2  —  c2)  ff,  suppose. 

Hence  the  equation  for  the  determination  of  f(a))  assumes 
the  form 


the  ordinary  equation  for  a  zonal  spherical  harmonic.    Hence 
we  may  write 


fjb  being  the  cosine  of  the  eccentric  angle  of  the  point  x,  y,  z, 
considered  with  reference  to  the  confocal  ellipsoid  passing 
through  it. 
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21.  We  have  thus  discussed  the  form  of  the  potential, 
corresponding  to  a  distribution  of  attracting  matter,  sym 
metrical  about 'the  axis.  When  the  distribution  is  not 
symmetrical,  but  involves  <f>  in  the  form  A  cos  crc/>  -f  B  sin  cr(£, 
we  replace,  as  before,  Pt  (/&)  by  T^  (//),  and  p.(p)  by  a 
function  t^(v)  determined  by  the  equation 


22.     As  an  application  of  these  formulae,  consider  the  fol 
lowing  question. 

Attracting  matter  is  distributed  over  the  shell  whose 
surface  is  represented  by  the   equation  —^  +      7.2  —  =  1>    so 

that  its  volume  density  at  any  point  is  Pt  (//-),  //,  being  the 
cosine  of  the  eccentric  angle,  measured  from  the  axis  of 
revolution  ;  required  to  determine  the  potential  at  any 
point,  external  or  internal. 

The  potential  at  any  internal  point  will  be  of  the  form 

CP.MP.V)  ........................  (i), 

and  at  an  external  point,  of  the  form 

QM  ........................  (2), 


where   (a2  —  62)^  v  =  the  semi-axis  of  the  figure  of  the  con- 
focal  ellipsoid  of  revolution  passing  through  the  point  (//,,  v). 

Now  the  expressions  (1)  and  (2)  must  be  equal  at  the 

surface  of  the  ellipsoid,  where  v  =  —  -  -  —  ,  . 

(a  —  b  )'* 

Hence 
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But  generally 

ft  (v)  =  Pt  (i/)  f dX . 

'-/•^WrCX"-!) 
Hence 

Q\      a      l  =  p[      «       ]f 
"  V  -  »5»)       "'  V  -  &>$  J  _"_  7 


We  may  therefore,  putting  (7'  =  ^7J.  J  -  -  —  ri  ,  write 

I  (a2  -  6s)  *J 


and  we  thus  express  the  potentials  as  follows  : 

i  W  Pt  W  C,  ].  gg|al[  at  an  internal  point,  , 
l(o  —  6§;*J 

r^ir    at  an  external  p°int- 


Or,  substituting  for  Q.  its  value  in  terms  of  Pit 

Vl  =  AP,  (fjL)  Pt  (v)  Pi  I  -  a- 
u    ' 


at  an  internal  point, 
F  = 


at  an  external  point. 

Now,  to  determine  A,  we  have,  Sa  being  the  thickness 
of  the  shell  at  the  extremity  of  the  axis  of  revolution, 
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.!_  A_  (*ZL-*L\ 

.  ij  \  drj         dr)  )1]  =  a 


a 


a*  -b*\dv       dv 


—5-  -1 

a2  -  62 


Hence,  if  p  =  Pt  (ft),  we  obtain 


a 
And  we  thus  obtain 


4^36  P.Oi)  P,  W  ft 


If  the  shell  be  represented  by  the  equation 
a?+y*      z 
-aT*?-1' 
it  may  be  shewn  in  a  similar  manner  that  we  shall  have 
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Vl  =  4aradaPi(p)pi 

V  = 


23.  We  may  apply  this  result  to  the  discussion  of  the 
following  problem. 

//  the  potential  of  a  shell  in  the  form  of  an  ellipsoid  of 
revolution  about  the  greatest  point  be  inversely  proportional 
to  the  distance  from  one  focus,  find  the  potential  at  any 
internal  point,  and  the  density. 

If  the  potential  at  P  be  inversely  proportional  to  the 
distance  from  one  focus  $,  and  //  be  the  other  focus,  we  have, 

HP+SP  =  2r),     HP-SP=2a>, 

.'.     SP=1}-<0. 

Hence  if  M  be  the  mass  of  the  shell,  F"2  the  potential  at 
any  external  point, 


77  —  a) 
M 


Now,  by  what  has  just  been  seen,  the  internal  potential, 
corresponding  to  Pi  (/LI)  Qi  (v),  is 

Q, 


P 


Hence,  if  F2  be  the  potential  at  any  internal  point, 
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And  the  volume  density  corresponding  to  Pt  (^}  Qi  (v)  is 


Hence  the  density  corresponding  to  the  present  distri 
bution  is 


p 

• 


If  F2  had  varied  inversely  as  JGTP,  we  should  have  had 


and  our  results  would  have  been  obtained  from  the  foregoing 
by  changing  the  sign  of  o>,  and  therefore  of  p. 

24.  Now,  by  adding  these  results  together,  we  obtain 
the  distributions  of  density,  and  internal  potential,  corre 
sponding  to 


7?-W         77 

or,  in  geometrical  language, 


SP.HP  ' 

=  M  multiplied  by  the  axis  of  revolution  of  the  confocal 
ellipsoid,  and  divided  by  the  square  on  the  conjugate  semi- 
diameter.  We  may  express  this  by  saying  that  the  potential 
at  any  point  on  the  ellipsoid  is  inversely  proportional  to  the 
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square  on  the  conjugate  semi-diameter,  or  directly  as  the 
square  on  the  perpendicular  on  the  tangent  plane. 

Corresponding  to  this,  we  shall  have,   writing  2k  for  i, 
since  only  even  values  of  i  will  be  retained, 


A;  being  0,  or  any  positive  integer. 

Again,  subtracting  these  results  we  get 


7]  —  O) 


=  M  multiplied  by  the  distance  from  the  equatoreal  plane, 
and  divided  by  the  square  on  the  conjugate  semi-  diameter. 

This  gives,  writing  2k  +  1  for  i, 

0      f      a      I 

lla'-Wm 


25.  In  attempting  to  discuss  the  problem  analogous  to 
this  for  an  ellipsoid  of  revolution  about  its  least  axis,  we  see 
that  since  its  foci  are  imaginary,  the  first  problem  would  re 
present  no  real  distribution.  But  if  we  suppose  the  external 
potential  to  be  the  sum  or  difference  of  two  expressions,  each 
inversely  proportional  to  the  distance  from  one  focus,  we 
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obtain  a  real  distribution  of  potential — in  the  first  case 
inversely  proportional,  to  the  square  on  the  conjugate 
semi-diameter,  in  the  latter  varying  as  the  quotient  of  the 
distance  from  the  equatoreal  plane  by  the  square  on  the 
conjugate  semi- diameter. 

It  will  be  found,  by  a  process  exactly  similar  to  that  just 
adopted,  that  the  distributions  of  internal  potential,  and 
density,  respectively  corresponding  to  these  will  be ; 

In  the  first  case 


00 


k  being  0,  or  any  positive  integer. 
In  the  second  case 


7c  being  0,  or  any  positive  integer. 

26.  We  may  now  resume  the  consideration  of  the  ellip 
soid  with  three  unequal  axes,  and  may  shew  how,  when  the 
potential  at  every  point  of  the  surface  of  an  ellipsoidal  shell 
is  known,  the  functions  which  we  are  considering  may  be 
employed  to  determine  its  value  at  any  internal  or  external 
point.  We  will  begin  by  considering  some  special  cases, 

F.  H.  10 
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by  which  the  general  principles  of  the  method  may  be  made 
more  intelligible. 

27.     First,  suppose  that  the  potential  at  every  point  of 
the  surface  of  the  ellipsoid  is  proportional  to  x  =  —  —  suppose. 

& 

In  this  case,  since  x  when  substituted  for  V,  satisfies  the 
equation  v27=  0,  we  see  that  F0-  will  also  be  the  potential 

at  any  internal  point.  But  this  value  will  not  be  admissible 
at  external  points,  since  x  becomes  infinite  at  an  infinite 
distance. 

Now,  transforming  to  elliptic  co-ordinates 


And  the  expression 

ro{(*+**)(v+<f)(v'+a*)}*  f  __ 
a  \     (a2-&2)(a2-c2)     j    J.  •      +a' 


*£ 


(^  +  a2)  (0/r  +  a2)  (i/r  +  6 
satisfies,  as  has  already  been  seen,  the  equation  y27=0,  is 
equal  to  F0-  at  the  surface  of  the  ellipsoid,  and  vanishes 

at  an  infinite  distance.     This  is  therefore  -the  value  of  the 
potential  at  any  external  point.    It  may  of  course  be  written 

tty 


a 


28.     Next,  suppose  that  the  potential  at  every  point  of 
the  surface  is  proportional  to  #s=F0f-,  suppose.     In  this 

OG 
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case,  as  in  the  lasfc,  we  see  that,  since  yz  when  substituted 
for  F,  satisfies  the  equation  v2^=  0,  tne  potential  at  any 

internal  point  will  be  F0  y— ;  while,  substituting  for  y,  z  their 

values  in  terms  of  elliptic  co-ordinates  we  obtain  for  the 
potential  at  any  external  point 


be  J*  (^  +  tf)  (\fr  +  c2)  [ty  +  a2)  (fy+b 


T 

J  a 


C2)  {(^r  +  a2)  (Vr  +  62).  (t  +  C2)]*  " 

29.     We  will  next  consider  the  case  in  which  the  po- 

2 

tential,  at  every  point  of  the  surface,  varies  as  a;2  =  F0  -5 

66 

suppose.  This  case  materially  differs  from  the  two  just  con 
sidered,  for  since  a?  does  not,  when  substituted  for  V,  satisfy 
the  equation  y2  V—  0,  the  potential  at  internal  points  cannot 
in  general  be  proportional  to  xz.  We  have  therefore  first  to 
investigate  a  function  of  x,  y,  z,  or  of  e,  v,  v  which  shall 
satisfy  the  equation  \72F=0,  shall  not  become  infinite  within 
the  surface  of  the  ellipsoid,  and  shall  be  equal  to  o?  on  its 
surface. 

Now  we  know  that,  generally 
(62  +  o>)  (c2  +  o>)  x*  +  (c2+  w)  (a2  +  a)  ya  +  («2  +  fl>)  (&2  +  a))  z* 

—  (a2  +  a))  (b2  +  a))  (c2  +  ft))  =  (e  —  a))  (u  —  CD)  (V  —  ft)). 

And,  if'^,  ^2  be  the  two  values  of  co  which  satisfy  the 
equation 

(62+ft))  (c2+ft))  +  (c2+ft))  (a2+co)  +  (a2+o>)  (62+ft))  =  0...(l), 

we  see  that 

V2(6-^)(i;-^)(u'-^)=0, 

and  ^(€-_^}-(v-^)(t/--^)-a 

And,  by  properly  determining  the  coefficients  AQ)  A^  A2, 
it  is  possible  to  make 


when  62cV  +  c2aV  +  aW  -  a26V  =  0. 

10—2 
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Hence,  the  expression  (2)  when  A0,  Alt  A2  are  properly 
determined  will  satisfy  all  the  necessary  conditions  for  an 
internal  potential,  and  will  therefore  be  the  potential  for 
every  internal  point. 

Now,  we  have  in  general 


-  (a"  +  0,)  (V  +  <y  (c2  +  tf  ,)  =  («  -  0J  (v  -  6}  (i/  -  fl.) 
(//  +  6}  (c2  +  02)  a?  +  (c2  +  02)  (a2  +  6}  f  +  (a2  +  02)  (62  +  ^  z* 

-  («2  +  «,)  (V  +  6)  (c2  +  0J  =  (e  -  0)  (v  -  0,  (v'  -  0,) 
and,  over  the  surface 

c?Vz*  -  aW  =  0. 


Hence,  S-  being  any  quantity  whatever,  we  have,  all  over 
the  surface, 

(If  +  ty  (c2  +  &)  a;2  +  (c2  +  *)  (a2  +  &)  ^  +  (a2 


and  therefore,  putting  ^-  =  —  a2, 


Hence,  the  right-hand  member  of  this  equation  possesses 
all  the  necessary  properties  of  an  internal  potential.  It 
satisfies  the  general  differential  equation  of  the  second  order, 
does  not  become  infinite  within  the  shell,  and  is  proportional 
to  x2  all  over  the  surface. 

"We  observe,  by  equation  (1),  that 
(V  +  o>)  (c2  +  to)  +  (c2+  a))  (a'-f  w)  +  (a2  +  to)  (62+a>)  =  3  ft-®)  ft-») 
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identically,  and  therefore,  writing  —  a2  for  o>, 

(a2  -  &2)  (a2  -  c2)  =  3  (a2  +  0X)  (a2  +  02). 

Hence,  over  the  surface  of  the  shell, 


and  we  therefore  have,  for  the  internal  potential, 

ejjv-ej  (I/-0J  .  (e-ej(v-oj(v-oj    \ 

• 


This  is  not  admissible  for  external  points,  as  it  becomes 
infinite  at  an  infinite  distance.  We  must  therefore  substi 
tute  for  the  factor  e  —  6l 

T      • 


*OQ 

(e  ~  ^  >•     - 


.  r 

'  Jo      - 


*+ 


with  a  similar  substitution  for  e-^2,  thus  giving,  for  the 
external  potential, 


FJCe-fl.xu-a.Kv-fl.)  r 
'  31  *,(V,)(««'+*,)    J« 

.  r 

'  Jo        - 


H«-g1)(u'-g1)  r 

,-*,)  («'+»,)    /.       - 


r          ^  4 

.«  [(f  +a!)(f  +J2)(^H-c2)J^  '  J 
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The  distribution  of  density  over  the  surface,  correspond 
ing  to  this  distribution  of  potential,  may  be  investigated  by 
means  of  the  formula 


or  its  equivalent  in  Art.  13  of  this  Chapter.   We  thus  find  that 
1    a     F0    r        fr-^fc'-fl,) 


.  r 
'  J.    - 


U 


30.  The  investigation  just  given,  of  the  potential  at  an 
external  point  of  a  distribution  of  matter  giving  rise  to  a 
potential  proportional  to  a?  all  over  the  surface,  has  an  in 
teresting  practical  application.  For  the  Earth  may  be  re 
garded  as  an  ellipsoid  of  equilibrium  (not  necessarily  with 
two  of  its  axes  equal)  under  the  action  of  the  mutual  gravi 
tation  of  its  parts  and  of  the  centrifugal  force.  If,  then, 
V  denote  the  potential  of  the  Earth  at  any  point  on  or  with 
out  its  surface,  and  fl  the  angular  velocity,  of  the  Earth's 
rotation,  we  have,  as  the  equation  of  its  surface,  regarded  as 
a  surface  of  equal  pressure, 


/.  V+  ^  fl2  (xz  +  ?/2)  =  a  constant,  IT  suppose. 

Hence,  if  a,  b,  c  denote  the  semi-axes  of  the  Earth,  we 
have,  for  the  determination  of  F,  the  following  conditions  : 
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5+%+S-  .............  •••••<»• 

V=  0  at  an  infinite  distance  ...............  (2), 

when 


The  term  II  will,  as  we  know,  give  rise  to  an  external 
potential  represented  by 

n  r        *+        +  r 

Je        +c+£2(+c2*  '  J0        +a 


The  two  terms  —-  HV,  —^H2?/2,  will  give  rise  to  terms 


which  may  be  deduced  from  the  value  of  F2  just  given  by 
successively  writing  for  F0,  -^HV,  and  —  H2&2,  and  (in 
the  latter  case)  putting  62  for  az  throughout.  We  thus  get 


P 

*'  J 


0 

6 


fl2/    a 


Jo.(&-0-}* 
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31.  Any  rational  integral  function  V  of  x,  y,  z,  which 
satisfies  the  equation  y2F=0,  can  be  expressed  in  a  series 
of  Ellipsoidal  Harmonics  of  the  degrees  0,  1,  2...i  in  x,  y,  z. 
For  if  V  be  of  the  degree  it  the  number  of  terms  in  V  will 
,  (*  +  !)(*  +  2)  (t  +  3)  AT 
be  ~  6~  -•  Now  the  condition  y*V=Q  is 

equivalent  to  the  condition  that  a  certain  function  of  a,  y,  z 
of  the  degree  i-2,  vanishes  identically,  and  this  gives  rise 

(t  —  1)  i  (*  +  l)         j.  . 
to  "       —  Q  —     —  conditions.     Hence  the  number  of  inde 

pendent  constants  in  V  is 

(t  +  1)  (i+2)(t+3)      (*-l)  *(/+!) 
6  ~Q~      ~~' 

or  (i+  I)2.     And  the  number  of  ellipsoidal  harmonics  of  the 

degrees  0,  1,  2...i  in  a?,  y,  z  or  of  the  degrees  0,  5,  1,  f  ...- 

L        2t     2 
in  e,  v,  v  ,  is,  as  shewn  in  Arts.  6  to  10  of  this  Chapter, 

1  +  3  +  5+..  .  +  2i  +  l, 

or  (i  +  1)2.     Hence  all  the  necessary  conditions  can  be  satis 
fied. 

32.  Again,  suppose  that  attracting  matter  is  distributed 
over  the  surface  of  an  ellipsoidal  shell  according  to  a  law  of 
density  expressed  by  any  rational  integral  function  of  the 
co-ordinates.  Let  the  dimensions  of  the  highest  term  in  this 
expression  be  i,  then  by  multiplying  every  term,  except  those 
of  the  dimensions  i  and  i  —  1  by  a  suitable  power  of 


we  shall  express  the  density  by  the  sum  of  two  rational  inte 
gral  functions  of  x,  y,  z  of  the  degrees  it  i  —  1,  respectively. 
The  number  of  terms  in  these  will  be 
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And  any  ellipsoidal  surface  harmonic  of  the  degree  it  i  —  2... 
in  x,  y,  z,  may,  by  suitably  introducing  the  factor 


be  expressed  as  a  homogeneous  function  of  x,  y,  &  of  the 
degree  i  ;  also  any  such  harmonics  of  the  degree  i—  1,  ^—3... 
in  x,  y,  z  may  be  similarly  expressed  as  a  homogeneous 
function  of  x,  y,  z  of  the  degree  i  —  1.  And  the  total  number 
of  these  expressions  will,  as  just  shewn,  be  (i+  I)2,  hence  by 
assigning  to  them  suitable  coefficients,  any  distribution  of 
density  according  to  a  rational  integral  function  of  x,  y,  z 
may  be  expressed  by  a  series  of  surface  ellipsoidal  harmonics, 
and  the  potential  at  any  internal  or  external  point  by  the 
corresponding  series  of  solid  ellipsoidal  harmonics. 

33.  Since  any  function  of  the  co-ordinates  of  a  point  on 
the  surface  of  a  sphere  may  be  expressed  by  means  of  a  series 
of  surface  spherical  harmonics,  we  may  anticipate  that  any 
function  of  the  elliptic  co-ordinates  v,  v  may  be  expressed  by 
a  series  of  surface  ellipsoidal  harmonics.  No  general  proof, 
however,  appears  yet  to  have  been  given  of  this  proposition. 
But,  assuming  such  a  development  to  be  possible  at  all,  it 
may  be  shewn,  by  the  aid  of  the  proposition  proved  in 
Art.  15  of  this  Chapter,  that  it  is  possible  in  only  one  way, 
in  exactly  the  same  way  as  the  corresponding  proposition 
for  a  spherical  surface  is  proved  in  Chap.  IV.  Art.  11. 

The  development  may  then  be  effected  as  follows.  De 
noting  the  several  surface  harmonics  of  the  degree  i  in  xy  y,  z, 

or  -  in  v,  v,  by  the  symbols   Vp\   F/2>,  .  .  .  F,<3i+1>,  and  by 

2 
F(v,  v)  the  expression  to  be  developed,  assume 

^   F(v,  v')=CQV0 


Then  multiplying  by  eV^  and  integrating  all  over  the 
surface,  we  have 


feF  (v,  v'}  F/*>  dS  =  CPJe  (F/*))1  dS. 
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The  values  ot  feF(v,  v'W>dS,  and  of  je  ( F/"))2  dS  must 

be  ascertained  by  introducing  the  rectangular  co-ordinates 
x,  y,  z,  or  in  any  other  way  which  may  be  suitable  for  the 
particular  case.  The  coefficients  denoted  by  C  are  thus 
determined,  and  the  development  effected. 


EXAMPLES. 

1.  Prove  that   (sin  ey  =  ^P0-~Pa  +  ^ P4. 

"Why  cannot  (sin  0)3  be  expanded  in  a  finite  series  of  spherical 
harmonics  ? 

0 
1          1          1  1  +  am « 

2.  Prove  that  1  +  ^  Px+  ^  P2  +  -  P3+  ...  =log J!  . 

sin  | 

3.  Establish  the  equations 

( w.2  —  1 )    ,  " 

ClfJi 

4.  If  /*  —  cos  Oj  prove  that  / 


and  also  that 


f  i  +  ra 

—— 


i  +  ra       /        fl  \ 

(_  1)'*-  ~——~  (cos2  1  )  +  ... 
(\m)2\i-m  \        2/ 


5.     Prove    that,  if    a    be  greater    than    c,   and   i   any   odd 
integer  greater  than  m, 


/i  /fjp  \  s 
(^~}d^  =  i(i 
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7.     Prove  that,  when  ^  =  *  1, 
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. 
a//."1       \i-rti    2m  I  771 


8.     Prove  that 


IP  P 

*l         *i      •••     *| 

P         P 
JT       JT 


P        P  P 

-1  i>     ^  i+i  •••     •«  jj 


is  a  numerical  multiple  of 


9.     Prove  the  following  equation,  giving   any   Laplace's   co 
efficient  in  terms  of  the  preceding  one : 


where  Cp  =  /x//  +  J  1  -  /x.2  Jl  -  ^  cos  (<o  -  to7)  and  C  is  zero  if  ?i  be 
even,  and 


(-  1)  2  -^HT 


if  w  be  odd. 


10.     If  ifjt  Jc  be  three  positive  integers  whose  sum  is  even, 
prove  that 


r 

Ji 


2A  ...  (k  +  i-j)          2.4  ...  (i+j-k) 


1.3  ...  (i  +  ^'  +  ^- 
Hence   deduce   the  expansion  of  P.P^  in,  a   series  of   zonal 
harmonics. 

11.  Express     x*y  +  y*  +  yz  +  y  +  z     as   a   sum   of  spherical 
harmonics. 

12.  Find  all  the  independent  symmetrical  complete  harmonics 
of  the  third  degree  and  of  the  fifth  negative  degree. 

13.'  Matter  is  distributed  in  an  indefinitely  thin  stratum  over 
the  surface  of  a  sphere  whose  radius  is  unity,  in  such  a  manner 
that  the  quantity  of  matter  laid  on  an  element  (8S)  of  the  surface 
is  S£  (  1  +  ax  +  ly  +  cz  +fx*  +  gy* 
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where  x,  y,  z  are  rectangular  co-ordinates  of  the  element  SS  re 
ferred  to  the  centre  as  origin,  and  a,  b,  c,  f,  g,  h  are  constants, 
Find  the  value  of  the  potential  at  any  point,  whether  internal 
or  external. 

14.     If  the  radius  of  a  sphere  "be  r,  and  its  law  of  density  be 
p  =  ax  +  1)y  +  CZ)  where  the  origin  is  at  the  centre,  prove  that  its 

potential  at  an  external  point  (£,  17,  £)  is  —  -     (a£  +  brj  +  c£)  where 


E  is  the  distance  of  (£,  17,  £)  from  the  origin. 

15.  Let  a  spherical  portion  of  an  infinite  quiescent  liquid  be 
separated  from  the  liquid  round  it  by  an  infinitely  thin  flexible 
membrane,  and  let  this  membrane  be  suddenly  set  in   motion, 
every  part  of  it  in  the  direction  of  the  radius  and  with  velocity 
equal  to  S0  a  harmonic  function  of  position  on  the  surface.     Find 
the  velocity  produced  at  any  external  or  internal  point  of  the 
liquid.      State   the  corresponding   proposition  in  the   theory   of 
Attraction. 

16.  Two  circular  rings  of  fine  wire,  whose  masses  are  M  and 
J/7,  and  radii  a  and  a',  are  placed  with  their  centres  at  distances 
£>,   b',  from  the  origin.     The  lines  joining   the   origin  with  the 
centres  are  perpendicular  to  the  planes  of  the  rings,  and  are  in 
clined  to  one  another  at  an  angle  0.     Shew  that  the  potential  of 
the  one  ring  on  the  other  is 


where  s- 


and  B'n  and  Qn  are  the  same  functions  of  V  and  a'  and  of  cos  0 
and  sin  0  respectively,  and  c  is  the  greater  of  the  two  quantities 
\/c*2  +  W  and  *Ja*  +  b'\ 

17.  A  uniform  circular  wire,  of  radius  a,  charged  with 
electricity  of  line-density  e,  surrounds  an  uninsulated  concentric 
spherical  conductor  of  radius  c;  prove  that  the  electrical  density 
at  any  point  of  the  surface  of  the  conductor  is 


l     5 

-1'  5- 


the  pole  of  the  plane  of  the  wire  being  the  pole  of  the  harmonics. 
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18.  Of  two  spherical  conductors,  one  entirely  surrounds  the 
other.      The   inner  has   a    given  potential,    the  outer  is  at  the 
potential  zero.     The   distance    between    their    centres   being  so 
small  that   its    square   may  be  neglected,  shew  how  to  fincf  the 
potential  at  any  point  between  the  spheres. 

19.  If  the   equation   of  the    bounding  surface  of  a   homo 
geneous  spheroid  of  ellipticity  e  be  of  the  form 


prove  that  the  potential  at  any  external  point  will  be 

~-~^P*> 

where  C  and  A  are  the  equatoreal  and  polar  moments  of  inertia 
of  the  body. 

Hence  prove  that  F  will  have  the  same  value  if  the  spheroid 
be  heterogeneous,  the  surfaces  of  equal  density  differing  from 
spheres  by  a  harmonic  of  the  second  order. 

20.  The  equation  R  =  a(l+ay)  is  that  of  the  bounding 
surface  of  a  homogeneous  body,  density  unity,  differing  slightly 
in  form  and  magnitude  from  a  sphere  of  radius  a;  a  is  a 
small  quantity  the  powers  of  which  above  the  second  may  be 
neglected;  and  y  is  a  function  of  two  co-ordinate  angles,  such 
that 


where  Y0,  Ft ...  ZQ,  Zl  ...  are  Laplace's  'functions.  Prove  that 
the  potential  of  the  body's  attraction  on  an  external  particle, 
the  distance  of  which  from  the  origin  of  co-ordinates  is  r,  is 
given  by  the  equation 
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21.  If  if  be  the  mass  of  a  uniform  hemispherical  shell  of 
radius  c,  prove  that  its  potential,  at  any  point  distant  r  from 
the  centre,  will  be 


__ 

l      2.4     3c3 


2.4.6V     2,4.6.8    7c5 


M     I 
°r        + 


2.4.6       r6     2.4.6.8        r8 

according  as  r  is  less  or  greater  than  c;  the  vertex  of  the  hemi 
sphere  being  at  the  point  at  which  p  =  1. 

22.  A.  solid  is  bounded  by  the  plane  of  xy,  and  extends  to 
infinity  in  all  directions  on  the  positive  side  of  that  plane. 
Every  point  within  the  circle  x*  +  y2  =  a2,  z~0  is  maintained  at 
the  uniform  temperature  unity,  and  every  point  of  the  plane  xy 
without  this  circle  at  the  uniform  temperature  0.  Prove  that, 
when  the  temperature  of  the  solid  has  become  permanent,  its 
value  at  a  point  distant  r  from  the  origin,  and  the  line  joining 
which  to  the  origin  is  inclined  at  an  angle  6  to  the  axis  of  z  will 
be 

r      I      r*     1.3       r5 

0  i        ~*"T7-t3     <i~~  o      T  •*  I     5  >    '  '  ' 

1  a     2      a?     2.4       a 

1.3...(2»-1)         £?'• 
^     '       2.4...2i        2i+1  a*i+1 
if  r  <  a,  and 


J^l^p*.  ,    iy  1.3..  .(2^-1) 

r*     2.4     3r4          ~(      >       2.4...2i 


if  r  > 


23.  Prove  that  the  potential  of  a  circular  ring  of  radius  c, 
whose  density  at  any  point  is  cos  m\j/,  cij/  being  the  distance  of  the 
point  measured  along  the  ring  from  some  fixed  point,  is 
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2n-  cos  m&  (sin  OY  {  z — -r- « — ~ -y— -  - 

<    (2.4.  6...2m  c//xw  r 


d 


+  .. 


wliere  r  is  greater  than  c.     If  r  be  less  than  c,  r  and  c  must  be 
interchanged. 

24.  A  solid  is  bounded  by  two  confocal  ellipsoidal  surfaces,  and 
its  density  at  any  point  P  varies  as  the  square  on  the  perpendicular 
from  the  centre  on  the  tangent  plane  to  the  confocal  ellipsoid 
passing  through  P.  Prove  that  the  resultant  attraction  of  such 
a  solid  on  any  point  external  to  it  or  forming  a  part  of  its  mass 
is  in  the  direction  of  the  normal  to  the  confocal  ellipsoid  passing 
through  that  point,  and  that  the  solid  exercises  no  attraction  on  a 
point  within  its  inner  surface. 
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